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As part of the European LAGUNA design study on a next-generation neutrino detector, we propose 
the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a multipurpose neu- 
trino observatory. The outstanding successes of the Borexino and KamLAND experiments demon- 
strate the large potential of liquid-scintillator detectors in low-energy neutrino physics. Low energy 
threshold, good energy resolution and efficient background discrimination are inherent to the liquid- 
scintillator technique. A target mass of 50 kt will offer a substantial increase in detection sensitivity. 

At low energies, the variety of detection channels available in liquid scintillator will allow for an 
energy- and flavor-resolved analysis of the neutrino burst emitted by a galactic Supernova. Due 
to target mass and background conditions, LENA will also be sensitive to the faint signal of the 
Diffuse Supernova Neutrino Background. Solar metallicity, time-variation in the solar neutrino flux 
and deviations from MSW-LMA survival probabilities can be investigated based on unprecedented 
statistics. Low background conditions allow to search for dark matter by observing rare annihilation 
neutrinos. The large number of events expected for geoneutrinos will give valuable information 
on the abundances of Uranium and Thorium and their relative ratio in the Earth's crust and 
mantle. Reactor neutrinos enable a high-precision measurement of solar mixing parameters. A 
strong radioactive or pion decay-at-rest neutrino source can be placed close to the detector to 
investigate neutrino oscillations for short distances and sub-MeV to MeV energies. 

At high energies, LENA will provide a new lifetime limit for the SUSY-favored proton decay mode 
into kaon and antineutrino, surpassing current experimental limits by about one order of magnitude. 
Recent studies have demonstrated that a reconstruction of momentum and energy of GeV particles 
is well feasible in liquid scintillator. Monte Carlo studies on the reconstruction of the complex event 
topologies found for neutrino interactions at multi-GeV energies have shown promising results. If 
this is confirmed, LENA might serve as far detector in a long-baseline neutrino oscillation experiment 
currently investigated in LAGUNA-LBNO. 
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1 Introduction 

Over the past decades, neutrinos have been firmly es- 
tablished as astronomical messengers. The feeble in- 
teraction strength of these elusive particles requires 
unusually large detectors, but on the other hand al- 
lows us to investigate processes in the deep interior of 
stars that are shrouded from view in other forms of 
radiation. Neutrino astronomy therefore complements 
observations in the electromagnetic spectrum, charged 
cosmic rays, and gravitational waves. In recognition 
of this importance, the pioneering first observations of 
solar and supernova (SN) neutrinos were honored with 
the Physics Nobel Prize in 2002. 

Even the first solar neutrino observations about 40 
years ago showed an apparent deficit that today is un- 
ambiguously explained by flavor oscillations. In this 
way neutrino astronomy has triggered an avalanche 
of fundamental discoveries, shedding completely new 
light on the inner properties of neutrinos with ramifi- 
cations both for the fundamental theory of elementary 
particles and the universe at large. 

With the standard three-flavor oscillation scenario 
established, neutrinos can be used as new messengers 
from astrophysical sources. In this context, a variety 
of far-reaching questions can be addressed, notably 

• Is the core-collapse SN paradigm correct? Are 
there substructures in the neutrino signal? 

• How large is the flux of the diffuse SN neutrino 
background (DSNB) and what is its spectrum? 

• What is the Sun's metal content? Are there 
time-variations in the solar fusion rate? 

• How large is the concentration of radioactive el- 
ements in the Earth and what is their contribu- 
tion to its heat flow? 

• What is the dark matter of the universe? 

On the other hand, neutrinos remain fundamental 
particle-physics messengers and large-scale detectors 
will shed new light on topics like 

• What is the value of #13? 

• Is the CP symmetry violated among leptons? 

• What is the neutrino mass hierarchy? 

• Do sterile neutrinos exist? 

• Are there non-standard neutrino interactions 
and how do they affect flavor oscillations? 

• Is baryon number conserved? 

These question will be addressed by observing low- 
energy neutrinos from SNe, the Sun, Earth, reactors, 
and radioactive sources, by neutrino beams and at- 
mospheric neutrinos in the GeV range, and finally by 



signatures of possible nucleon decays in the detector 
material itself. 

The small event rate of neutrino interactions or the 
search for extremely rare processes requires a large 
target mass. In the past, large-volume unsegmented 
detectors have played a dominant role in this field. 
Originally triggered by the search for nucleon decay, 
the Kamiokande and later Super-Kamiokande water 
Cherenkov detectors provided crucial measurements 
of solar, atmospheric, SN and beam neutrinos. It 
was only the huge target mass of 50 kt that allowed 
Super-Kamiokande to accrue enough statistics to mea- 
sure precisely the deformation of atmospheric neutrino 
spectra caused by flavor oscillations. In a parallel de- 
velopment, liquid-scintillator detectors on the kiloton 
scale explored neutrino fluxes at energies below 5 MeV. 
In particular, the KamLAND measurements of reactor 
neutrino oscillations tightly confine the mass-squared 
difference of solar neutrino mixing, while Borexino 
confirmed solar neutrino oscillations at sub-MeV en- 
ergies. Both detectors provided first evidence for the 
faint geoneutrino signal originating from radioactive 
elements embedded in the Earth's crust. 




FIG. 1: Artist's view of the LENA detector: The detector 
tank is 100 m in height and 30 m in diameter. See Fig. [2] 
for details. 
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Based on this success, we propose a next-generation 
neutrino observatory LENA ( "low energy neutrino as- 
tronomy"). It is foreseen as an unsegmented liquid- 
scintillator detector of 50 kt target mass (Fig. [I]), com- 
bining the advantages of the low-energy threshold and 
background discrimination capabilities of Borexino 
and KamLAND with the size of Super-Kamiokandc. 

LENA will profit from the virtues of the scintilla- 
tor technique that were impressively demonstrated by 
KamLAND and Borexino. 

• Good energy resolution below 10 MeV. 

The light yield is at least 200 photoelectrons per 
MeV, corresponding to about 3% energy resolu- 
tion at 5 MeV. 

• Low detection threshold. A neutrino energy 
of 1.8 MeV is the threshold for inverse beta de- 
cay. For electron scattering, the threshold can 
be at a recoil energy as low as 200 kcV, the limit 
arising from the intrinsic background of radioac- 
tive 14 C in the scintillator. 

• Excellent background discrimination. The 

final-state neutron of inverse beta decay pro- 
vides a clear coincidence signature for v e de- 
tection. Pulse shape analysis allows for an ef- 
ficient discrimination against fast neutrons and, 
for detecting elastic neutrino-electron scattering, 
against alphas and even positrons. 

• Radio purity. Years of development and expe- 
rience in Borexino have advanced the techniques 
for scintillator purification, identifying the most 
efficient methods. 

• Self shielding. A large monolithic detector 
shields its central detection volume against ex- 
ternal backgrounds. 

In most respects, the performance is competitive with 
a water Cherenkov detector of several times its size. 

LENA will be a true multi-purpose facility. High- 
statistics measurements of strong neutrino sources like 
a galactic core-collapse SN, the Sun or the Earth's in- 
terior will resolve energy spectra and their time evolu- 
tion in unprecedented detail. Reactor neutrinos enable 
a high-precision measurement of the "solar" neutrino 
mixing parameters. In addition, a strong radioactive 
neutrino source can be placed close to the detector 
to investigate flavor oscillations at short distances and 
sub-MeV energies. At the same time, the search for 
very rare events becomes possible because the excel- 
lent background rejection allows to identify a hand- 
ful of events out of several years of data. Thus, the 
faint flux of the predicted Diffuse Supernova Neutrino 
Background (DSNB) is well within reach. Likewise, 
observation of rare annihilation neutrinos allows for 
indirect dark matter search. 



This rich low-energy program is complemented by 
several physics objectives at GeV energies. LENA 
will further advance the search for proton decay and 
thus baryon number violation. The new lifetime sen- 
sitivity for the proton decay mode into kaon and an- 
tineutrino, favored by supersymmetric theories, will 
surpass current experimental limits by about one or- 
der of magnitude. Moreover, recent studies indicate 
that a large- volume liquid-scintillator detector can re- 
solve both momentum and energy of GeV particles 
with a precision of a few percent. Monte Carlo sim- 
ulations of the complex event topologies of charged- 
current neutrino interactions show promising results 
for the reconstruction capabilities. These techniques 
may offer the opportunity to use LENA as far detec- 
tor in long-baseline neutrino oscillation experiments, 
either for an accelerator-produced neutrino beam or 
atmospheric neutrinos. 

LENA is one of three options discussed within 
the LAGUNA and the forthcoming LAGUNO-LBNO 
design studies that are sponsored by the European 
Union under the 7th Framework Programme. This 
design study aims at the eventual construction of 
a large-volume neutrino observatory in a European 
underground laboratory based on scintillator, water 
Cherenkov (MEMPHYS), or liquid argon (GLACIER) 
techniques. Due to the high level of expertise built up 
in several European and international research groups 
and dedicated R&D activities over the past years, the 
liquid-scintillator technique can be regarded as suffi- 
ciently mature to allow for an early start of detector 
realization. Based on recent feasibility studies, the 
LENA construction time is estimated to about eight 
years. 

This paper lays out the science case for LENA and 
the current state of R&D and detector design activi- 
ties. It is meant to support and justify a proposal for 
the construction of a next-generation large neutrino 
observatory based on the liquid-scintillator technique. 
It provides a work of reference for future discussions 
and decision making. The paper starts out with the 
state of R&D and detector design activities in Sec.[2j 
followed by an outline of the technical detector prop- 
erties (Sec. [3]). The core physics objectives are in the 
low-energy domain where a liquid-scintillator detector 
can play out its unique capabilities particularly well. 
The main low-energy topics revolve around solar, su- 
pernova, reactor and geo neutrinos that are described 
in Sec. [3] In addition, LENA has convincing capabil- 
ities in the range of GeV energies, where the search 
for nucleon decay and flavor oscillation physics with 
accelerator-produced beams and atmospheric neutri- 
nos form the main topics that are discussed in Sec.[6j 
The paper concludes with a brief summary and out- 
look in Sec.0 







2 Detector design 



Design, construction, and operation of the LENA de- 
tector will be a challenging endeavor. However, there 
are two neutrino detectors in operation that already 
anticipate scale and techniques of the LENA project: 
The Super-Kamiokande detector is of almost the same 
volume, featuring similar requirements concerning de- 
tector cavern, photocoverage and number of channels. 
On the other hand, the enormous amount of R&D that 
led to the tremendous success of the Borexino experi- 
ment can be re-applied for LENA, covering questions 
concerning the liquid scintillator, the radiopurity of 
the used materials and their purification, requirements 
for photosensors and electronic read-out and so forth. 
Based on this foundation, but also on the laboratory 
and design activity carried out in the last few years es- 
pecially for LENA, the following section describes the 
current design draft for LENA. 

Fig.[2]shows a schematic overview of the current LENA 
design: 

Cavern 



height: 115 m, diameter: 50 m 
shielding from cosmic rays: -4,000 m.w.e 

Muon Veto 



plastic scintillator panels (on top) 
Water Cherenkov Detector 
3,000 phototubes 
100 kt of water 
reduction of fast 
neutron background 

Steel Cylinder 

height: 100 m, diameter: 30 m 
70 kt of organic liquid 
30,000 - 50,000 phototubes 

Buffer 



thickness: 2 m 

non-scintillating organic liquid 
shielding from external radioactivity 



Nylon Vessel 

separating buffer liquid 
and liquid scintillator 



Target Volume 




height: 100 m, diameter: 26 m 
50 kt of liquid scintillator 

FIG. 2: Schematical view of the LENA detector [T]. 

Laboratory. The detector will be constructed in a 
dedicated cavern, about 115 m in height. The shape 
will depend on the laboratory site: Pyhasalmi and 
Frejus will be described as exemplary sites in Sec. |2.1| 
The aspired rock shielding above the detector corre- 
sponds to 4000mwe, a requirement fulfilled by both 
sites. 

Tank. The liquid-scintillator will be contained in a 
cylindric steel or concrete tank of 100 m height and 
30 m diameter. Several design options are discussed 
in Sec.|2.2| Inside the tank, the volume is divided by 



a thin nylon vessel into the buffer volume shielding 
external radioactivity and the target volume. 

Liquid scintillator. The target volume is 26 m 
in diameter and 100 m in height, corresponding to 
5.3xl0 4 m 3 . Depending on the exact composition of 



the liquid scintillator (Sec 2.3), the target mass ranges 
from 45 to 53 kt. The buffer volume is filled with an 
inactive liquid, which should have a similar density as 
the scintillator in order to minimize buoyancy forces 
on the nylon vessel. The scintillator transparency will 
be monitored by a system of reference LEDs viewed 
by PMTs across the tank diameter. Light yield cal- 
ibration at low energies will be performed by intro- 
ducing radioactive sources at different positions inside 
the detection volume, while cosmic muons tracks re- 
constructed by the top muon veto provide information 
on the high-energy detector response. 

Nylon vessel. A cylindrical nylon tube of 26 m di- 
ameter and a surface of ~10 4 m 2 is needed to separate 
the active target from the inactive buffer volume. The 
vessel will be very thin (^100 urn), and fastened by 
two flanges to the tank lid and bottom. Intermediate 
support by nylon strings reaching out to the tank walls 
will help to carry the weight and to keep the cylindri- 
cal shape. Compared to Borexino, the construction of 
the vessel might prove easier as - compared to a sphe- 
ciral vessel - a cylinder needs much less glued seams. 
Nevertheless, clean and radiopure fabrication of a ves- 
sel of these dimensions will be challenging. Currently, 
studies are on-going to integrate a small buffer vol- 
ume within the encapsulation of each individual PMT 
which might allow to omit the vessel, filling the whole 
tank with active scintillator. 

Photomultipliers. The intended photosensitive cov- 
erage is 30 % of the inner tanks walls. This requires 
e. g. ~45 000 eight-inch photomultiplier tubes (PMTs). 
Currently PMTs with a photocathode diameter be- 
tween 5 and 10 inch are the most likely solution. Re- 
flective light-concentrators mounted on the PMTs will 
be used to reduce the number of PMTs (Sec 2.4). 



PMT response and timing can be calibrated by a sys- 
tem of optical fibers fed by ns-long LASER pulses. 

Readout electronics. A further option to reduce the 
large number of channels is to group several PMTs into 
a PMT array, featuring a common high voltage sup- 
ply, signal digitization and readout channel. Possible 
solutions for the read-out electronics are discussed in 
Secl2~5l 



Muon veto. Cosmic muons crossing the main detec- 
tor will be identified by layers of plastic scintillator 
panels, Resistive Plate Channels (RPCs), or limited 
streamer tubes mounted above the upper lid of the 
detector. A dense instrumentation featuring several 
layers would offer the possibility to aid the reconstruc- 
tion of muon tracks in the scintillator. On the outside 



7 



of the tank, the interspace to the cavern walls is filled 
with water (at least 2 m in width) shielding the in- 
ner detector from external radiation coming from the 
rock and from muon-induced neutrons. The outer tank 
walls can be equipped with PMTs to identify cosmic 
muons passing the detector by their Cherenkov light. 

2.1 Laboratory sites 

No underground laboratory existing today is suffi- 
ciently large to host the LENA detector. This im- 
plies that a cavern of appropriate size must be ex- 
cavated, along with additional shafts and tunnels to 
house the auxiliary systems for filling and operation 
of the detector. To minimize the associated costs, it 
seems reasonable to construct the detector adjacent to 
an underground infrastructure already existent, either 
an underground science laboratory or a deep mine. In 
the following, two exemplary sites in Europe are pre- 
sented that would suit the depth and infrastructure 
requirements of LENA. These places have been iden- 
tified in the course of the FP7 LAGUNA design study 
which will publish its conclusive results in 2011. 

2.1.1 Pyhasalmi 

The Pyhasalmi mine is located close to the geographic 
center of Finland, near the town of Pyhajarvi. The dis- 
tance from CERN is 2288 km. The mine is the deepest 
in Europe, the bottom level at ^1450 m. The prod- 
ucts are copper, zinc and pyrite, and operation will 
last at least until 2018. The mine already hosts a 
small underground laboratory, the Finnish Center for 
Underground Physics in Pyhasalmi (CUPP). The fea- 
sibility study for LENA at Pyhasalmi was carried out 
by the Finnish company Rockplan Ltd. 1 [2J. 

Geology. The characteristics of the rock surrounding 
the mine are well known due to the exploratory work 
performed by the mining company. The cavern will 
be constructed adjacent to the deepest level of the 
mine, about 500 m from the central mine shaft. At 
this depth, the rock will be very hard, dry and at 
23° C. The seismic activity in the region is very low. 

Background levels. The air content of radon 
at the deepest level is 20Bq/m 3 , the muon flux is 
l.lx 10 -4 /m 2 s. The closest nuclear power plant is 
350 km from the mine. The expected reactor v e 
background has been calculated to 1.9xl0 5 /cm 2 s 0. 
However, Finland plans to construct two additional 
reactors within the next decades. 

Excavation. The excavation will begin from the 
deepest mine level, creating two connecting tunnels 



1 Kalliosuunnittelu Oy Rockplan Ltd, Asemamiehenkatu 2, 
00520 Helsinki (Finland). 



from mine to cavern, a vertical shaft close to the lab- 
oratory, and the detector cavern itself (Fig. [3]). To ac- 
commodate the high vertical and horizontal stresses 
acting on the final cavern, its walls will be curved and 
the ceiling domed, as depicted in Fig. [2] Moreover, 
the plan view will be elliptical, the semi-major axis 
aligned to the direction of main stress. Overall, a vol- 
ume of 200 000 m 3 will be excavated, leaving room for 
an extensive water buffer. 

Infrastructure. The laboratory can profit from the 
already available underground infrastructure of the 
working mine (power, ventilation, transport). In addi- 
tion to the main shaft, a road tunnel spiraling from the 
surface to the deepest level of the mine will allow to 
bring large building elements to the detector cavern. 
There exists also the possibility to share to a certain 
extent the equipment and machines for underground 
excavation with the mining company. The transport 
of liquid scintillator to the laboratory to the mine will 
be possible both by road truck and freight trains as 
the mine is directly connected to the Finnish railway 
network. 




FIG. 3: LENA at Pyhasalmi (artistic impression by Rock- 
plan Ltd.). 



2.1.2 Frejus 

The Laboratoire Souterrain de Modane (LSM) is lo- 
cated adjacent to the Frejus road tunnel in the French- 
Italian Alps, connecting Modane (F) and Bardonec- 
chia (I). Originally, a new laboratory nearby has been 
discussed in the context of the MEMPHYS detector. 
However, the FP7 LAGUNA design study has shown 
that Frejus will suit well the requirements of LENA. 
Lately, a laboratory hosting both detectors in a com- 
mon infrastructure has been discussed (Fig. [4]). The 
feasibility study was carried out by the Suisse com- 
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FIG. 4: LENA and MEMPHYS at Frejus (artistic impression by Lombardi Ltd.). 



pany Lombardi Ltd. 2 [I]. 

Geology. The characteristics of the rock surrounding 
the road tunnels have been investigated thoroughly 
during its excavation in the 1970s. In spite of its 
ductile behaviour, the calc-schist formation is of good 
quality for building, relatively dry and at a tempera- 
ture of 30° C. Seismic activity is present but not dan- 
gerous. 

Background levels. The air content of radon 
was measured to 15Bq/m 3 in the LSM. Due to the 
large rock overburden of the Frejus mountain, cor- 
responding to 4 800 mwe, the muon flux is very low, 
5xl0 -5 /m 2 s. However, Frejus is close to the nuclear 
power plants of France, the closest at Bugey is merely 
130 km from the laboratory. The expected reactor v e 
background has been calculated to 1.6xl0 6 /cm 2 s [3]. 

Excavation. The excavation of the large detector 
caverns will be made in various stages, using a prelim- 
inary support of anchors and shotcrete. Once exca- 
vated, the cavern walls will be sealed by a strong layer 
of concrete, more than 1 m in width, in order to com- 
pensate for plasticity of the rock. The LENA cavern 
will be cylindrical with vertical walls, corresponding to 
an excavation volume of 100 000 m 3 . Two additional 
caverns would hold the MEMPHYS detector (Fig.[l]). 

Infrastructure. The laboratory can profit from the 
already available underground infrastructure of the 
road tunnel (ventilation). Currently, a safety tunnel is 
being excavated close to the already existing road tun- 
nel. Excavation works and transport of materials will 
mainly use this safety tunnel to minimize interference 



2 Lombardi SA Engineering Limited, Via R. Simen 19, 6648 
Minusio (Switzerland). 



with road traffic. Liquid scintillator will be supplied 
by road trucks. 

2.2 Detector tank 

The Rockplan prefeasibility study on the LENA de- 
tector tank resulted in four options, two out of steel 
and two out of concrete [5] . A concrete tank will be 
much more resistive to the compression generated by 
the water-scintillator density difference. However, it 
is also significantly more radioactive. To obtain the 
same fiducial volume for low-energy neutrinos, the di- 
ameter of the tank would have to be increased by 1- 
2 m. The cost saving due to the low price of concrete 
will roughly compensate the additional expenses for 
organic solvent. 

Conventional Steel Tank. A conventional tank re- 
quires a sizable amount of steel (driving the costs), 
and consists of many structural elements that would 
have to be brought separately into the laboratory and 
to be joined during underground construction. A ba- 
sic steel mainframe will be erected first, followed by 
load-bearing plates and a final stainless steel sheet 
welded on (Fig.[5| . There is also only one load-bearing 
layer separating scintillator and water. However, such 
a tank could be built straightforward and will be ro- 
bust. 

Sandwich Steel Tank. This tank will consist of 
thin walled sandwich elements, featuring a very high 
strength-to-weight ratio and providing a multiple-layer 
defense against liquid leaks (Fig. [6]). The elements can 
be prefabricated, reducing costs and allowing for ex- 
tensive quality control. They will be lifted into place 
and welded together. There is also the opportunity to 
equip the interior of the elements with thermal insula- 
tion or cooling pipes, or to use it for active leak preven- 
tion. However, this tank will require a lot of welding 
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liquid scintillators as active material and their inter- 



play with the detector hardware. Finally in Sec. 2.3.3 



FIG. 5: Conventional steel tank (Rockplan Ltd. 




FIG. 6: Sandwich steel tank (Rockplan Ltd.) 



(bearing the risk of radioactive contamination), and 
the mechanical design for tangential pressure will be 
challenging. 

Sandwich Concrete Tank. To assure water tight- 
ness, the concrete tank will be a steel-concrete-steel 
plate sandwich, about 30 cm in width. The external 
steel plates are connected to each other with welded re- 
bar. In construction, rings of steel plates will be lifted 
in place, the concrete being cast in between. Due to 
the slow curing of the concrete, construction will take 
a long time. Finally, an additional thin stainless steel 
sheet will be laser-welded on the inside for purity. 

Hollow Core Concrete Tank. Based on the Sand- 
wich Concrete Tank, hollow cores are now added to 
the concrete layer of the tank. This increases mechan- 
ical strength, and allows to install a cooling system or 
active leak prevention. However, this option is up to 
now little used in tank construction. 



2.3 Liquid scintillator 

The power of the liquid scintillator technology has 
been demonstrated in the past by successful neutrino 
experiments like Borexino [6] or KamL AND [7] . Large 
target masses, high energy resolution and a low en- 
ergy threshold are beneficial characteristics that en- 
able real-time detection of rare low-energy events. As 
liquid scintillator is the central component of the de- 
tector, this chapter will cover the main properties of 



the most promising scintillator mixtures are presented. 
At present, LAB as solvent with the admixture of PPO 
and Bis-MSB as solutes is favored. 

2.3.1 Scintillator properties 

Light output and quenching. Organic scintillators 
are excited by charged particle radiation or ultraviolet 
(UV) light. In the molecular dccxcitation process, UV 
light is emitted. Charged particles which cross a scin- 
tillating medium ionize and excite molecules on their 
track. However, ionization and radiationless deexcita- 
tion processes lead to a loss of fluorescence efficiency. 
In general, processes reducing the efficiency of energy 
to light conversion are known as quenching. 
Ionization and excitation densities are high for large 
energy deposition per unit length, which is the case 
for heavy particles, such as protons or as. This affects 
not only fluorescence efficiency, but also the scintilla- 
tion pulse shapes and can thus be used for particle 
identification. 

Emission spectra. The emission spectra of a single- 
component scintillator has a significant overlap with 
its own absorption spectra. This results in multiple ab- 
sorption and reemission processes where an important 
part of the information gets lost. In order to prevent 
additional losses in the energy conversion efficiency, 
usually one or multiple organic solutes are added. 
The solvent transfers its excitation energy mainly non- 
radiatively by dipole-dipole interaction to the solute 
(also called wavelength shifter or fluor) emitting a 
higher wavelength region (usually blue light) at which 
the solvent is transparent. 

Scintillation pulse shape. For excited states of the 
scintillator molecules, there are several processes to 
decay: photon emission, radiationless electronic relax- 
ation, inter-system crossing processes (i.e. transition 
between singlet and triplet states), and energy 
transfer by collision to other molecules. As several 
deexcitation modes are possible, scintillating pulse 
shapes commonly show more than one radiative decay 
constant. The shape of the scintillation pulse can be 
described by the sum of several exponential functions: 



(2.1) 



In this case, r t denotes the decay constants and rij 
the amplitudes of the corresponding decay processes. 
These constants are typical parameters for each scin- 
tillator material; for most organic scintillator the fast 
decay component dominates the emission. 

The amplitudes of the time components depend on 
the energy deposition per unit length. Consequently, 
the pulse shape can be used for particle discrimina- 
tion of a-particles or neutron-induced proton recoils 
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from electron signals and thus provides a fundamental 
method for background rejection. 

Attenuation length. As the scintillation photons 
propagate through the medium, absorption and scat- 
tering processes can occur. These processes strongly 
depend on the emitted wavelength; in general, the 
transparency of the medium increases with the wave- 
length. The main parameters for the description of the 
light propagation are the absorption length, the scat- 
tering length, and the self-absorption length of the so- 
lute. For short wavelengths, the solute self-absorption 
dominates. For large-volume particle detectors, long 
absorption and scattering lengths are required. As ab- 
sorption processes decrease the total number of pho- 
tons which arrive at the photo-sensors, the effective 
light yield of the detector is reduced. However, a high 
effective light yield is desirable, since it is directly con- 
nected to the energy resolution and energy threshold 
of the detector. 

Scattering length. Scattering processes change the 
direction of the scintillation photons, thereby elongat- 
ing the photon path from the scintillation event to 
the photosensors. Due to the increase in photon time 
of flight and the corresponding smearing of the photon 
arrival time patterns, scattering has a potentially dete- 
riorating effect on timing and pulse shape discrimina- 
tion. Laboratory studies have shown that the scatter- 
ing lengths of the investigated solvents are of the order 
of 20 to 30 m and therefore of the same order as the 
envisaged detector diameter of LENA [8] . The domi- 
nating processes are Rayleigh scattering off the solvent 
molecules and absorption-reemission processes on or- 
ganic impurities. The application of Winston cones 
reduces the fraction of scattered photons registered by 
the PMTs (SecpO]) . 



While MC studies for low event energies show that 
scattering has to be taken into account for correct ver- 
tex reconstruction and pulse shape discrimination, it 
has only a subdominant effect on the time resolution of 
the detector (Sec. 3.2 1. The latter is mainly governed 
by the fast fluorescence decay constant of the initial 
scintillation light [5]. The situation at GeV energies 
is comparable: As discussed in Sec. |6.1| scattering has 
only a small effect on proton decay sensitivity [TU] . 
while it has to be included in track reconstruction al- 
gorithms pO to obtain correct results. 



Radiopurity. As solvents of organic liquid scintil- 
lators are hydrocarbons, intrinsic radioactivity of the 
scintillator originates mainly from the 14 C f3 decay. 
The 14 C (3 background rate by far surpasses all neu- 
trino signals at energies below the endpoint of the /3 
spectrum at 156 keV. While the original solvent pro- 
duced in distillation plants is usually rather pure at the 
refinery, the surface contamination of the transporta- 
tion and experimental containers may dissolve in the 
liquid scintillator. However, experiments like Borexino 
have demonstrated the feasibility of ultrahigh radiop- 
urity levels in liquid scintillators, with contamination 



levels of 238 U at the order of 10 _17 g/g. The achieved 
concentrations of 232 Th and 40 K are even at the level 
of 10~ 18 g/g. The scintillator purification for LENA 
aims at the same level of radiopurity. 

2.3.2 Influence on detector design 

The properties of a liquid scintillator mixture di- 
rectly sets constraints on the technical design of the 
detector and vice versa. In order to optimize the detec- 
tor performance, one has to take a look at the impact 
of the scintillator properties on the detector geometry, 
its demands on the photosensors, on health and on 
handling issues. 

Geometry. The absorption length of the liquid scin- 
tillator is the main defining parameter for the detector 
geometry. It strongly affects the effective light yield of 
the scintillator detector, which accounts for the energy 
resolution and the energy threshold. With the most 
foreseen liquid scintillator compounds having absorp- 
tion lengths of 10-20 m (see Sec 2.3.31, it is unfeasible 
to build an unsegmented spherical 50kton detector. 
Still, the preferred geometry of the detector is a cylin- 
drical shaped tank of 15 m radius, which ensures a 
good effective light yield and volume to surface ratio. 
With a liquid-scintillator density close to 1 g/cm 3 , the 
corresponding height of the cylinder is about 100 m in 
case of a 50kton detector. 

Photosensors. The photosensors are the link be- 
tween scintillation light and data acquisition. There- 
fore, their properties play a crucial role for the perfor- 
mance of the detector. Especially the spectral sensi- 
tivity of the photosensors should match the emission 
spectra of the liquid scintillator, or at least be max- 
imal between 320-450 nm. For conventional bialkali 
photomultiplier tubes this wavelength range is quite 
common (Sect. 2.4). 



Liquid handling. Liquid handling comprises mainly 
the filling, pumping and storing of the liquid scintil- 
lator starting from its production. In this context 
radiopurity, purification, chemical compatibility, and 
safety issues have to be considered. 

The purification of the scintillator is necessary for 
reaching a low-level radioactive contamination. For 
example, this can be done by water extraction and 
prevention of exposing the scintillator to cosmic rays, 
i.e. fast transport on surface and underground storage. 
For such a large detector as LENA, it seems most feasi- 
ble to put purification plants on surface. At all times, 
the detector has to be kept away from oxygen con- 
tamination, which would induce a degradation of the 
scintillator. Thus, ultra-clean nitrogen gas must be 
used to flush pipes and tanks. Moreover, any material 
that comes into contact with liquid scintillator has to 
be tested on its chemical compatibility. Possible ma- 
terials are teflon or passivated stainless steel. 

Chemically, the solvent of an organic liquid scintil- 
lator is a hydrocarbon. Thus, flammability is the most 
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important concern. The Hazardous Materials Identifi- 
cation System (HMIS) classifies the danger in handling 
the liquid (from - safe to 4 - dangerous) concerning 
flammability, reactivity and health. The HMIS ratings 
for scintillators in flammability range from 1-3, reac- 
tivity (almost all 0) and health (0-1) are in general 
not a big concern. The flash point 3 of the material is 
another characteristic parameter for its flammability. 

2.3.3 Candidate scintillator mixtures 

All scintillator options consist of a mixture of a sol- 
vent and one or more solute powders, both scintillat- 
ing organic compounds. This section is a compilation 
of properties of the most promising scintillation com- 
pounds and mixtures for LENA. Prospective scintilla- 
tor mixtures have to provide an emission spectra with 
A > 400 nm in order to guarantee large absorption and 
scattering lengths (see Sect. 2.3.1). In addition, this is 
the region where the photomultipliers are most sensi- 
tive. Other critical parameters of the mixture are the 
scintillation pulse shape and light yield, which are im- 
portant for the energy and time resolution of the detec- 
tor. Detailed studies on these parameters for various 
organic liquid scintillation mixtures have been carried 
out in [TP] . 

The mixtures under consideration are both linear- 
alkyl-benzene (LAB) and phenyl- o-xylylethane (PXE, 
with an admixture of non-scintillating dodecane) with 
PPO+Bis-MSB or PMP as fluors. Tabs.0] and [H] sum- 
marize the main properties of solvents and fluors, while 
Tab . |III| provides an overview on the main properties of 
the resulting mixtures. Amongst this selection, LAB 
with PPO (2 g/l) and Bis-MSB (20 mg/t) seems most 
promising. 

The first solvent under consideration, LAB, has 
first come to attention in the R&D studies for SNO+ 
[2"2Tl2"4"] . It is very appealing due to its high trans- 
parency, high light yield, and its low cost. Moreover, 
it is a non-hazardous liquid with relatively high 
flash point of 140°C. The second scintillator solvent 
option PXE has already been tested in the Counting 
Test Facility (CTF) of Borexino [25]. Results look 
very promising as it shows a high light yield, good 
transparency and a high flash point of 145° [H] . 
Dodecane (C12) is a non-scintillating mineral oil 
which is a possible admixture to PXE. It is highly 
transparent and increases the total number of free 
protons in the mixture. This is a key issue, since free 
protons play a key role in the proton decay search 
and in the inverse j3 decay detection channel. 

The emission and absorption spectra of the scin- 
tillating solvents have a significant overlap. In 



Solvent 


LAB 


PXE 


C12 


Physical and Chemical Data 






[TTHT3] 


OY\f i Vr\\c£\] TOTTTlllla 

V_/J..LV_j±±±±\_'Oj± ±±±1_1.±<H.' 


^18^30 


^16^18 


C12H26 


Molecular weight [g/moll 

O LO/ J 


241 


210 


170 


Density p \ks/£] 

J { L Of J 


0.863 


0.986 


0.749 


Viscosity [cps] 


4.2 




1.3 


Flash point [°C] 


140 


167 


83 


Molecular density [10 27 /m 3 ] 


2.2 


2.8 


2.7 


Free protons [10 28 /m 3 ] 


6.6 


4.7 


7.0 


Carbon nuclei [10 28 /m 3 ] 


4.0 


4.2 


3.2 


Total p/e~ [10 29 /m 3 ] 


3.0 


3.2 


2.6 


HMIS Ratings 






[IDCES] 


Health 


1 


1 


1 


Flammability 


1 


1 





Reactivity 











Optical Properties (n,L,t s @J^30nm) 




Refractive index n 


1.49 


1.57 


1.42 


Absorption maximum [nm] 


260 


270 




Emission maximum [nm] 


283 


290 




Attenuation length L [m] 


-20 


12 


>12 


Scattering length £ s [m] 


25 


22 


35 



TABLE I: Overview of the solvent parameters of PXE, 
LAB and Dodecane (C12). The information on physical 
parameters, HMIS (Hazardous Material Identification Sys- 
tem) rating, and refractive index are cited from material 
safety and product specification sheets of the producers 
[11H13] . The HMIS rating quantifies the danger in han- 
dling the liquid (from - save to 4 - dangerous). The molec- 
ular, proton and Carbon densities were computed using 
this information. The absorption and emission maxima of 
PXE and LAB are from |10j . Attenuation and scatterings 
lengths at 430 nm have been experimentally determined in 

Ei in nana. 



3 The flash point is the lowest temperature at which a liquid 
generates sufficient vapor to form a flammable mixture with 
air. Above this temperature, a spark is sufficient for ignition. 



order to shift the wavelength to a region where the 
solvent is transparent, one or more organic solutes 
are commonly added. There are two very promis- 
ing candidates as wavelength shifters for LENA. 
First, a mixture of 2,5-diphenyl-oxazole (PPO) 
and l,4-bis-(o-methyl-styryl) -benzene (bis-MSB) 
or l-phenyl-3-mesityl-2-pyrazoline (PMP). Their 
properties are summarized in Tab.[Tl| lOi |2"6"]. 

The fluor PMP [57] is a primary solute with a large 
Stokes shift of about 120 nm, resulting in a marginal 
overlap of the absorption and emission spectra and 
thus a small self-absorption. The primary fluor PPO 
has a Stokes shift of —60 nm and is therefore usually 
used in combination with a secondary fluor, Bis-MSB. 
It is added in small quantities (mg/£) to avoid self- 
absorption at large wavelengths. The PPO absorption 
spectra significantly overlaps with the emission spectra 
of both the solvents PXE and LAB. Bis-MSB absorbs 
where the primary fluor PPO emits and produces a 
further shift of — 60 nm. 
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Solute 


PMP 


PPO 


Bis-MSB 


Chemical formula 


C18H20N2 


C15H11NO 


C24H22 


Absorption maximum 


294 nm 


303 nm 


345 nm 


Emission maximum 


415 nm 


365 nm 


420 nm 



TABLE II: Properties of the fluors PMP (1-phenyl- 
3-mesityle-2-pyrazolin), PPO (2,5-diphenyl-oxazole) and 
Bis-MSB (l,4-bis-(o-methylstyryl-benzene). Absorption 
and emission maxima are taken from 1 101 1261 



2.3.4 Summary and Outlook 

The final choice of a suited liquid-scintillator mix- 
ture for the LENA experiment must be made most 
carefully. Different physics goals put different con- 
straints on the scintillator properties, making the op- 
timization process difficult. All available solutes and 
solvents cannot be considered independently, but have 
to be tested in specific mixtures. Studies on possible 
scintillation mixtures are not yet concluded. 

The most promising state-of-the-art mixtures are 
listed in Tab. lIIII The solutes under consideration have 
been LAB and PXE, as solvents PPO, PMP and Bis- 
MSB have been investigated. As can be seen, the 
combination of PPO and Bis-MSB provides a good 
light yield for both solutes LAB and PXE. Due to its 
high Stokes' shift, taking PMP as solvent makes the 
use of a secondary fluor needless, however the light 
yield is reduced distinctly Recently, also the solvents 
2,5-diisopropyl-naphtalene (DIN) and n-paraffin drew 
attention and first results concerning decay time con- 
stant and weight look competitive. 

Comparing the fastest decay processes of the differ- 
ent mixtures, one has to take into account both time 
constant and weight relative to the total light emit- 
ted. Obviously, using PMP as fluor, makes the mix- 
ture much slower. In the case of PPO and Bis-MSB 
as fluors, both LAB and PXE show good timing prop- 
erties. PXE is slightly preferable since it is faster and 
the weight of the fastest decay is higher, too. 

With respect to the scintillation light yields, LAB 
and PXE show comparable results. However, the light 
propagation properties of the liquid-scintillator mix- 
ture strongly influence the light collection efficiency of 
the detector. In terms of effective light yield, LAB 
(attenuation length L of ^20 m) is preferable to PXE 
(L = 12m). 



2.4 Light detection 

The photosensors used for the detection of the scin- 
tillation light play an important role for the possi- 
ble physics yield of LENA: The light detection effi- 
ciency affects energy resolution and detection thresh- 
old, while the timing influences event reconstruction 
and particle identification. In the following, we present 
a review of the most important sensor parameters, 
focussing on bialkali Photomultiplier Tubes (PMTs) 



and the existing possibilities to enhance their perfor- 
mance. The properties of alternative light sensors, 
Silicon Photomultipliers and Hybrid Phototubes, are 
shortly addressed. 

2.4.1 Photosensor requirements 

The crucial parameters of photosensors can be split up 
in four aspects: sensor performance, environmental 
properties, availability until start of construction 
and cost-performance ratio. The most important 
properties are: 

Photo detection efficiency (PDE). The PDE is 

determined by the quantum efficiency (QE) of the 
photocathode, the collection efficiency for photoelec- 
trons (pe) as well as backscattering losses of pe at the 
first dynode. In LENA, the baseline value for PDE is 
20 % at 420nm. 

Spectral response. The sensitivity of the photosen- 
sors must match the spectrum of the scintillation light 
arriving at the detector walls. A bialkali photocath- 
ode is arguably the best choice in case of a PMT as 
the maximum PDE corresponds well to the effective 
scintillator emission spectrum. 

Optical coverage (OC) denotes the fraction of the 
detector walls covered with photosensors. Together 
with the detection efficiency of the PMTs, the OC 
determines the overall efficiency of light collection. 
Light-collecting reflective concentrators 33 might be 
used to increase this area beyond the active area, i.e. 
the photocathodes. The present design foresees an OC 
of 30 %. Primarily the energy resolution, but also both 
time and spatial resolution of LENA will depend on 
the product of PDE and OC. 

Time jitter. The timing uncertainty for single pho- 
toelectrons (spe) for individual sensors is vital for the 
overall timing and position resolution of the detec- 
tor. For PMTs, it is given by the transit time spread 
(TTS). 



Afterpulses (AP). AP are spurious pulses occur- 
ring in correlation to primary pulses. Fast after- 
pulses appearing within several tens of ns after the 
primary pulse might impede the resolution of fast dou- 
ble peak structures (e.g. for proton decay, Sect. 6.1). 
APs induced by ions produced by the initial electron 
avalanche occur with delays of several 100 ns to us. 
These ionic APs affect the veto of cosmogenic back- 
grounds by lowering the detection efficiency of muon- 
induced neutrons 34 . Usually the probability of ap- 
pearance is quoted as percentage of primary pulses. 



Dark count (DC). The number of dark counts (per 
photosensitive area) has to be low since it affects po- 
sition and energy resolution by introducing fake hits. 
In extreme cases, it might cause triggers to random 
coincidences. 
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Solvents 


Wavelenj 


gth shifters 


Y[%] 


ni [%] 


n [ns] 






2 nd 


m 


\MM 


1MM 


PXE + 2 g/ 


'£ PMP 




79.1 ± 3.1 


95.9 ± 0.02 


4.15 ± 0.02 


+ 2g, 


/^PPO + 


20 mg/l Bis- 


-MSB 102.0 ± 3.3 


85.3 ± 1.4 


2.61 ± 0.05 


LAB + 2g, 


>t PMP 




83.9 ± 3.0 


85.1 ± 0.9 


8.53 ± 0.15 


+ 2g, 


/£PPO + 


20 mg/i Bis- 


-MSB 99.7 ± 3.2 


77.7 ± 0.8 


5.21 ± 0.005 


DIN + 1.5 j 


;/£ PPO 






86.2 ± 0.2 


6.95 ± 0.02 



TABLE III: Summary of the scintillation properties of different scintillation mixtures. The solvents are PXE and LAB, 
the dissolved wavelength shifters are PMP and a combination of the fiuors PPO and Bis-MSB. The relative light yield 
Y refers to the mixture of PXE with 2g/£ PPO [14], ni and n are the weight and the decay constant of the fast pulse 
component, respectively, as measured in a small cylindrical cell of 2.5 cm in diameter and 2 cm in length. 
Measurements on PXE with PPO show that adding small quantities of the secondary fluor Bis-MSB hardly affects the 
light yield and fast pulse shape component of the scintillation mixture. The values of the scintillation mixtures LAB + 
2g/£ PPO + 20mg/^ Bis-MSB are expectations based on measurements of LAB + 2g/^ PPO 10 . For DIN, the light 
yield measurement is presently carried out. 



Property 


LENA 


Borcxino 


Double Chooz 


SNO+ 


TTS(FWHM) [ns] 


3.0 


3.1 


a few 


<4 


Early pulses a 


<1% 






<1.5% 


Late pulses a 


<4% 


<4% 




<1.5% 


QE for A pcak 


>21% 


>21% 


>20% 




Apeak M 


420 


420 


400 




Optical coverage 


30% 


30% 


13.5% 


54% 


Winston cones 


yes 


yes 


no 


yes 


— > effective area 


xl.75 


x2.5 




xl.75 


Dynamic range b 


spe— >0.3pe/cm 2 






Gain PMTs 


3-10 6 


L10 7 


1-10 7 


L10 7 


spe p/V 


>2 


>1.5 


>2(typ.) 


>1.2 


DC per area [Hz/cm 2 ] 


<15 


<62 




< 25 


Ionic AP (0.2-200 us) a 


<5% 


<5% 




< 1.5% 


Fast AP (5-100 ns) a 


<5% 








Pressure resistance [bar] 


>13 






>3 C 


238 U content [g/g] 


<3-10~ 8 


<3-10~ 8 


<1.2-10" 7 




232 Th content [g/g] 


<L10- 8 


<L10- 8 


<9-10- 8 




nat K content [g/g] 


<2-10~ 5 


<2-10~ 5 


<2-10~ 4 




Detector lifetime [yrs] 


>30 






10 c 



TABLE IV: Overview of the parameter limits for photosensors in LENA and other running or upcoming liquid-scintillator 
detectors 6, 30 32 . spe: single photoclectron, TTS: transit time spread for spe, QE: quantum efficiency, PDE: photo 
detection efficiency, A pea k: peak wavelength of spectral response, pe: photoelectrons, p/V: peak-to-valley ratio, DC: dark 
count, AP: afterpulses. 



Dynamic range. The dynamic range of the response 
of the whole detector has to extend from low energy 
events with mostly only a spe per sensor up to events 
depositing several GeV of energy, corresponding to 
hundreds of pe for the sensors closest to the event. 

Gain, single electron resolution (SER). The am- 
plification gain must be sufficiently large and the single 
electron pulse height resolution sufficient (correspond- 
ing to a large peak-to-valley ratio) in order to distin- 
guish reliably between noise and spe signals. 

Pulse shape, timing effects. Short rise and fall 
times of the spe voltage pulses, in the order of a few 
ns, are advantageous for tracking and position recon- 
struction. Effects shifting the detection time, like early 



pulses, prepulses, and to some extent also late pulses, 
have a detrimental influence on the reconstruction. 

Photosensitive area. For a given optical coverage, 
the granularity of the detector increases with lower 
sensor areas, improving the spatial resolution of the 
photon arrival pattern. Moreover, for smaller areas a 
smaller dynamic range of individual sensors is suffi- 
cient. 

Environmental properties encompass radioactive 
purity of the materials used in the sensors, pressure 
resistance of the sensors placed near the tank bot- 
tom, their long-term reliability over 30+ years and 
susceptibility to magnetic fields. The availabil- 
ity until start of construction will define the candidate 
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sensor types. To calculate the cost-performance ratio, 
it is necessary to know all sensor parameters and their 
effect on detector performance. 

Parameter constraints. Minimum requirements for 
the photosensors to be used in LENA can be formu- 
lated based on the experimental experience gathered 
in the Borexino and Double Chooz experiments, and 
relying on MC studies carried out for LENA. Tab.|TV| 
lists the preliminary limits for LENA in comparison to 
the requirements of other liquid-scintillator detectors. 
MC simulations to further refine the values for LENA 
are ongoing. 

2.4.2 Bialkali photomultipliers 

A detailed description of the functionality and prop- 
erties of photomultiplier tubes (PMTs) can be found in 
[351 136] . In view of the requirements on photosensors, 
head-on hemispherical or plano-convex bialkali PMTs 
with low-background borosilicate glass are presumably 
the subtype best suited for LENA. 

a) Survey of available PMT series 

PMTs are the most natural choice of photosensors for 
LENA, fulfilling all technical requirements (pricing, 
availability and environmental properties) . A compre- 
hensive study to identify the most promising commer- 
cially available PMT models is ongoing. Also, the next 
generation of PMTs featuring high quantum efficiency 
(HQE) photocathodes are being evaluated. So far, the 
characterization has been carried out based on three 
test setups: 

Borexino PMT testing facility. This setup at the 
LNGS was originally used to screen ~2 300 PMTs for 
Borexino 37J. For LENA, it has been used to char- 
acterize a selection of PMTs manufactured by Hama- 
matsu, ranging from 3" to 10" in diameter and fea- 
turing regular bialkali photocathodes: R6091 (3"), 
R6594-ASSY (5"), R5912 (8") andR7081 (10"). Fig.[7] 
shows the transit time distribution measured for the 
8" PMT with this setup: Prepulses and early pulses 
that potentially compromise the arrival times of the 
first photons, and therefore position and track recon- 
struction, as well as late pulses are discernible. 

INFN Milano. This test stand is based on a pi- 
cosecond 405 nm laser (Edinburgh Instruments EPL- 
405) as light source and an 8-bit 1GHz National In- 
struments FADC. Here, Hamamtsu PMTs featuring 
HQE Superbialkali photocathodes (QE«35 % at peak 
wavelengths) have been characterized. PMTs with di- 
ameters from 5" to 10" were tested. 



a Probability of occurrence per primary pulse 
b Assuming xl.75 Winston cones; estimate valid for large sen- 
sor sizes only 

c In ultrapure water, at a water pressure of 3 bar and exposed 
to earthquakes from mining activity 
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FIG. 7: Pulse timing effects in PMTs. Naming conven- 
tion as used in Borexino publications. In this example a 
Hamamatsu R5912 (8") with +1425V applied was mea- 
sured (threshold 0.2 pe) with the Borexino testing facility. 
A close description of the various features of the PMT re- 
sponse can be found in [381 139) . 



Universitat Tubingen. The test stand uses 
fast LED drivers and both FADCs and standard 
TDC+ADC read-out. Measurements of TTS and SER 
were done for 3", 8", 10" and 20" PMTs, includ- 
ing HQE tubes of Hamamatsu R7081 (10"), Photonis 
XP5301 (3") and XP5312 (3"). 

TU Munchen A fourth test stand using a modi- 
fied EPL-405 ps laser and a fast pulsed LED as light 
sources and a 10-bit 8 GHz Acqiris DC282 FADC is 
currently being set up. 

The results of these measurement are currently being 
analyzed and will lead to a first preselection of PMT 
series. 

b) Number of PMTs 

Currently, the LENA reference design aims at an opti- 
cal coverage of 30 %. Light concentrators will be used 
in order to enhance the light collection efficiency. A 
concentrator increasing the effective collection area by 
a factor of 1.75 seems realistic, reducing the necessary 
active photocathode area to 17%. For 8" PMTs of 
normal QE, this corresponds to 63 000 pieces. HQE 
photocathodes could be employed to further reduce 
this number. The scaling of the necessary number is 
demonstrated in Tab.lVl 

In general, the use of smaller-sized PMTs offers 
some advantages: With smaller diameter, the transit 
time spread, the requirements on the dynamic range 
and the rate of ionic afterpulses decrease in general. 
A higher sensor granularity would be beneficial for 
position reconstruction and tracking. On the other 
hand, smaller diameters increase the number of sen- 
sors and electronic channels drastically, which favors 
larger PMTs from the economic point of view. Further 
simulations of the detector behavior are necessary to 
determine the dependence of the physics potential on 
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the PMT diameter and to derive the optimum cost- 
performance ratio. 



Cathode Unarmed 


Light 


HQE-PMTs 


diameter 


PMTs 


concentrators 


-l-concentrators 


3" 


979 100 


574 900 


410 600 


5" 


341 900 


200 700 


143 400 


8" 


114 600 


67300 


48 100 


10" 


85 500 


50 200 


35 800 


12" 


55 900 


32 800 


23 500 


20" 


22 400 


13100 


9 400 



TABLE V: Total number of PMTs necessary to achieve 
an OC of 30 % in the inner detector (ID) and 1.1 % in the 
outer detector of LENA (as in Borexino), depending on the 
photocathode diameter and the use of light concentrators 
increasing the effective detection area by a factor of 1.75 
in the ID. Additionally, the use of HQE photocathodes 
of 35 % peak QE (assuming a collection efficiency of 80 % 
without photoelectron backscattering losses) can further 
reduce the number of PMTs, in this case the OC in the ID 
amounts to 21%. 



2.4.3 Optimization of light detectors 

There are various aspects beyond the selection of 
an optimum PMT that enter the design of the LENA 
detector instrumentation. In the following, a short 
overview of these issues is presented. 

Winston cones. Compound parabolic concentrators 
also known as Winston cones are non-imaging light 
concentrators focussing incident light onto a flat or 
curved surface Mounted to a PMT, they increase 
its effective light collection area (see Fig. [8]). The op- 
timum concentrator shape can be constructed math- 
ematically by rotating the tilted positive branch of a 
parabola around the z-axis. The length of the result- 
ing concentrator determines the area increase ea but 
also the field of view of the PMT: Light incident at 
larger angles than the cutoff angle 9 C to the z-axis in 
a first approximation is not collected by the PMT but 
reflected back into the detector. Winston cones have 
been widely used in neutrino experiments (Tab. |IV[ ). 

Assuming a fiducial volume of 1 1 m radius in LENA 
(e.g. for solar neutrinos, Sect. 4.3), the optimum value 

1.71. 



of 9 C is 50° for 8" PMTs, corresponding to e A 
First MC simulations have shown that for equal OC 
and £a = 2.0, the average pe yield decreases only by 
a few percent. However, the spatial dependence of the 
pe yield increases substantially. Further studies are 
required. 

HQE photocathodes. A further way of increasing 
photocollection efficiency is the use of HQE photocath- 
odes. Conventional PMTs feature a peak QE of not 
more than 25 %. Since a few years, ET Enterprises and 
Hamamatsu Photonics are commercializing PMTs fea- 
turing HQE photocathodes. From Hamamatsu, two 
new types of HQE cathodes are available, and have 



FIG. 8: Winston cone type light concentrator used in 
Borexino, increasing the effective photon collection area 
by a factor of 2.5 (courtesy of Borexino). 



already been tested in the Milano and Tubingen se- 
tups: Super BiAlkali (SBA) photocathodes feature a 
peak QE of 35%, while Ultra BiAlkali (UBA) cath- 
odes reach a peak QE of up to 43%. SBA cathodes 
are available for PMTs of up to 10" diameter, while 
for UBA cathodes the maximum diameter is 2" at the 
moment. However, 12" SBA PMTs as well as larger 
UBA PMTs are expected for the next years. 

Currently, it is not evident whether HQE PMTs 
lower the cost per detected pe. The considerably larger 
price per PMT must be set in relation to the lower 
amount of cabling and electronic channels needed. 
Also, performance aspects like the dynamic range per 
channel or increased afterpulsing probabilities must be 
taken into account. 

Pressure resistance. The pressure tolerance re- 
quired for PMTs is given by the hydrostatic pres- 
sure at the bottom of the tank, depending on the 
scintillator density. The minimum requirement is a 
resistance to 9.8 bar for LAB and 11.1 bar for PXE 



(Sect. 2.3). Another 3 to 4 bar might be necessary to 
withstand an implosion shock wave as it occurred in 
Super-Kamiokande. Simulations will be needed to de- 
termine precise numbers. 

However, the pressure tolerance of available PMTs 
amounts to only ~7bar. The weak spots of the PMT 
glass bulbs are located at the sharp curvatures at the 
PMT neck and base. One option is to increase the 
thickness of the glass casing, which would increase the 
weight and the amount of radioactivity introduced by 
the PMTs, and is deemed to be demanding by manu- 
facturers. 

A promising alternative is the use of pressure- 
withstanding encapsulations housing the PMTs at at- 
mospheric pressure. A widely used approach in HE 
neutrino physics are spherical glass bulbs. However, 
the currently preferred solution is an adaptation of the 
Borexino Outer Detector PMT encapsulations. Their 
design consists of a conically shaped metal housing 
enclosing the PMT body combined with a transpar- 
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ent window 5 in front of the photocathode (Fig. [9]). 
This design also allows for an easy integration of the 
light concentrators and the u-metal shielding. Design 
work based on finite element simulations has recently 
started and will result in a prototype encapsulation for 
pressure testing. 




FIG. 9: Illustration of a LENA PMT encapsulation. 



2.4.4 Alternative photosensor types 

In spite of the broad range of photosensors cur- 
rently available or under development, the number of 
options for use in LENA is fairly limited. Consider- 
ing the requirements regarding time resolution, avail- 
ability within the next decade and long-term reliabil- 
ity, only Silicon Photomultipliers and Hybrid Photo- 
tubes are promising alternatives to conventional bial- 
kali PMTs. While their performance is mostly supe- 
rior, it remains uncertain whether the price per photo- 
sensitive area is low enough and if their performance 
is sufficient in all respects. 

Silicon Photomultipliers (SiPMs), also known as 
Multi-Pixel Photon Counters (MPPCs) are arrays of 
avalanche photo diodes, operated in limited Geiger 
mode [H]. Typically there are 100-1000 pixels per 
mm 2 , which are electrically connected in parallel, so 
the total signal is proportional to the number of cells 
hit by one or more photons. 

SiPMs have several advantages over conventional 
PMTs: The QE of SiPMs can reach 70 to 80%, with 
spectral response very close to the PPO emission spec- 
trum and PDEs over 70 %. spe time resolutions as low 
as 60 ps have been observed [12], dominated by the jit- 
ter of photons being detected in different cells of the 
array. With gain similar to conventional PMTs (10 5 



5 In the Borexino design, this window is a thin PET foil. In 
LENA, this window will have to absorb the external pressure, 
so a more resistive material like acrylic might be used. 



to 10 6 ), their SER is much smaller, which allows ex- 
act counting of small numbers of photons at tempera- 
tures of 10-15°C. Regarding environmental properties, 
SiPMs are insensitive to magnetic fields, which would 
allow to magnetize the LENA detector. Furthermore 
their radioactive contamination can be expected to be 
extremely low and they are very slim, possibly allow- 
ing to reduce the buffer thickness or even omit the 
buffer, which would permit to increase the target vol- 
ume substantially or greatly reduce costs. In addition, 
the pressure tolerance of SiPMs probably is very high. 
Also, their bias voltage is below 100 V, allowing for a 
much simpler voltage supply system. 
On the other hand, there are several disadvantages: 
The size of the active area is very small, currently 5 x 
5 mm 2 for the largest commercially available SiPMs. 
With larger sizes the number of channels decreases 
and the PDE increases. It is reasonable to assume 
that lxl cm 2 SiPMs will become available in the near 
future and many SiPMs could be combined into local 
clusters to save channels. However, if one wants to 
preserve the photon counting ability, the active area 
per channel is limited to about 1cm 2 . Another prob- 
lem is the dark count, which can approach MHz fre- 
quencies, disturbing the reconstruction of events and 
increasing with active area. Also, temperature stabi- 
lization is necessary to prevent drifts in gain. Finally, 
it is not yet clear, whether the price per area will be 
competitive to PMTs. 

Hybrid Phototubes (HPTs). The basic design of 
HPTs combines a large area hemispherical photocath- 
ode for photoelectron conversion, a HV field accel- 
erating these electrons towards a small luminescent 
screen and a small diameter PMT (1") reading out the 
scintillation signal. HPTs have been used successfully 
in the Lake Baikal Neutrino Telescope, featuring the 
QUASAR-370 of 15" diameter produced in the 1990s 
43 45 . However, they are currently out of produc- 
tion. 

Compared to conventional PMTs, HPTs would fea- 
ture a range of advantages: The angular acceptance is 
2tt, the transit time spread is typically less than 1 ns, 
prepulses, ionic AP, fast AP and late pulses occur only 
on very low levels, the DC is very low and the suscep- 
tibility to magnetic fields is greatly decreased. How- 
ever, it is currently unknown, whether the dynamic 
range of HPTs is sufficient for use in LENA, as per 
detected photon about 25 pe hit the small PMT. One 
solution would be to use less dynodes as it was done in 
the TUNKA Air Cherenkov Array experiment for the 
QUASAR-370G 05], which on the other hand could 
affect the discrimination from noise. 

2.4.5 Conclusions 

The benchmark value for the photon detection effi- 
ciency of LENA is 6%, arising from the product of 
optical coverage (OC=30%) and PMT photo detec- 
tion efficiency (PDE=20%). 
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Bialkali PMTs currently appear to be the most 
promising choice for photosensors, basically fulfill- 
ing all requirements concerning detection performance 
and long-term reliability, and low in cost per photoac- 
tive area. PMT diameters from 3" to 12" are consid- 
ered. It is foreseen to equip the PMTs with Winston 
cones, the factor of area increase ranging from 1.6 to 2, 
and u-metal shielding. If available and cost-effective, 
HQE PMTs might be used to reduce the necessary 
OC. Furthermore, a pressure-absorbing encapsulation 
will be necessary to protect the PMTs at the bottom 
of the detector from implosion. 

Based on the interim results of the performance 
tests, further tests on the most promising series of 
Hamamatsu Photonics and ET Enterprises will be 
conducted with larger numbers. This study will also 
extend to new releases as the 12" PMT by Hamamatsu 
(R11780) and its upcoming HQE version. The middle- 
term aim is the development of a prototype optical 
module consisting of PMT, encapsulation, u-metal and 
light concentrator by the end of 2014. 



Parameter 


T TP 


rirj 


Number of channels (8" PMTs) 


45 000 


45 000 


Time resolution (pulse on-set) 


<1 ns 


<1 ns 


Dynamic range 


0-30 pe 


0-300 pe 


Dead time per channel 


< 100 ns 


FADCs 


Channel buffer size 


-100 


FADCs 


Number of FADC channels 


X 


10 000 


Sampling rate 


X 


500MS/s 


Voltage resolution 


X 


2x8 bit 



TABLE VI: Requirements for the electronic read-out in 
LENA. The table lists both the minimum requirements 
for low-energy (LE) neutrino detection and the optimal 
configuration for high-energy (HE) beam physics. 



In the following we describe two very different op- 
tions for the readout electronics: a full FADC option, 
and a custom ASIC option. We believe that the first 
would be preferable, but the second could well meet 
the basic scientific requirements and be less expensive. 
For these reasons we include both. 



2.5 Read-out electronics 



2.5.2 Full FADC readout 



In view of the huge number of photosensors needed for 
the next generation detector the only solution techni- 
cally available today arc PMTs. With a properly de- 
signed read-out electronics, PMTs can guarantee the 
required time resolution, charge resolution and dy- 
namic range. Two different options for readout elec- 
tronics are described, one relying on FADC readout 
of all channels, the other featuring customized ASIC 
boards inside the detector servicing small arrays of 
PMTs. 

2.5.1 Minimum requirements 

The broad and rich physics scope of LENA puts rather 
severe requirements on electronics and data acquisi- 
tion. The relevant performance parameters are sum- 



marized in Tab. VI Particularly, the goal to perform 
both low energy solar neutrino physics (especially pp 
neutrinos) and high energy beam neutrino physics 
force a very large dynamic range on the signal am- 
plitude because the system must work both on single 
photoelectron mode and with very large signals cor- 
responding to several hundreds of photoelectrons per 
channel. Moreover, a precise determination of the rel- 
ative time of each PMT hit with sub-ns resolution is 
necessary for the spatial reconstruction of the events, 
for the pulse shape analysis and in order to disentan- 
gle the different event topologies that are associated 
to neutral and charged current interactions of elec- 
tron, muon and tau neutrinos. Zero dead time is a 
must to avoid the loss of interesting time correlated 
events, and a very flexible triggering system is highly 
desirable, taking into account again the very broad 
scientific scope and the long expected lifetime of the 
experiment, which may lead to unanticipated physics. 



The scientific requirements can be met by the use of 
suitable FADCs with the right pulse height resolution, 
sampling speed, and the combination of on-board zero 
suppression, software trigger and careful synchroniza- 
tion of the boards. 

At the time of writing (many years before the be- 
ginning of data taking) there are several products on 
the market that basically meet these requirements; 8- 
bits, 2 GS/s FADC boards are already available, and 
better products are very close to be. This solution is 
still quite expensive today, but it is very likely that it 
will not be at the time when LENA will need it. 

The system could work as follows: each PMT should 
be connected to a simple linear front end that should 
perform high-voltage decoupling, some amplification 
and some shaping. The details of this rather standard 
front end electronics must not be worked out here. 
It is probably cost-effective to imagine this front end 
built into a custom chip mounted close to the PMTs. 
The output signal of this front end should be sent to 
the FADC. The FADC board will sample it at the 
correct speed (1 GS/s seems appropriate but 2 GS/s 
is feasible) with 10-12 bit precision (14 bit is already 
available and probably more than adequate) , and store 
into a local memory. An on-board fast FPGA will 
perform the zero suppression, keeping the data only 
around valid PMT pulses, and storing also the time 
stamp of the pulse. Data from each pulse will then be 
stored into another internal memory ready for read- 
out. 

Each FADC board shares a common distributed 
clock with all other boards, so that sampling is syn- 
chronous throughout the whole detector. The acquisi- 
tion software will read continually the time stamps 
(time stamps only at first), perform the triggering 
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functions and decide whether the 'event' should be 
read out or not. In case the read-out is needed, a cus- 
tom designed protocol between the FPGA onboard, 
the FADC and the readout software will retrieve all 
samples, and the event will be stored on disk. 

This architecture should guarantee the maximum 
possible information (the fully digitized shape of each 
PMT pulse for the whole detector in a programmable 
time window around the event), and zero dead time. 
Also, the use of the time stamps for triggering should 
keep the complete data flow to a very reasonable level. 
This fact is easy to prove: the PMT activity is largely 
dominated by dark noise; physics events of any kind 
contribute a much lower amount of data. A typical 
large PMT has a room temperature dark noise of the 
order of 1 KHz. Assuming 2 kHz to have some contin- 
gency, the expected total activity of LENA is of the 
order of 10 8 Hz (conservative). This means 400Mb/s 
of time stamps to handle, a very reasonable number. 
The event triggering rate will be dominated by 14 C 
and will reach 10 4 Hz with a few hundreds of PMT 
hits (again conservative); the expected data flow will 
be of the order of 100 Mb/s— 1 Gb/s, again an easy 
number to handle. 

Sufficient onboard memory is required to buffer the 
sampled signals and timestamps. By employing com- 
mercially available memory aimed at the consumer 
market, available FADC modules already offer an on- 
board memory of 512 MB or more per channel at rea- 
sonable rates. If assuming a dynamic range of 12 bit 
and a sampling rate of 2 GHz with a time window of 
300 ns saved in the buffer, less than 1 kB of memory 
is used per trigger. With an average trigger rate of 
2 kHz, 250 seconds of live data can be buffered before 
data loss occurs. Data can be read continuously from 
the memory while new events are stored at the same 
time, thus leading to a dead-time free data acquisition. 
To prevent data loss in the rare cases where the trig- 
ger rate can be up to ten times the dark noise for the 
duration of several seconds (e.g. a galactic Supernova 
sufficiently near), normally only a fraction of the total 
buffer size should be used before readout is initiated. 
A further possibility is to set a programmable thresh- 
old for the trigger rate. If it is exceeded for several 
seconds, the onboard logic will stop to store the full 
pulse shape to prevent a buffer overflow; instead, it will 
save only the timestamp and the energy deposition de- 
rived by integrating over the sampled pulses. Based 
on this, complete event reconstruction is still possible, 
while the memory usage of each trigger is reduced to 
^12 byte per channel. In this mode, more than 400 
seconds of live data could be buffered at a trigger rate 
of 20 kHz, even if only 100 MB are reserved for this op- 
eration mode. So, a FADC-based DAQ would record 
the essential information even under the extreme con- 
ditions of a galactic SN neutrino burst, and should be 
seriously considered. 

It is not crucial here to decide whether the FADC 
board should be commercial or custom. Most likely, 



commercial boards will be cheap enough to be prefer- 
able, but this is not a decision for today. It is however 
important to distribute the FADC boards close to the 
PMTs, in order to minimize cable length. Data from 
the FADC chassis will be collected via optical fibers 
which are of course cheap, simple and reliable. The 
distance between front end and FADC should be min- 
imized to avoid spurious noise. 

This architecture is good both for the main detec- 
tor and for the muon veto, which of course makes the 
system simpler. 

A final note on dynamic range: Even with a large 
number of bits, it may be difficult to guarantee a good 
linearity for the signal range from single to hundreds 
of photoelectrons. In this case it should be considered 
as an option to have a front end with doubled output 
and two different gains. This will increase the number 
of channels (not all channels must be duplicated) but 
would guarantee very good performance. A trade-off 
between cost and performance is mandatory on this 
point, and a decision must be taken after a full simu- 
lation. 




FIG. 10: On the left the demonstrator of the PMm 2 R&D 
program that is going to be tested with its electronics sys- 
tem in the MEMPHYNO prototype (right). 



2.5.3 Custom ASIC read-out 

The coverage of a large area with PMTs at a "low" 
cost can be met by a readout electronics integrated 
circuit (called ASIC) for groups of PMTs. The devel- 
opment of such electronics is the aim of a dedicated 
French R&D program, called PMm 2 @B]. This R&D 
program was initially stimulated by the MEMPHYS 
Water Cherenkov project [471 14"5] , 

PMm 2 intends to realize a new electronics board 
dedicated to a grouped acquisition of a matrix of 
16 PMTs. Each matrix will have a common board 
(PARISROC) for the distribution of high voltage and 
for the signal readout. The circuit under development 
allows to integrate for each group of PMTs: a high- 
speed discriminator on the single photoelectron signal, 
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the digitization of the charge (on a 12-bit ADC) to pro- 
vide numerical signals, the digitization of time (on a 
12-bit TDC) to provide time information, a channcl- 
to-channel gain adjustment and a common high volt- 
age. DAQ system, trigger and mechanical integration 
of the matrix is currently under development in a joint 
effort by teams of the Laboratoire de l'Accelerateur 
Lincarc (LAL), the Institut de Physique Nucleaire in 
Orsay (IPNO), and the APC Paris. 

To test the system with real physical signals, the 
Water Cherenkov prototype detector MEMPHYNO is 
presently under construction at the APC Laboratory 
in Paris [49 . The aim is to install a complete 4x4 
array of PMTs and the complete electronics and ac- 



quisition readout chain in a cubic water tank of 2 m 
edge length. A muon hodoscope based on four scintil- 
lator planes will provide the trigger for cosmic muons 
as well as muon track information. Based on this, 
MEMPHYNO will evaluate the system trigger thresh- 
old, the track reconstruction performance in water and 
the properties of the PMTs. Fig.[lO] displays both 
the PMm 2 matrix and a conceptional view of MEM- 
PHYNO. 
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3 Detector performance 



The expected performance of the final LENA detec- 
tor can be extrapolated to a large extent from the 
currently running large-volume liquid-scintillator de- 
tectors, first of all of course Borexino and KamLAND. 
However, a precise assessment of the accuracy achieved 
in the reconstruction of neutrino events and back- 
ground levels present in the detector requires detailed 
Monte Carlo simulations. The necessary input on the 
properties of scintillator, PMTs and detector electron- 
ics is presented in Sec. [2] 

We start out with a set of baseline parameters for 
the detector components in Sec. |3.1| which are later 
on implemented in the phenomenological studies of 
Sees 4][6 The precision of vertex and energy recon- 
struction is discussed in Sec. |3.2| A special emphasis is 
given to the reconstruction of extended particle tracks 
at higher energies in Sec. |3.3| motivated by its im- 
portance for long-baseline oscillation physics at GeV 
energies. A discussion on the backgrounds affecting 
the neutrino and rare event searches is presented in 
the corresponding phenomenological sections. 



3.1 Baseline parameters 

The performance of the LENA detector will depend 
first of all on the quality of the liquid scintillator used, 
but also on the granularity of and active area covered 
by photonsensors as well as the time resolution and 
dynamic range of the readout electronics. Therefore, 
all these aspects have to be included in a Monte Carlo 
simulation of neutrino and background events in the 
detector. For this purpose, a GEANT4-based frame- 
work has been set up for LENA, using the geometry 
described in Sec. [2| an d a set of baseline parameters 



presented in Tab. VII Several steps have to be run 



through in order to obtain the detector response: 

Detector environment. Central to the simulation is 
the definition of the materials and dimensions making 
up the detector. In accordance with the design pre- 
sented in Sec. [2] the scintillator volume is fixed to 26 m 
diameter and 96 m height, surrounded by a mantle of 
inactive buffer medium with a width of 2 m. While 
earlier studies presume PXE as solvent, more recent 
works use LAB as baseline. Corresponding material 
properties like the number densities of electrons and 
protons are listed in Tab. [T] PMTs are assumed as 
photo detectors, their photocathode planes located at 
14.5 m radius. At 15 m a steel tank of 10 cm separates 
buffer and outside water volume, adding another 2 m 
of (active) shielding. 

Event vertex. While most event vertices at low en- 
ergies are rather simple, mostly one charged leptons 
and at times a 7-quantum or free nucleon in the fi- 
nal state, high-energy neutrino interactions may cre- 
ate additional particles like pions and often break up a 
Carbon nucleus acting as interaction target. GEANT4 
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TABLE VII: Baseline parameters of the detector perfor- 
mance entering the LENA simulation. For target masses 
and interaction center densities, values are given separately 
for LAB and PXE. 



is well suited to follow the path of these particles in 
the target volume, creating a realistic picture of the 
energy deposition inside the scintillator. 

Light emission. The number of photons produced 
per energy deposit, i.e. the light yield of the scintil- 
lator, corresponds to roughly Y p w 10 4 photons per 
MeV. The exact value may vary by the order of 10 % 
depending on the scintillator solvent and fluor con- 
centrations (cf. Tab. III|. Moreover, there is a strong 



dependence of the light output on the dE/dx of the 
ionizing particle, resulting in a quenching for heavy 
particles (protons, a) compared to light ones (elec- 
trons). This effect is taken into account by the Birks' 
formalism, using Ub = O.Olcm/MeV []. 

Emission profile, the light is not emitted instanta- 
neously, but following a time profile described by a su- 
perposition of exponential decays. The simulation uses 
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a description based on three components, which are 
close to the ones of the favored LAB+PPO+bisMSB 
cocktail. Timing properties depend mostly on the fast 
fluorescence time T\ . The relative weight of these com- 
ponents depends on the ionizing particle, protons and 
a's featuring enhanced secondary components com- 
pared to electrons. 

Light propagation. While the scintillation light is 
produced in a relatively broad span of wavelengths 
that is defined by the emission spectrum of the sec- 
ondary fluor bisMSB, an effective description of the 
light transport through the liquid can assume that the 
light is shifted to the wavelength range around 430 m. 
As the scintillator transparency increases strongly 
with wavelength, most of the light is transmitted and 
detected towards the long-wavelength part of emission 
spectrum and PMT sensitivity. Photon absorption, 
rcemission as well as Rayleigh scattering are repro- 
duced by the code. The corresponding parameters are 
chosen conservatively at the lower end of the expecta- 
tion for LAB, which is also a good approximation for 
the case of PXE . 

Light detection. The baseline design for LENA as- 
sumes a photocoverage of 30% and a photosensitiv- 
ity of 20%. In many cases, it is sufficient to reduces 
the light yield to 6 % and count all photons reaching 
the detector walls as detected. However, the baseline 
scenario assumes 13,000 unarmed 20"-PMTs evenly 
distributed over the inner surface. This choice can 
be regarded as conservative, as a greater number of 
smaller PMTs resulting in the same photocoverage 
will only increase the granularity and resolution for 
photon arrival times. Other configurations like, e.g., 
8" -PMTs equipped with light concentrators) can be 
implemented, taking also the limited photon accep- 
tance angle into account. Arrival times of individual 
photons are smeared by a Gaussian of 1 ns width (Iff), 
corresponding to the PMT specifications presented in 
Sec. [231 

Electronic readout. Up to know, the effects of elec- 
tronics are described only coarsely. However, a base- 
line design comparable to the one of Borexino is of- 
ten applied, featuring a time resolution better than 
the PMT time jitter for the arrival time of the first 
photon per channel and integrating the number of all 
subsequent photons in the same channel to provide 
an integral charge. For the moment, a similar scheme 
is used for most of the simulations aiming at particle 
track reconstruction at sub-GeV energies. However, 
studies using the equivalent of a FADC re cordi ng each 
individual channel are also on-going (Sec. |2.5| ). 

Based on this chain, the photoelectron (pe) yield per 
deposited energy, Y pe (r) can be obtained as a function 
of the vertex position inside the target volume. Un- 
like current spherical detectors were the light yield is 
approximately uniform over the whole detection vol- 
ume, a significant dependence is expected for LENA 
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FIG. 11: Map of the photoelectron yield Y pe (r) that 
features a considerable dependence on the radial and z- 
positions of the event vertex. 



mainly because of the larger detector dimensions and 
the cylindrical shape. Fig. [TT] shows Y pe in the height- 
radius plane for one fourthof the detector. For the 
baseline values shown in Tab. |VII[ the mean value for 
uniformly distributed events is (Y pe ) = 238pe/MeV. 

This conservative value will be significantly ex- 
ceeded if the liquids prove to be more transparent. For 
an attenuation length of 15 m, well within the range 
feasible for LAB, the mean yield increases to (Y pe ) = 
340pe/MeV. Yields as large as 450pe/MeV might be 
obtained in optimistic scenarios, almost equalling the 
pe yield achieved in Borexino. 

3.2 Low-energy vertex reconstruction 

The reconstruction of the primary vertex position as 
well as the reconstruction of the event time and the 
visible energy 6 of low energy events is a prerequisite 
for all further analysis. 

Accounting for the spacial and time resolution of 
LENA, events with visible energies below 10 MeV can 
be considered as point-like in space-time. Hence, they 
are described by a set of five coordinates 

X = (r,t ,T) 

where r is the position of the event, to is the time when 
the event occurred and T is the visible energy. These 
parameters are determined by a negative logarithmic 
likelihood fit to the hit times of the first photons as 
well as the number of photons detected on each PMT. 

The main fit is seeded by a series of previous anal- 
ysis steps. First, as an estimate for r, the charge 
barycenter r& is determined by an approximate fit to 



6 The visible energy of a particle in the scintillator is defined 
as the energy an electron has to deposit in the scintillator to 
create the same number of photons. 
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the integrated number of photoelectrons collected on 
each PMT. Using r& as input, the distribution of TOF- 
corrected first-photon arrival times is used to gain a 
seed for t . Finally, an approximate value for T is 
obtained from the total charge seen in the detector, 
taking into account the spatial inhomogeneity of the 
photoelectron yield. 

The probability density function (PDF) used in the 
main fit considers the time resolution and the finite di- 
mensions of the PMTs, the expected dark count rate, 
the fluorescence times of the scintillator as well as ab- 
sorption and scattering in the scintillator. Currently, 
the PDF assumes LENA to be a cylinder of infinite 
height. Hence it does not describe events close to the 
lid/bottom of the detector correctly. 

The fit based reconstruction has been tested with 
several datasets of point like events generated with 
the LENA Monte-Carlo simulation 50J . It produces 
stable fit results without significant systematic shifts 
for \z\ < 25m and T G [200keV, lOMeV]. To speed up 
the simulation as well as the reconstruction, 20" PMTs 
(without Winston cones) have been used throughout 
the study. PMT resolution is included by a Gaussian 
transit time spread of Ins (la), dark noise is neglected. 

shows 
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Figure 12 



the results for 12 000 electron 
events of 1 MeV, distributed randomly in the detec- 
tion volume (\z\ < 25m). As LENA has a cylindrical 
shape, the resolution in z-direction, a z = (9.99 ±0.05) 
cm, differs from the resolution in x- and y-direction, 
&x,y = (8.25 ± 0.05) cm. Furthermore, as no Winston 
cones are used in the simulation, the photoelectron 
yield increases with rising radii. Hence, the resolution 
at the border of the active volume is better than the 
resolution in the center of the detector. The difference 
is of order 20% at T = 1 MeV. 

Averaged over the whole detector, the number of 
detected photoelectrons is N pe = 220 at 1 MeV. The 
resulting energy resolution follows the expected behav- 
ior of (AE/E) = (N^)- 1 / 2 . Also the spatial resolu- 
tion improves with energy, but the energy dependence 
is more complicated. 

At low energies, the dark noise of the PMTs has 
a degrading influence on the vertex reconstruction. 
For an integral dark rate of 45 kHz (roughly 4 kHz per 
PMT) , the resolution is decreased on the percent level. 
While this effect plays almost no role above 1 MeV, 
it becomes much more prominent at lower energies, 
reaching approximately 30 % for electron events of 200 
keV deposited in the center of the detection volume. 

In the final experiment, the ambiguities of position 
and energy reconstruction will be resolved by a calibra- 
tion campaign based on radioactive sources. Similar 
campaigns have been conducted in both KamLAND 
and Borexino detectors (e.g. [ST]). The point recon- 
struction in liquid scintillators has proven to be a very 
powerful tool in analyses. The future development of 
the presented algorithm will be focused on including 
events close to the lid and the bottom of the detector. 
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FIG. 12: Overall performance of the vertex reconstruc- 
tion at 1 MeV: The distance of the reconstructed to the 
true event position is shown in the upper histogram. The 
lower panel shows the deviation of the reconstructed en- 
ergy. Both distributions are fitted with the appropriate 
functions, the resulting resolutions are indicated in the cor- 
responding panels. The simulated events were randomly 
distributed in a slice of the active volume with \z\ < 25 m, 
the effects of dark noise are neglected. 



3.3 GeV event reconstruction 

The reconstruction of particle momenta is a prereq- 
uisite for the analysis of atmospheric and accelera- 
tor neutrinos. It was realized only recently that liq- 
uid scintillator detectors — opposed to general opinion 
— feature this capability, provided the particle track 
length exceeds several tens of centimeters. This sec- 
tion reflects the state-of-the-art of Monte Carlo simu- 
lations that investigate energy and angular resolution 
both for single and multiple-particle events. 

3.3.1 Introduction 

The analyses of beam and atmospheric neutrinos re- 
quire a neutrino detector to be capable of reconstruct- 
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ing both energy and momentum of the incoming neu- 
trino. Depending on the exact task, it may also be nec- 
essary to identify the flavor (or antiflavor). At higher 
energies, charged current interactions will excite reso- 
nances and start to scatter inelastically, creating not 
only a lepton but also pions and heavier hadrons in the 
end state. Moreover, background signals due to beam 
contaminations or flavor-insensitive neutral current in- 
teractions must be identified and rejected. Detectors 
that fulfill these requirements are usually highly seg- 
mented or feature excellent tracking capabilities (like 
liquid argon time projection chambers). 

However, at energies not exceeding a few GeV, event 
vertices and backgrounds are less complex, and low- 
energy neutrino experiments become viable candidates 
for a far detector. This is most impressively demon- 
strated by the Super-Kamiokande experiment, that 
found neutrino oscillations both by analyzing atmo- 
spheric neutrinos and the K2K neutrino beam, and 
is currently serving as far detector in the T2K beam 
experiment searching for #13 [52TI54] . 

At first glance, it seems unlikely that unsegmented 
liquid-scintillator detectors might be used in the same 
way: An imminent feature of water Cherenkov detec- 
tor is the directional information coded in the orien- 
tation of the Cherenkov cone. Using this information 
is possible even for particles close to the Cherenkov 
threshold. Opposed to that, scintillation light emit- 
ted by low-energy events is distributed isotropically, 
bearing no directional information at all for the quasi- 
pointlike events. 

However, high-energy particles deposit their energy 
over macroscopic distances. In liquid scintillator, they 
will create a track of ionization extending for tens 
of cm or even meters, leading to a distortion of the 
spheric light front emerging from the track. As il- 
lustrated in Fig. |13| the superposition of spherical 
light waves emitted along the particle track creates a 
light front which resembles the Cherenkov light cone, 
adding a spherical backward running front to the v- 
shaped forward front. 

The possibility to exploit the inherent directionality 
has been neglected for a long time, as the deformation 
of the light front is too small for low-energy neutrinos. 
Only in the rejection of cosmic background, namely 
in the muon track reconstruction algorithms of Kam- 
LAND and Borexino, the arrival time patterns pro- 
jected by the light front on the PMTs are exploited 
[3H [SS]- Based on this, the orientation of the track 
can be reconstructed with an astounding accuracy: 
The Borexino track reconstruction achieves an angu- 
lar resolution of 3° for muons crossing the scintillator 
volume [3"4"] . 

Only recently, the basic possibility to reconstruct 
the track direction by exploiting the peculiar shape of 
this light front has been brought to attention [551 ISZ] : 
Since then, tracking algorithms based on the Monte 
Carlo simulations of GeV neutrino events have con- 
firmed the basic notion that energy and momen- 




FIG. 13: Construction of the first photon surface (blue) by 
superposition of spherical waves (red) created by a particle 
transversing the scintillator (black) 



turn of the end-state particles and finally of the in- 
cident neutrino can be resolved. Moreover, the stud- 
ies demonstrated that the accuracy of the reconstruc- 
tion could in principle exceed the performance of wa- 
ter Cherenkov detectors due to the much larger light 
yield. 

In the following, two different approaches to the 
reconstruction of GeV neutrino events in LENA are 
presented: The first one is a GEANT4-based evalu- 
ation of tracking single electrons and muons at sub- 
GeV energies, the second investigates the possibility 
to resolve more complex interaction vertices in the 
1—5 GeV range, based on a specifically written pro- 
totype code [57] , 

3.3.2 Tracking in the sub-GeV range 

The reconstruction of a particle track must rely on 
the projection of the Fermat surface depicted in Fig. [13] 
on the surface composed by the PMTs mounted to the 
detector walls. As a start point, the patterns of first 
photon arrival times and integrated charge per PMT 
can be exploited. 

Fig. 14 shows an example for a single particle 
track: A 500 MeV muon traveling from the center 
of the detector towards the wall. The event was 
created using a Geant4 based simulation of the 



LENA detector. Depicted in figure 14 1) is the charge 



distribution which features only a slight asymmetry 
due to the displacement of the track's center of charge 
with respect to the symmetry axis of the detector. 
Nevertheless, the charge signal of the PMTs can be 
used to obtain the track's barycenter. This allows 
removing most of the dependence of the photon arrival 
times on the track position by a time of flight (TOF) 
correction with respect to the barycenter. The 
resulting distribution is depicted in figure 14 d). 
Only the first 11ns (from -8 ns to +3ns of TOF 
corrected hit time) arc shown to enhance clarity. The 
observed distribution is clearly anisotropic and can 
be used to get a rough estimate on the track direction. 

A more precise reconstruction of the track is 
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FIG. 14: A 500 MeV muon in LENA. On the left, the color coded information is the charge seen by each PMT, while the 
hit time of the first photon at each PMT is shown on the right, applying a time of flight correction with respect to the 
charge barycenter of the track. 



achieved by determining the track parameters using 
a negative logarithmic likelihood fit to the integrated 
charge and the first hit times of each PMT. This is 
done employing the continuous slowing down approx- 
imation i.e. neglecting any kind of statistical fluctua- 
tions of the track. The number of parameters required 
to characterize a track is therefore reduced to seven: 
The kinetic energy of the particle, the start point of 
the track, the track direction and the time when the 
particle was created. 

Fig. [15] shows the results for single 300 MeV muons 
traveling from the center of the detector towards the 
wall. The resolution obtained for the start point of 
the track is in the order of a few centimeters and the 
start time of the event can be determined with sub- 
nanosecond accuracy. The obtained angular resolution 
is in the order of a few degrees. The results obtained 
for electrons are in the same order of magnitude but 
tend to be slightly worse compared to muons due to 
the higher statistical fluctuations of electron tracks. 
Track reconstruction yields useful results for kinetic 
energies down to 100-200 MeV for single muons and 
down to ^250 MeV for single electrons. The perfor- 
mance for muons at energies of order 100 MeV is lim- 
ited as muons are no longer minimum ionizing which 
leads to very short track lengths of a few ten centime- 
ters. The low energy limit for electrons on the other 
hand is due to the increasing statistical deviations of 
the track from the straight line. 

Quasielastic CC j/^-events can be reliably distin- 



guished from quasielastic J'e-events, by tagging the de- 
cay of the muon produced in the former case. As the 
spectrum of the Michel electrons exceeds the energies 
of background processes, most notably the energies re- 
leased due to neutron capture on hydrogen or carbon, 
a cut on the visible energy of the muon decay efficiently 
eliminates background processes mimicking a muon 
decay. This principle has been verified by simulat- 
ing 7600 quasielastic v e interactions. As non of these 
events was wrongly assigned a r^-event, the probabil- 
ity to wrongly identify a quasielastic ^ e -event as a Vur 
event is smaller than 4- 10 -4 at 95% C.L.. At the same 
time, the efficiency is roughly 70%. This method to 
reject background is of course insufficient for inelastic 
interactions where pions are produced. These events, 
though suppressed at beta-beam energies, require fur- 
ther studies to devise a reliable discimination. 

3.3.3 Tracking in the 1—5 GeV range 

Here we consider a charged current (CC) neutrino 
event producing a charged lepton. In the lower part of 
this energy range the scattering is usually quasielas- 
tic, but at higher energies single and multi-pion events 
dominate. Typically a recoil nucleon is emitted but 
intranuclear collisions and excitations may absorb en- 
ergy and momentum. 

The reconstruction of the 1-5 GeV event is done 
with a prototype code that tries to find the test event 
giving the best fit with the recorded PMT data of the 
"true" event. The code simulates the scintillation light 
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x = (-0.3± 1.8) cm y = (-0.1± 1.9) cm z = (0.0± 2.4) cm 
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FIG. 15: Results obtained by reconstructing 300 MeV muons created in the center ol the detector and traveling in 
negative x direction (500 events). The upper row shows the results for the start point of the track, the lower row shows 
the reconstructed start time (left), the angular deviation of the reconstructed track from the Monte Carlo truth (center) 
and the kinetic energy of the muon (right). 



emission from all the secondary particles but no ter- 
tiary particles nor nuclear physics. 

We found that the single lepton tracks can be re- 
constructed very well. Also 2 sufficiently long (> 
O(50 cm)) angularly separated tracks can be distin- 
guished easily. Events consisting of 3 tracks can be 
reconstructed if the tracks are clear, long (> 0(1 m)) 
and well separated, though tracks of a few 10 cm in 
length remain unseen. Almost parallel tracks are al- 
ways hardest to distinguish. 

Events with 4-5 tracks are very challenging, and 
can be reliably reconstructed only in special cases, like 
well-separated tracks longer than a meter, with addi- 
tional signals from particle decays. The lepton track 
itself can be distinguished in almost every case. 

For the studied energy range from 1 to 5 GeV, the 
identification of the lepton flavor shows no ambigu- 
ities. No misidentification occurred throughout the 
testing for straight leptons with any number of sec- 
ondary particles. The muon and electron signals are 
very different and there is no way to confuse clear lep- 
ton tracks at GeV energies unless there are some rare 
tertiary processes. Also the statistical fluctuations in 
photon emission or detection cannot cause such errors 
with any reasonable probability. 

The position of the interaction vertex can be de- 
fined within a few centimeters for most event cate- 
gories. The length of the longest track (muon) can be 
measured at 0(10 cm) accuracy or even better. The 
length scale of the electron shower is found at lower 
accuracy. 



The angular resolution for a long muon track is bet- 
ter than one degree. For electron showers it is a few 
degrees. For the additional tracks (proton, pions, gam- 
mas) the angular resolution is weaker, typically some 
tens of degrees at one meter track lengths, while for 
tracks shorter than 0(50 cm) the directions remain 
undefined. 

Successive studies point to 0(1%) deviations be- 
tween the energies of the reconstructed and the true 
events. However, the error varies remarkably with 
different event topologies and the error distribution 
seems to be far from Gaussian. 

For most applications we may assume 5 % energy 
resolution throughout the regime from 1 to 5 GeV. 
This value includes already the nuclear physics un- 
certainties of 1-2%. For events where the neutrino 
direction is known — i.e. known neutrino beams - 
additional kincmatical information is available to im- 
prove the energy resolution. 

The recognition of neutral current background (y + 
X —> v + X* + 7r) has not been fully demonstrated. 
It is evident that a substantial fraction of neutral cur- 
rent events could be confused with v e elastic scattering 
events (ir° — > 77 as v e , n~ as 1^). 

The high-energy performance depends primarily on 
the time resolution of the detector. This suggest to use 
a fast scintillator with short fluorescence times, and 
phototubes of low time jitter. Multi-particle track- 
ing will also profit from pulse shape information for 
individual channels, suggesting the use of FADCs for 
recording the signals of individual PMTs or PMT ar- 
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rays (Sec. 2.5 This is particularly important for event 
recognition and background rejection. 

3.3.4 Conclusions 

Based on the algorithms presented above, energy 
and momentum reconstruction for GeV neutrinos in 
LENA seem well feasible. However, several other as- 
pects must be investigated to obtain a definite result 
on the final sensitivity of LENA for a long-baseline 



experiment. The imminent next step is for sure the 
investigation of neutral current and chargcd-current 
pion backgrounds. Nevertheless, the idea of determin- 
ing the direction of a particle track at degree accu- 
racy in liquid scintillator would have seemed outra- 
geous even a few years ago, let alone the reconstruc- 
tion of multiparticle vertixes. Once again, the ver- 
satility of the liquid scintillation technique has been 
demonstrated. 
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4 Astrophysical neutrino sources 



LENA's core science program is in the low-energy 
range with neutrino energies up to a few tens of McV. 
Most of the relevant sources are based on nuclear reac- 
tions defining this energy scale, in particular the Sun, 
Earth, and power reactors. The same energy range 
is covered by the quasi-thermal emission of neutrinos 
by collapsing stars. The science goals reach from a 
better understanding of astrophysical sources and the 
Earth to the investigation of neutrino properties based 
on flavor oscillations. We begin with neutrinos emit- 
ted from extraterrestrial sources: Core-collapse super- 
novae (SNe) create intense bursts of neutrinos, and in 
particular the next galactic SN will provide more than 
ten thousand neutrino interactions in LENA (Sec. |4.l[ ). 
In addition, the detector will be sensitive to the faint 
signal generated by the diffuse flux from all past SNe 
(Sec. 4.2). Also neutrinos from thermonuclear fusion 
reactions within the Sun will provide a high-statistics 



signal (Sec. 4.3). Finally, the annihilation signature of 



dark-matter particles is studied in Sec. 4.4 



4.1 Galactic Supernova neutrinos 

Measuring neutrinos from the next galactic super- 
nova (SN) is at the frontier of low-energy neutrino 
physics and astrophysics. LENA provides a high- 
statistics neutrino signal — roughly twice that of Super- 
Kamiokande — that can confirm, refute or extend the 
standard paradigm of stellar core collapse and deter- 
mine detailed neutrino "light curves" and spectra. Ad- 
ditionally, LENA's superior energy resolution and var- 
ious flavor-sensitive detection channels are particularly 
advantageous for identifying flavor oscillation effects 
that are sensitive to the unknown mixing angle #13 
and the neutrino mass hierarchy. 

4.1.1 Basic picture 

Core-collapse SNe are the spectacular outcome of 
the violent deaths of massive stars, including the spec- 
tral types II, lb and Ic [55J [53] ■ The early universe 
aside, it is only here that neutrinos do not stream 
freely in spite of their weak interactions and actu- 
ally dominate the dynamics and energetics. The basic 
picture of core collapse is supported by the neutrino 
observation from SN 1987A [5UH55] . This historical 
measurement and the early solar neutrino observations 
remain the only astrophysical sources detected in neu- 
trinos. A high-statistics neutrino observation of stellar 
core collapse is at the frontier of low-energy neutrino 
astronomy, providing an unprecedented wealth of as- 
trophysical and particle-physics information [66H70) . 

A core collapse anywhere in the Milky Way and its 
satellites (such as the Magellanic Clouds) provides a 
detailed neutrino light curve and spectrum. The dis- 
tance distribution is rather broad with an average of 
around 10 kpc [TTj. At this distance, a SN produces 



around 10 4 events in LENA from the dominant inverse 
beta decay reaction v e +P n+e + . Many existing and 
near-future detectors will pick up tens to hundreds of 
events [7HE3], whereas statistics comparable to LENA 
is provided only by Super-Kamiokande. Moreover, the 
high-energy neutrino telescope IceCube at the South 
Pole will register roughly 10 6 uncorrelated Cherenkov 
photons in excess of background, providing superior 
sensitivity to fast signal variations that are suggested 
by recent multi-dimensional simulations [741476] . 

Galactic SNe occur a few times per century as im- 
plied by SN statistics of external galaxies [77H75] . the 
historical record [80, 81 , and the galactic abundance 
of the unstable isotope 26 Al measured with the INTE- 
GRAL gamma-ray observatory [82 . The low-energy 
neutrino sky has been systematically watched since 30 
June 1980 when the Baksan Scintillator Telescope took 
up operation. Only SN 1987A was detected over thirty 
years, beginning to provide non-trivial constraints on 
hypothetical "invisible" core-collapse phenomena [53] . 
Still, the neutrinos from about a thousand galactic 
SNe are on their way and observing one of them is a 
once- in-a- lifetime opportunity. 

Readiness for a galactic SN burst is an essential de- 
tector capability. Reaching Andromeda (M31) and 
Triangulum (M33) at 750 kpc, the next large galaxies 
in the local group, requires megaton-class detectors for 
tens of events. Multi-megaton detectors would detect 
a few neutrinos from SNe out to a few Mpc [84j. One 
could systematically build up an average SN neutrino 
spectrum, but such a project is for the more distant 
future. On a shorter term, another realistic opportu- 
nity to detect SN neutrinos is the diffuse SN neutrino 
background (DSNB) from all past SNe (Sec.[4~2|. 



LENA has about twice the signal statistics of the 
Super-Kamiokande water Cherenkov detector. More 
importantly, it has superior energy resolution, a lower 
threshold, and distinguishes inverse beta decay from 
other channels by recognizing the final-state neu- 
trons. (Dissolved gadolinium in water Cherenkov de- 
tectors [85], currently studied in the EGADS project 
at Super-Kamiokande |86j . will also provide neutron 
tagging.) LENA's excellent energy resolution is a huge 
advantage for recognizing Earth effects in SN neutrino 
flavor oscillations (Sec. 4.1.6). Moreover, LENA is 
complementary to water Cherenkov detectors by in- 
cluding 12 C as a target nucleus and by sensitivity to 
elastic scattering on free protons [87] . 

An increased neutron rate in LENA can signify neu- 
trino emission by thermal processes in the progenitor 
during its last weeks of pre-SN evolution [55J ■ While 
this effect requires the star to be close, the red su- 
pergiant Betelgeuse at about 200 pc [89] is a possible 
candidate. The neutrino burst after collapse would 
trigger about 10 7 events in LENA. The data acquisi- 
tion system must be able to handle such a case without 
being blinded by neutrinos. 
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4.1.2 Supernova astrophysics 

Supernovae and the related, though much rarer, 
long cosmic gamma-ray bursts are the strongest astro- 
physical sources of low-energy neutrinos. These core- 
collapse events are the final stages of the evolution of 
massive stars and as such play a central role in stellar 
and nuclear astrophysics [90]. Besides being the birth 
sites of neutron stars and stellar-mass black holes, they 
are probably the origin of about half of the chemical 
elements heavier than iron. 

While the basic concept of stellar core collapse and 
neutron-star formation was confirmed by the histori- 
cal measurement of neutrinos from SN 1987A, our un- 
derstanding of the detailed processes driving the core's 
evolution and ultimately causing the SN blast remains 
incomplete and has little empirical underpinning |59j . 
Observations of the bright electromagnetic spectacle 
that accompanies stellar death provide only indirect 
information about the initiating mechanism: the cen- 
ter of the explosion is obscured by several solar masses 
of intransparent, gaseous ejecta. 

The current theory of stellar explosions strongly re- 
lies on numerical modeling that requires empirical sup- 
port. While the produced heavy elements somewhat 
probe the conditions around the origin of the explo- 
sion, only neutrinos and gravitational waves can es- 
cape directly from the densest regions and thus are 
unique messengers from the very center. Their time- 
dependent signal features carry detailed and comple- 
mentary information of the evolving thermodynamical 
state and of the dynamical motions in the compact 
remnant assembling at the heart of the dying star. 

Neutrino measurements from a galactic SN together 
with the detection of a gravitational-wave burst will 
trigger a breakthrough in our understanding of some of 
the most important questions in stellar astrophysics: 
What are the conditions in collapsing cores of mas- 
sive stars? Is there a significant amount of rotation? 
Do magnetic fields play an important role? How can 
massive stars succeed to reverse their catastrophic in- 
fall to a powerful explosion? What are the properties 
of hot nuclear matter? What are the mass, radius, 
and binding energy of the newly- formed neutron star? 
Does the compact remnant undergo a phase transi- 
tion to a more compressed quark-matter state or even 
collapse to a black hole? Are SN explosions and new- 
born neutron stars the long-sought formation sites of 
the heaviest neutron-rich elements that arc made by 
the rapid-neutron capture process (r-process)? 

4.1.3 Expected neutrino signal 

In spite of many open questions, today's numerical 
SN models may well provide a reasonable first guess 
of the signal characteristics. Spherically symmetric 
simulations have recently provided robust explosions 
for small progenitor masses of 8-10 Mq, the class of 
electron-capture SNe (or O-Mg-Ne-core SNe) [9"2"H9T] . 
For more massive stars, leading to the conventional 



iron-core SNe, strong deviations from spherical sym- 
metry caused by large-scale convection and the stand- 
ing accretion shock instability (SASI) are probably 
important, but full-fledged 3D simulations with suf- 
ficiently sophisticated neutrino transport are only be- 
ginning to come into reach. 

The expected neutrino signal consists of three main 
phases (Fig. 16 1, testing different aspects of SN theory 
and neutrino flavor oscillations. 

1. Few tens of ms after bounce: Shock break-out 
and deleptonization of the outer core layers, 
emission of the "prompt v e burst." Emission 
of other flavors only begins and that of v e is at 
first suppressed. Largely independent of progen- 
itor mass and equation of state [55] ■ 

2. Accretion phase, few tens to several hundred ms, 
depending on progenitor mass and other param- 
eters. Shock stalls at 100-200 km, neutrino emis- 
sion is powered by infalling material. Fluxes of 
v e and v e much larger (a factor of two is not 
unrealistic) than those of the other flavors. Pro- 
nounced hierarchy (E Ve ) < (E^J < (E Ulc ) with 
v x representing any of v^ T and v^r- Large-scale 
convection and SASI mode build up, leading to 
strong time variation of the neutrino signal. 

3. Cooling phase, up to 10-20 s. Neutrino flux 
powered by cooling of the deep core on a 
diffusion time scale. Approximate luminosity 
equipartition between all species and only a mild 
(Ep s } I (E Vx ) hierarchy. Larger v e number flux 
due to de-leptonization. 

Of course, completely different signatures can arise if 
new phenomena occur. Examples are a late time QCD 
phase transition or black hole formation. 

A broad range of possible spectral properties of the 
neutrino signal have been considered in the literature, 
but until recently the only numerical model of multi- 
flavor SN neutrino emission from bounce to long-term 
cooling was provided by the Livermore group |99j . It 
is only during the past year that their pioneering work 
has been superseded by modern long-term simulations. 
For the first time hydrodynamic simulations coupled 
with modern neutrino Boltzmann solvers in ID have 
been carried all the way to proto-neutron star cool- 
ing. The Basel group has evolved progenitors with 
different masses up to 10 s after bounce ]§T\. The 
Garching group has published a similar simulation for 
an electron-capture SN [93] . 

The emerging picture suggests smaller average en- 
ergies than often assumed and much less pronounced 
spectral hierarchies, particularly during the cooling 
phase. Time- integrated values in the range (E^ c ) — 
12-14 MeV, and somewhat larger for Vu,ti l°°k rea ~ 
sonable and are in agreement with the SN 1987A ob- 
servations and with analytic [100] and Monte-Carlo 
studies of neutrino transport [101] , We use such rel- 
atively modest energies to gauge our expectations for 
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FIG. 16: Neutrino signal of a core-collapse SN for a 10.8 M© progenitor according to a numerical simulation of the Basel 
group All quantities are in the laboratory frame of a distant observer. In this spherically symmetric simulation the 
explosion was triggered artificially by implementing enhanced neutrino energy deposition. Left: Prompt neutrino burst. 
Middle: Accretion phase. Right: Cooling phase. 



LENA. Of course, it is the very purpose of SN neutrino 
observations to measure the neutrino flux character- 
istics independently of theoretical predictions and it 
remains quite possible that typical SNe produce much 
larger or very different signals. Moreover, one expects 
large variations between different SNe, depending for 
example on different rates and amounts of accretion. 

4.1.4 Detection channels in LENA 

The purpose of a high-statistics SN neutrino obser- 
vation is to measure time-dependent features. How- 
ever, a first impression of the detector capabilities is 
gained from integrated detection rates. Detailed spec- 
tral studies will be important and so we assume a 
range of different source characteristics. To this end 
we treat the SN schematically as a black-body source 
for all neutrino species. We assume a total emitted 
energy of -Etot = 3 x 10 53 erg, equipartitioned among 
all neutrino species, and Maxwell-Boltzmann spectra 
with (E v ) = 12, 14 and 16 MeV. Of course, in a real- 
istic SN one expects flavor-dependent differences that 
will be used, for example, to search for flavor oscilla- 
tions. 

LENA's golden detection channel is inverse beta de- 
cay v e + p — > n + e + . The produced neutron thermal- 
izes and wanders in the detector until it is captured 
by a proton, n+p —> d + -f (2.2 MeV) after an average 
time of ^250 us. The large homogeneous detection vol- 
ume ensures efficient neutron capture and 7 detection. 
Therefore, these events are tagged by the delayed co- 
incidence between the prompt positron and the 7-ray 
from neutron capture. 

Three charged-current (CC) reactions measure v e 



and v e fluxes and spectra while three neutral-current 
(NC) processes, sensitive to all flavors, give informa- 
tion on the total flux. Typical event rates for a generic 



SN at a distance of 10 kpc are reported in Table VIII 
for our representative cases of (E v ). A LAB-based 
scintillator and a fiducial mass of 44 kt is assumed, 
providing about 3.3 x 10 33 protons. The 12 C reactions 
have a high kinematical threshold (E > 15 MeV). The 
resulting steep energy dependence in principle pro- 
vides information on the neutrino spectra. 

It is particularly difficult to detect the Vu, v T , 
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TABLE VIII: Expected event rate in LENA for a SN at 
a distance of 10 kpc, where v stands for a neutrino or 
antineutrino of any flavor. The NC rates are summed over 
all flavor channels. Our three representative values for {E v ) 
are taken to be equal for all flavors. The Birks constant 
was taken to be 0.010 cm/MeV, and a threshold of 0.2 MeV 
was assumed to calculate total number of events for the 
vp — > pv channel. 
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and v T flavors and measure their energies, because, 
unlike v e and D e , they have only NC interactions. On 
the other hand, observing these flavors is essential to 
disentangling flavor mixing and correctly estimating 
the total energy emitted in neutrinos. Directly ob- 
serving two spectral components due to flavor mixing 
in the CC data, one for v e and another for the non 
electronic flavors, will only be possible if the average 
energies are significantly different. That possibility is 
not strongly favored by current theory. Promisingly, 
one of the main strengths of LENA is its low threshold 
that should allow us to observe neutrino-proton elas- 
tic scattering events, which have spectral information 
and a substantial yield from these neutrinos [S3 U02| . 
Other NC channels typically lack spectral informa- 
tion and have lower yields. The expected spectrum 
of elastic neutrino-proton scattering events for typi- 
cal SN parameters is shown in Fig. [17] The spectrum 
has been smeared with the energy resolution given in 
the previous section, and the threshold is taken to be 
0.2 MeV to stay above the 14 C background. With a 
few thousand observed events, one could reconstruct 
the non electron flavor neutrino spectra with almost 
the same precision as that of v e |102j . The prompt v e 
burst alone will produce around 50 events, depending 
on the mixing scenario, by electron scattering and 90, 
independent of oscillations, by proton elastic scatter- 
ing (using the Basel model of Fig. 16). 



4.1.5 Astrophysical lessons 

A gravitational-wave signal provides information on 
non-radial deformation and non-spherical mass mo- 
tions [103] . whereas a high-statistics neutrino signal 
allows us to follow directly the different stages of 
core collapse without additional assumptions (Fig. 16 1 . 
The prompt v e burst is a robust and uniform land- 



mark structure of all theoretical predictions. Because 
of LENA's capability of distinguishing NC and CC 
events, it offers a unique possibility of identifying this 
feature. For example, one could estimate the SN dis- 
tance in the plausible case that the optical display 
is hidden behind the dense gas and dust clouds of 
a star- forming region [98] . Moreover, one could use 
the prompt v e burst in LENA for coincidence mea- 
surements with the gravitational wave burst that may 
arise at core bounce. Using the prompt v e burst 
could provide an even sharper coincidence than can 
be achieved with the onset of the v e signal in Super- 
Kamiokande [104] and IceCube [105] . Moreover, the 
prompt v e burst could help to find the SN direction by 
neutrino triangulation [106] . although the recoil elec- 
tron signal in a water Chercnkov detector provides su- 
perior pointing capabilities [1061 1107] . 

The magnitude of the z/ e and v e accretion luminosi- 
ties after core bounce (Fig. 16 middle) depends on the 



mass infall rate and thus on the progenitor-dependent 
structure of the stellar core, with more massive cores 
producing higher luminosities |108| 1109] . Luminos- 
ity variations during this phase [74H76] . accompanied 
by sizable gravitational-wave emission at several hun- 
dred Hz [7411110] would confirm the presence of violent 
hydrodynamic instabilities stirring the accretion flow 
around the assembling neutron star. Such activity and 
a several hundred millisecond delay of the onset of the 
explosion are expected within the framework of the de- 
layed neutrino-driven mechanism. A pronounced drop 
of the v e and v e luminosities, followed by a close simi- 
larity to those of heavy-lepton neutrinos, would finally 
signal the end of the accretion phase and the launch 
of the outgoing SN blast wave. The cooling signature 
of a nascent neutron star is characterized by a mono- 
tonic and gradual decline of the neutrino emission. It 
would be prolonged if additional energy was released 
by phase transitions in the nuclear matter. Exotic sce- 
narios might feature a secondary v e burst [111] or an 
abrupt end of neutrino emission if the collapse to a 
black hole occurcd 112 . 



LENA could provide even more information: Due 
to its superior energy resolution it could help to dis- 
entangle source-imposed spectral features from those 
caused by neutrino-flavor conversions. Moreover, de- 
tecting significant numbers not only of v e but also 
of v e and heavy-lepton neutrinos (Table VIII) would 



yield at least time-averaged spectral information for 
different emission channels. Conceivably one could ex- 
tract information on the neutron-to-proton ratio in the 
neutrino-processed SN outflows, presently also a sen- 
sitive result of numerical modeling of a multitude of 
complex processes. The relative abundance of neu- 
trons and protons determines the conditions for nu- 
cleosynthesis and are set by competing v e and v e cap- 
tures, which in turn depend delicately on the relative 
fluxes and spectral distributions of these neutrinos. A 
LENA measurement of a SN burst may offer the only 
direct empirical test of the possibility for r-processing 
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in the SN core, except for an extremely challenging 
in-situ measurement of r-process nuclei in fresh SN 
ejecta. 

4.1.6 Particle physics and neutrino properties 

On the particle-physics side, the high-statistics ob- 
servation of a SN neutrino burst can provide cru- 
cial particle-physics lessons. Numerous results derived 
from the sparse SN 1987A data can be refined. One 
can also probe more exotic scenarios. Spin-flavor con- 
versions caused by the combined action of magnetic 
fields and matter effects can transform some of the 
prompt v e burst to u e , leading to a huge inverse- 
beta signal |113) . Such an observation would provide 
smoking-gun evidence for neutrino transition magnetic 
moments. Non-radiative decays would also produce a 
v e — » v e conversion during the prompt burst |114j . 

Perhaps of greatest interest are flavor oscillations. 
Neutrinos propagating through the SN mantle and en- 
velope encounter a large range of matter densities, 
allowing for Mikheyev-Smirnov-Wolfenstein (MSW) 
conversions driven by the atmospheric neutrino mass 
difference and the small mixing angle #13. Therefore, 
in principle a SN neutrino signal is sensitive to the two 
as yet unknown neutrino mixing parameters: #13 and 
the ordering of the neutrino masses that could be in 
the normal (NH) or inverted hierarchy (IH). 

Our understanding of SN neutrino oscillations has 
recently undergone a change of paradigm by the in- 
sight that the neutrino-neutrino refractive effect is cru- 
cial. These collective (or self-induced) flavor conver- 
sions occur within a few hundred km above the neu- 
trino sphere; see Ref. 115 for a review of the recent 
torrent of literature on this topic. The most important 
observational consequence is a swap of the v e and v e 
spectrum with that of v x and v x in certain energy in- 
tervals |116j . The sharp spectral features at the edges 
of these swap intervals are known as "spectral splits." 
Their development depends on the neutrino mass hi- 
erarchy as well as on the ordering of the flavor fluxes 
at the source. Therefore, the split features can depend 
on time in interesting ways |117H121] . 

The main problem to detect oscillation features is 
that one can not rely on detailed theoretical pre- 
dictions of the flavor-dependent fluxes and spectra. 
Therefore, model-independent signatures are crucial. 
One case in point is the energy-dependent modulation 
of the neutrino survival probability caused by Earth 
matter effects that occur if SN neutrinos arrive at the 
detector "from below" |122j . The appearance of Earth 
effects depends on the flux and mixing scenario (121] - 
Therefore, its detection could give hints about the pri- 
mary SN neutrino fluxes, as well as on the neutrino 
mass hierarchy and the mixing angle #13. 

The excellent energy resolution of LENA is a par- 
ticular bonus for discovering small energy-dependent 
flux modulations caused by Earth effects [123] . but 
of course depends on seeing the SN shadowed by 



the Earth. In a far-northern location such as the 
Pyhasalmi mine in Finland the shadowing probability 
for a galactic SN is about 58%, against an average of 
50% for a random location [71]. A particularly inter- 
esting scenario consists of a large volume scintillator 
detector in the north to measure the geo-neutrino flux 
in a continental location and another one in Hawaii 
to measure it from the oceanic crust. The probabil- 
ity that only one of them is shadowed exceeds 50% 
whereas the probability that at least one is shadowed 
is about 80%. Therefore, Pyhasalmi and Hawaii are 
complementary both for observing geo-neutrinos and 
Earth matter effects in SN neutrinos. 

Additional signatures of flavor conversions can be 
imprinted by matter effects of the shock fronts in the 
SN envelope |124j . The number of events, average en- 
ergy, or the width of the spectrum may display dips 
or peaks for short time intervals [125] 1126] . Such sig- 
natures yield valuable information about shock-wave 
propagation, the neutrino mass hierarchy and #13. 
However, realistic chances to detect shock features re- 
main unclear. The flavor-dependent spectral differ- 
ences in the anti-neutrino channel are probably small 
during the cooling phase. Moreover, strong turbu- 
lence in the post-shock regions could affect these sig- 
natures |127| . 

4.1.7 Summary 

A worldwide network of neutrino and gravitational- 
wave detectors, constituting the SuperNova Early 
Warning System (SNEWS) [128], will provide early 
warning and detailed multi-messenger measurements 
of the next nearby SN. A high-statistics neutrino ob- 
servation, even from a single SN, will go a long way 
to answering many fundamental questions about the 
role of neutrinos for the astrophysics of core collapse 
and may shed new light on the properties of neutrinos 
and other particles. LENA will play an exceptional 
role due to its low energy threshold, excellent energy 
resolution, and multi-channel signatures that will al- 
low one to disentangle flavor-dependent properties of 
the neutrino signal and to identify subtle modulations 
imprinted by Earth effects. LENA may be the only 
facility that is able to spot the prompt v e burst and 
thus the earliest and largely model-independent signa- 
ture of stellar death, even yielding an estimate of the 
SN distance. 

4.2 Diffuse Supernova neutrinos 

The diffuse SN neutrino background (DSNB) from 
all core collapse events in the universe is a guaran- 
teed neutrino flux from cosmological distances. DSNB 
v e can be detected at energies above 10 MeV, where 
the reactor neutrino background vanishes and atmo- 
spheric neutrino backgrounds are small and likely con- 
trollable. LENA provides about twice the counting 
rate of Super-Kamiokande, and together they could 
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collect 5-10 events per year. Measurement of the av- 
erage D e emission spectrum will help test models of 
SNe, variation in emission, and neutrino properties. 

4.2.1 Basic picture 

A great and varied scientific return is expected from 
the observation of a nearby SN (Sec. 4.1), but such 
events are rare in the Milky Way. The guaranteed 
DSNB flux provides a way to detect SN neutrinos 
without a fortuitous burst |129H147j . DSNB signals 
depend on three ingredients. First, the cosmic core 
collapse rate, about 10 per second in the causal hori- 
zon; this is determined by astronomical measurements 
that are already precise and quickly improving. Sec- 
ond, the average SN neutrino emission, which is ex- 
pected to be comparable for all core collapses, in- 
cluding those that fail and produce black holes (for 
which it may be even larger, as discussed below); 
this is the quantity of fundamental interest. Third, 
the detector capabilities, including the energy depen- 
dence of the cross section and detector backgrounds; 
Super-Kamiokande and LENA should be able to de- 
tect DSNB v e . 

Detecting the DSNB is important even if a Milky 
Way burst is observed. DSNB v e will provide a unique 
measurement of the average neutrino emission spec- 
trum to test SN simulations. Comparison to data from 
SN 1987A and an eventual Milky Way SN will test the 
variation between core collapses. While the statistics 
of DSNB events will be low, like those of SN 1987A, 
this data will more effectively measure the exponen- 
tially falling tail of the spectrum at high energies. The 
DSNB is also a new probe of stellar birth and death: 
its energy density is comparable to that of photons 
produced by stars, but the DSNB is unobscured and 
has no known competition from astrophysical sources. 
Finally, the DSNB data will test flavor mixing and 
more exotic particle properties. 

The importance of running LENA to detect the 
DSNB should not be underestimated. If Super- 
Kamiokande does not add gadolinium, or does but 
encounters technical problems, LENA could be the 
only experiment to detect the DSNB. If both exper- 
iments are successful, their data, based on detection 
by inverse beta decay above 10 MeV, would be sim- 
ilar. Having two independent experiments would be 
very valuable, as this will be a challenging measure- 
ment. In addition, LENA provides about twice the 
counting rate of Super-Kamiokande. Collecting statis- 
tics at a combined three times higher rate than Super- 
Kamiokande alone could have a decisive impact on the 
physics that could be extracted. 

4.2.2 DSNB signals 

The DSNB event rate spectrum follows from a line 
of sight integral for the radiation intensity from a dis- 
tribution of distant sources. After integrating over all 
angles due to the isotropy of the DSNB and the trans- 
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On the right hand side, the ingredients are ordered 
as described above. The first is the comoving cosmic 
core-collapse rate, in units Mpc -3 yr _1 ; it evolves with 
redshift. The second is the average time-integrated 
emission per SN, in units MeV -1 ; redshift reduces 
emitted energies and compresses spectra. The third 
is the number of targets times the detection cross 
section; this does not need to be under the inte- 
gral. The last term is the differential distance, where 
\dtfdz\~ 1 = H (l + z)[Q A + ^m(l + zf] 1 / 2 ; the cos- 
mological parameters are taken as H — 70 km s _1 
Mpc -1 , J7a = 0.7, and fi m = 0.3. (The cosmology and 
the SN rate derived from star formation rate data are 
really one combined factor proportional to the ratio 
of the average luminosity per galaxy in SN neutrinos 
relative to stellar photons.) The left hand side is the 
DSNB spectrum in visible energy _E v is; the relation to 
neutrino energy E v depends on cross section and de- 
tector specifics. We next consider details of the three 
main ingredients. 

Cosmic SN rate. The cosmic core collapse rate is 
precisely known [14814150] . The redshift range relevant 
for the DSNB depends on energy, with lower energies 
probing higher redshifts. For detected energies above 
10 MeV, most DSNB neutrinos are emitted at red- 
shifts z < 1, where the astronomical data are most 
precise. The best determinations of the core collapse 
rate come from predictions based on measured star 
formation rates and related observables such as the ex- 
tragalactic background light |149) . As massive stars 
are short-lived, the redshift evolution of the core col- 
lapse and star formation rates must be the same. The 
relative normalization depends on just the minimum 
mass for core collapse, about 8M |151| : the predicted 
rate depends only weakly on the assumed stellar initial 
mass function because star formation data primarily 
sample massive stars. Comparable neutrino fluxes are 
expected for ordinary SNe and those that are faint, 
obscured, or even failed [92" 1 [T52HT54] . so the DSNB 
does not depend much on the outcomes, though it may 
be larger than assumed here. Measured SN and pre- 
dicted core collapse rates are in reasonable agreement, 
and the data will quickly improve |150l 11551 1156] . 

The local core collapse rate is Rsn(z = 0) = 
(1.25 ± 0.25) x 10~ 4 Mpc" 3 yr" 1 [149J. The evolu- 
tion of the comoving rate, roughly the rate per galaxy, 
has a strong and clear rise of one order of magnitude 
between z — and z = 1 and then a slow and eventu- 
ally steepening decline at higher redshift. Taking into 
account the variation of the uncertainties with red- 
shift, the uncertainty on the DSNB due to that on the 
core collapse rate is presently ±40% [149] . This will 
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decrease quickly with new data, so that the focus of 
DSNB measurements will be on the neutrino emission 
parameters [157] . 

SN neutrino emission. While we have some in- 
formation about neutrino emission from SN 1987A 
and SN simulations, detecting the DSNB is neces- 
sary to measure the average emission per SN. It is 
typically assumed that the total energy in neutrinos 
is 3 x 10 53 erg, that each flavor carries 1/6 of this, 
and that the spectra are quasi-thermal with tempera- 
tures of several MeV. But the total energy, its partition 
among flavors, and spectral distributions and average 
energies may be different or show more variation than 
expected. Uncertainties include those due to the col- 
lapse mechanism 158-164 , the effects of progenitor 
mass, rotation, and magnetic fields |165H167j . the neu- 
tron star equation of state |168H171j . and effects due 
to neutrino properties (Sec. 
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The neutrino emission is parameterized here with 
a Maxwcll-Boltzmann thermal spectrum, tp{E u ) = 
Etat [E 2 /(2T 3 )] exp(-E u /T v ), where the total energy 
and temperature (average emitted energy (E v ) = 3T„) 
are for v e after neutrino flavor oscillations, which occur 
in the SN and not en route. Following Sec. 4.1. 3| the 
nominal expectation for v e might be Btot = 0.5 x 10 53 
erg and T v = 4 MeV, with large uncertainties. Fur- 
ther, we do not know if SN 1987A was a typical core 
collapse or if present SN simulations are correct. Mea- 
surements of the DSNB are needed to help decide. 

The DSNB signal may be larger than assumed here, 
due to unusual core collapse outcomes that could be 
disproportionately important due to their larger-than- 
average neutrino emission. The most interesting pos- 
sibility is prompt black hole formation, as this is ex- 
pected to have a nonzero rate even in standard scenar- 
ios |1721I173| . and present constraints allow even larger 
rates QH] [1501 HSU [HO Q23J. Even though the neu- 
trino emission can be cut off, it is expected to be en- 
hanced before that, such that the time- integrated total 
and average neutrino energies can be larger than usual 
[9"21 I152H154] . Other possibilities include emission 
from the hot, magnetized corona of a proto- neutron 
star or accretion disk [175] or from fallback |176j . The 
DSNB will thus be especially valuable for probing out- 
comes that may not occur for a Milky Way core col- 
lapse and the corresponding extreme physical condi- 
tions in such collapses [1561 1177] . 

Detector capabilities. The detection channel in 
LENA and Super-Kamiokande is inverse beta decay, 
v e + p — > e + + n, while other DSNB neutrino in- 
teractions have smaller detectable rates |178j . The 
positron kinetic energy is close to that of the neutrino, 
T e ~ E v — 1.8 MeV, with a nearly isotropic distribu- 
tion. The low-energy neutron will thermalize and then 
register its presence by radiative capture. The time 
and space coincidence between positron and neutron 
suppresses detector backgrounds. 

LENA has 2.9 x 10 33 free protons in 44 kt of scin- 
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FIG. 18: DSNB signal spectra in LENA, with labeled lines 
indicating assumed u e emission spectra, which are of the 
Maxwell-Boltzmann form with the indicated temperatures 
and which have a total energy of 0.5 x 10 52 erg. Below 
10 MeV, the reactor neutrino background (spectrum not 
shown) is overwhelming. At higher energies, atmospheric 
neutrino backgrounds should be small and controllable. 



tillator (Super-Kamiokande has 1.5 x 10 33 in 22.5 kt 
of water). The cross section is <t(E v ) ~ 9.42 x 
10~ 44 cm 2 (Ev/MeV— 1.3) 2 at lowest order; we use the 
corrections to the cross section and kinematics from 
Refs. |179]I180] . In LENA, the visible positron energy 
is its kinetic energy plus the annihilation energy with 
an electron, E v - ls = T e + 2m e c 2 (Super-Kamiokande 
defines visible energies via the positron total energy, 
E V i S — T e +m e c 2 ). The effects of energy resolution on 
the DSNB spectrum are negligible. 

Figure [18] shows the expected DSNB signal spec- 
trum in LENA, following the above details and as- 
suming perfect detection efficiency. A range of SN 
i> e emission spectra, parameterized by changes in T v , 
are shown. There are also uncertainties due to E to t 
and the assumed spectrum shape. (Astrophysical 
uncertainties — those due to the SN rate alone — are ne- 
glected because they are small and quickly decreasing.) 
The minimal allowed case is close to Etot = 0.5 x 10 53 
erg, T = 4 MeV, as this DSNB prediction is com- 
parable to that obtained from a direct non-parametric 
reconstruction of the high energy SN 1987A data [142J 
(see also Ref. [1371 1141] for previous analyses of 
SN 1987A data applied to the DSNB). Maximal al- 
lowed cases (not shown), based on just the 2003 Super- 
Kamiokande limit, are close to E to t ~ 1 X 10 53 erg, 
T v = 6 MeV, or E tot - 2 x 10 53 erg, T v = 4 MeV [133] . 

The v e emission parameters can thus be directly 
measured from the DSNB spectrum if the atmospheric 
neutrino backgrounds can be controlled. The spectra 



in Fig. 18 contain 70, 55, and 35 events in the range 
10-30 MeV for ten years of LENA running, corre- 
sponding to a statistical uncertainty of 12-17%. Data 
from Super-Kamiokande will also help, especially if it 
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begins running with gadolinium soon. The spectrum 
shape will help break the degeneracy between E tot and 
T v . In effect, E tot will be probed best by the lower en- 
ergy data, where the T v dependence is weakest, and T v 
will be probed best by the higher energy data, where 
the falloff is exponential. Precise data will also test 
new physics scenarios |181[ I182j . 

Theory and data from a future Milky Way SN will 
be needed to relate these v e emission parameters with 
flavor mixing effects included to other SN parameters 
such as the total energy emitted by all flavors and the 
spectra before neutrino flavor mixing. 

4.2.3 LENA detector backgrounds 

As both positron and neutron emerging from an in- 
verse beta decay are observable in liquid scintillator, 
single-event backgrounds can be suppressed very ef- 
fectively. The only remaining backgrounds are due to 
other sources of & ei namely reactor and atmospheric 
u e that are irreducible. They limit the DSNB detec- 
tion window to 10-30 MeV, the exact range depend- 
ing on the detector site [9]. Other backgrounds are 
signals mimicking the fast coincidence signal: Cosmo- 
genic /3n-emitting isotopes, fast neutrons from the sur- 
rounding rocks, and NC interactions of higher-energy 
atmospheric neutrinos. However, these backgrounds 
can be identified, either by their production, location 
or pulse shape. These reducible backgrounds, strate- 
gies for their rejection and the accompanying loss in 
DSNB detection efficiency are outlined below. 

The delayed coincidence signal of inverse beta decay 
enables LENA to easily reject what is the predominant 
backgrounds for DSNB detection in water Cherenkov 
detectors, i.e. low energy muons produced in CC re- 
actions of atmospheric neutrinos, solar neutrinos and 
spallation products of cosmic muons. 

Cosmogenic /3n-emitters. are unstable isotopes 
produced by cosmic muons crossing the detector. 
They mimic the u e coincidence by the prompt emission 
of the electron followed by the emission of a neutron 
from an excited state of the daughter nucleus. For- 
tunately, only 9 Li (T1/2 — 178 ms) has a large enough 
Q- value to add to the background. If no cuts are ap- 
plied, the 9 Li rate is of the same order of magnitude 
as the DSNB signal. Due to its short half-life, 9 Li is 
easily associated to its parent muon: It is sufficient 
to veto a cylinder with 2 m radius around each muon 
track for 1 s (^5 7i/ 2 ), to reduce the residual 9 Li rate 
to about 2%, while the introduced dead time corre- 
sponds only to ~ 0.1% of the total measurement time 

Fast neutrons, are produced in the surrounding rock 
by cosmic muons that pass the detector undetected. 
There is a chance for these neutrons to propagate into 
the target volume. The coincidence is mimicked by 
a prompt signal due to elastic scattering of protons, 
while the delayed signal is caused by the capture of 
the stopped neutron. The fast neutron background 



rate was analyzed by Monte-Carlo simulations. As 
most neutrons will stop at the verge of the scintillator, 
a fiducial volume cut greatly reduces the rate. For 
10 m fiducial radius, ^0.2 fast neutron events per year 
are expected [SO] . However, proton recoils of neutrons 
feature a different typical pulse shape than positron 
signals. This allows for an alternative approach to 
distinguish neutron events more effectively, reducing 
the fast neutron background to 0.12 events per year in 
the nominal fiducial volume of 12 m radius. 

NC reactions of atmospheric neutrinos, might 
prove to be the most dangerous background: Besides 
the intrinsic background of atmospheric P e 's, atmo- 
spheric neutrinos at higher energies knock out neu- 
trons from 12 C in the scintillator. Neutron scattering 
off protons or particles emitted in the de-excitation 
of the remaining nucleus cause a prompt signal, while 
the neutron is later captured, mimicking the signal 
coincidence. MC simulations point towards a back- 
ground rate about 10-20 times higher than the ex- 
pected DSNB signal 183 . Several strategies have 
been devised to cope with this background: A pos- 
sible way is to search for the delayed /3 + decay of the 
residual 11 C that remains after the neutron knock-out. 
If the 11 C nucleus is created in its ground state, this 
is a very effective strategy, reducing the background 
by a factor of 2. However, if an excited n C state is 
created, it will mostly de-excite via proton, neutron 
and alpha emission. In this case, the only way of dis- 
crimination is a pulse shape analysis of the prompt 
signal. This rejection profits not only from the in- 
trinsic differences in light emission of protons and a- 
particles compared to the z? e -indudced positrons but 
e.g. also from the spatially extended energy deposition 
of high-energetic 7-quanta. The discrimination power 
as well as the remaining DSNB detection efficiency is 
currently evaluated in MC studies. Preliminary results 
indicate that in spite of a painful loss in efficiency, a 
signal-to-background ratio greater than unity can be 
obtained 50 . However, further MC studies as well 
as laboratory measurements on proton quenching in 
liquid scintillator are needed to quantify this result. 

While the latter background is absent in water 
Cherenkov detectors, their inability to detect the de- 
layed neutron signal makes them vulnerable to solar 
neutrinos, the decay of invisible muons and all kinds 
of spallation products. Nevertheless, much improved 
Super-Kamiokande sensitivity to DSNB D e is expected 
by the time LENA comes into operation. Their 2003 
limit is already strong and will improve with further 
data |184j . If gadolinium is added, Super-Kamiokande 
will reject detector backgrounds above 10 MeV and 
will cleanly collect a few DSNB signal events per 
year [85]. 

4.2.4 Summary 

The DSNB is a very promising astrophysical neu- 
trino source, with at most a factor of a few improve- 
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merit in flux sensitivity required for a first detection. 
This will directly probe the neutrino emission per core 
collapse via the measured v e spectrum above about 
10 MeV. This spectrum averages over all core-collapse 
outcomes, including some which may be relatively rare 
but which may be of disproportionate importance due 
to larger-than-usual neutrino emission. The main ad- 
vantages of LENA are its large size, native ability 
to detect neutrons to tag v e + p — > e + + n events, 
and low detector backgrounds and consequent low en- 
ergy threshold. LENA may make a first detection of 
the DSNB and would significantly increase statistics 
over Super-Kamiokande alone, leading to more deci- 
sive probes of the average neutrino emission per core 
collapse, a key comparison point for SN models and a 
Milky Way SN. 

4.3 Solar neutrinos 

Solar neutrino research is a mature field that has ac- 
cumulated a long series of outstanding achievements. 
Originally conceived as a powerful tool to investigate 
the Sun's deep interior, solar neutrinos provided the 
first indication for neutrino flavor oscillations and thus 
contributed in crucial ways to discover and analyze 
this profound phenomenon. Solar experiments provide 
sensitivity to Amf 2 and especially #12, with the fas- 
cinating prospect of a possible positive indication re- 
garding the value of the subleading # 13 angle, connect- 
ing solar and atmospheric sectors of the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) neutrino mixing ma- 
trix. These achievements are the starting point of the 
LENA solar neutrino program because high-statistics 
measurements will resolve energy spectra and possible 
time variations in unprecedented detail. 

4.3.1 Introduction 

The Standard Solar Model (SSM). The effort to 
develop a model able to reproduce fairly accurately 
the solar physical characteristics, as well as the spec- 
tra and fluxes of the several produced neutrino com- 
ponents, was led for more than forty years by the late 
John Bahcall; this effort culminated in the synthesis 
of the so called Standard Solar Model (SSM) [185], 
which represents a true triumph of the physics of 20 th 
century, leading to extraordinary agreements between 
predictions and observables. 

However, over the last years such previous excel- 
lent agreement has been seriously compromised by the 
downward revision of the solar surface heavy-element 
content from Z = 0.0229 [T5S] to Z = 0.0165 [W] . 
leading to a severe discrepancy between the SSM and 
the helioseismology results. Resolution of this puz- 
zle would imply either to revise the physical inputs of 
SSM or to modify the core abundances. 

In 2009 a complete revision of the solar photospheric 
abundances for nearly all elements have been done 
188 . This revision includes a new three dimensional 



hydrodynamical solar atmosphere model with an im- 
proved radiative transfer and opacities. The obtained 
results give a solar metallicity Z — 0.0178. In |189j the 
three different sets of solar abundances: GS98 186 , 
AGS05 [W] and AGS09 [HE] have been used originat- 
ing two different, low metallicity or high metallicity, 
versions of the solar model. 

Neutrino oscillations and the MSW effect. So- 
lar neutrino oscillations are governed by Ato^ and 
O12 of the PMNS mixing matrix. At low energies, 
v e survival probabilities are described well by vac- 
uum oscillations. However, at energies above ~1 MeV, 
matter effects first pointed out by Mikheyev, Smirnov 
and Wolfenstein (MSW) [1201 EEE] enhance the con- 
version v e — \ f M , r , leading to a further suppression 
of the v e rate observed in terrestrial detectors. By 
now, this MSW-LMA oscillation scenario is well con- 
firmed by solar neutrino experiments for vacuum- and 
matter-dominated regimes. However, the vacuum- 
matter transition region from 1 to ~5 MeV remains to 
be explored and might hold evidence for non-standard 
neutrino physics. 

4.3.2 Experimental status 

For almost 40 years, solar neutrino detectors have 
accumulated a large amount of data. Beyond the con- 
firmation of thermonuclear fusion as the solar energy 
source, the comparison of the experimental results to 
the accurate predictions of the SSM on solar neu- 
trino flux and spectrum led to the establishment of 
the LMA-MSW oscillation scenario in the solar sec- 
tor. 

In the 70's and the 80's the Chlorine radiochemical 
experiment at Homestake [192] and the Cherenkov de- 
tector Kamiokande [193] played the fundamental role 
to establish on solid rock basis what became known as 
the Solar Neutrino Problem (SNP), i.e. the persisting 
discrepancy between the measured and predicted solar 
neutrino flux. 

In the early 90's the scene was dominated by the sec- 
ond generation of radiochemical experiments based on 
Gallium, GALLEX/GNO and SAGE [IMfH^ . which 
not only reinforced the physics case of Homestake and 
Kamiokande, thus further strengthening the SNP, but 
also marked the first detection ever of the overwhelm- 
ing flux of the pp neutrinos. This milestone result 
represented the first direct proof of the nuclear burn- 
ing mechanism as the actual stellar energy generating 
engine. 

Later, the final assault to the SNP was launched by 
the three real time experiments Super-Kamiokande, 
SNO (whose measurements were decisive to solve the 
SNP) and Borexino, which were complemented in their 
effort by the reactor neutrino experiment KamLAND: 

Super-Kamiokande. Designed to detect the 
Cherenkov light emerging from the elastic scattering 
of the incoming neutrinos off the electrons of the wa- 
ter acting as detection medium, Super-Kamiokande 
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started its operation in 1990 at the Kamioka mine in 
Japan. It is a gigantic detector, with its 50 ktons of 
water viewed by more than 10 000 20- inch phototubes. 

Over 15 years of measurement, the experiment re- 
turned consistent data on the 8 B-^ flux and spec- 
trum [T9"T1 HM\ above a detection threshold of 5 MeV. 
For the latest phase III of data taking, a flux 
of (2.32±0.04 stQt ±0.05 sast )xl0 6 cm- 2 s- 1 has been 
determined from neutrino-electron scattering |199j . 
However, the low-energy upturn in the 8 B v e spec- 
trum that is predicted by the MSW-LMA solution has 
not been observed. The new phase IV of data taking 
will feature a lower threshold of ^4 MeV and so will 
further explore the vacuum-matter transition region. 

SNO. Located in a mine in Ontario, SNO exploited 
the same Cherenkov technique of Super-K, with the 
difference that the detection medium was heavy water. 
The neutral and charge current neutrino reactions on 
deuterium provided the experiment with a powerful 
tool enabling to measure concurrently the total 8 B all 
flavor neutrino flux (neutral current) and the electron 
neutrino only flux (charge current). The SNO result 
verified unambiguously that the SNP was due to the 
conversion between different neutrino flavors, as im- 
plied by the MSW paradigm. 

The experiment provided also an accurate mea- 
sure of the total 8 B flux. SNO progressed through 
three steps (pure heavy water, salt and 3 He coun- 
ters) that returned consistent result. The most re- 
cent low energy threshold (LETA) joint analysis of 
the phase I and II data results in a 8 B flux of 
(5.140lg;l|g(sia*)lg;i??( a i/at))-10 8 cm- a B- 1 [200], in 
good agreement with both high- and low-metallicity 
predictions of the SSM. The detector is now empty, 
ready to be filled with liquid scintillator for the future 
SNO+ data taking phase. 

Borexino. While the two Cherenkov experiments fo- 
cused their investigations to the high energy portion 
of the 8 B spectrum, Borexino [201] at Gran Sasso low- 
ered for the first time the research range of a real time 
solar neutrino experiment down to few hundreds keV, 
based on the much larger light output of the scintilla- 
tion technique. This allowed to test the LMA vacuum 
oscillations below 1 MeV, in contrast to the matter- 
dominated regime probed by Super-Kamiokande and 
SNO. 

Achieving the ultra-low radioactive background con- 
ditions required for the detection of the 7 Be-i/s (0.862 
MeV) poses an enormous technological challenge. 
However, the necessary techniques for purification of 
the scintillator and the selection and assembly of low- 
background materials were developed in an extensive 
R&D program, culminating in the clear detection of 
the 7 Be-7v recoil electrons. The corresponding evalua- 
tion of the v e survival probability was in good agree- 
ment with MSW-LMA and SSM predictions, obtain- 
ing a value of (5.18±0.51)xl0 9 cm~ 2 s _1 for the total 
7 Be flux. 
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FIG. 19: Global solar and KamLAND data oscillation 
analysis |2U4j . 

As further proof of the powerful flexibility of the 
scintillation technique, Borexino performed a mea- 
surement of the 8 B-/y flux above 3 MeV, achieving 
the currently lowest threshold in a spectral mea- 
surement. The result corresponds to a flux of 
(2.4±0.4 stat ±0.1 sa;jt ) xl 6 cm- 2 s" 1 [202], in excellent 
agreement with the Super-Kamiokande measurements. 

KamLAND. The current understanding of the ex- 
perimental solar neutrino results in term of the neu- 
trino oscillation paradigm heavily relies also on the 
outcomes of th e Ka mLAND [203] reactor neutrino ex- 
periment (Sec. |5.2[ ). By comparing the theoretically 
expected antineutrino spectrum from a number of nu- 
clear power plants at different distances to the experi- 
mental site, with the measured spectrum, KamLAND 
was able to detect in the latter the clear imprinting 
of the oscillation effect, thus independently ruling out 
the other possible explanations for the solar neutrino 
deficit. 

The collective analysis of the data from all the so- 
lar neutrino experiments performed so far, plus those 
coming from KamLAND, puts stringent limi ts o n the 
Am 2 2 and #12 oscillations parameters. Fig.[l9j from 
reference |203j . shows the allowed region in the pa- 
rameters space, stemming from a two flavor oscilla- 
tion analysis. In a first approximation, the strongest 
constraint on Am 2 2 comes from KamLAND, while the 
limit on the mixing angle derives from the solar data. 

4.3.3 LENA observables and capabilities 

Despite the impressive successes accumulated in this 
field in the past, still additional and important insights 
can be expected from the detection of solar neutrinos. 
With a first measurement of pep and CNO neutrino 
fluxes, Borexino and the upcoming SNO+ experiment 
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Source Channel EW [MeV] m fid [kt] Rate [cpd] 



pp 


ve -} ev >0.25 


30 


40 


pep 


0.8-1.4 


30 


2.8xl0 2 


7 Be 


>0.25 


35 


l.OxlO 4 


8 B 


>2.8 


35 


79 


CNO 


0.8-1.4 


30 


1.9xl0 2 


8 B 


ia C >2.2 


35 


2.4 



TABLE IX: Expected solar neutrino rates in LENA. The 
estimates are derived from the existing Borexino analy- 
ses [2011 1202) as well as expectation values for the respec- 
tive energy windows (EW) for observation [9, 205, 206 . 
The quoted fiducial masses (mad) in LAB are based on a 
Monte Carlo simulation of the external 7-ray background 
in LENA [207] . 



will probe oscillations in the MSW transition region 
and solar metallicity, respectively. Even a direct mea- 
surement of the fundamental pp-i^ might be within 
reach of Borexino. 

However, the high-statistics data collected by a gi- 
gantic scintillation detector like LENA would allow 
a precise determination of SSM neutrino rates and 
MSW-LMA oscillation probabilities. The benchmark 
experience to be taken as reference is that of Borexino, 
to date the only liquid scintillator experiment which 
has successfully accomplished the detection of low- 
energy solar neutrinos. While the expected perfor- 
mances of a large 50 kt detector will not equal those 
of the smaller 0.3 ktons Borexino detector, especially 
because of the likely inferior photoelectron yield, the 
extremely high neutrino event rate resulting from the 
huge target volume will enable not only a detailed 
study of the features of the neutrino spectrum, but 
will allow also a thorough investigation of even small 
time modulations possibly embedded in the recorded 
flux. 

Again building on the Borexino experience, we may 
anticipate the need of a smaller fiducial volume, com- 
pared to other measurements, in order to cope suc- 
cessfully with the external gamma rays background, 
mainly from the photomultiplier tubes. MC simula- 
tions point to a fiducial mass of ~30 kt for pep, CNO 
and low-energy 8 B-^ detection, while the fiducial mass 
adopted for 7 Be-^s and high-energy (E > 5 MeV) 8 B- 
vs will be 35 kt or more. 



Table IX lists the expected rates in 30 kt for the 
neutrinos emitted in the pp chain and the CNO cycle, 
using the most recent solar model predictions. This 
evaluation refers to a detection threshold set at about 
250 keV, a lower threshold being severely hindered by 
the intrinsic abundance of 14 C in the scintillator com- 
ponents. 

Spectral measurements. Based on [S] and a fiducial 
mass of 30 kt, about 40 pp neutrino-induced electron 
backscattering events per day are expected above the 
threshold: despite the non-negligible rate, it is hard to 
anticipate the capability to distinguish them from the 



huge tail of the 14 C, especially taking into account the 
limited resolution that can be expected in this energy 
range. On the other hand, 7 Be detection will occur 
with enormous statistics, the prediction amounting to 
almost 10 4 7 Be recoils events per day in 35 kt. In the 
assumption that background levels similar to those of 
Borexino will be achieved, such a high statistics will 
permit a measurement of the 7 Be flux with accuracy 
unprecedented in neutrino physics. In particular, ac- 
curate search for temporal variations in the detected 
rate will be possible (see below). 

Roughly in the energy range between 1 and 2 MeV, 
detection of CNO and pep solar neutrinos can occur. 
A major background to this effort will be formed by 
cosmogenic n C beta decays, induced by the muon- 
induced knock-out of neutrons from 12 C. The 11 C pro- 
duction rate is mainly a function of the rock over- 
burden shielding the detector. If LENA will be oper- 
ated at the intended depth of 4000 mwe (meters wa- 
ter equivalent), the ratio of the CNO/pep-v signal to 
n C background rate would be 1:8, a factor 3 better 
than for example at the depth of Gran Sasso. About 
500 CNO-^s per day will provide valuable informa- 
tion on solar metallicity, especially if the contributions 
from the individual subfluxes can be disentangled. In 
Borexino, the /3-decays of 210 Bi originating from 210 Pb 
dissolved in the scintillator prove to be the most severe 
background after n C subtraction. As 210 Pb is rather 
long-lived, a special effort will be needed to keep the 
initial contamination of the scintillator as low as pos- 
sible. Furthermore, the measurement of the pep neu- 
trino flux could be exploited for a precision test of the 
v e survival probability in the MSW-LMA transition re- 
gion. The transition onset can be probed via the low 
energy portion of the 8 B neutrinos spectrum, through 
the accurate detection of the expected spectral "up- 
turn". The detection threshold for 8 B-i/s might be 
lowered even further than in Borexino, as the back- 
ground due to the penetrating 2.6 MeV 7-rays from 
external 208 T1 decays can be avoided by adjusting the 
fiducial volume. 

Finally, the charged-current reaction on 13 C should 
be mentioned. The channel is only accessible to v e s 
and virtually background-free due to the delayed coin- 
cidence of the 13 N back decay. The reaction threshold 
is 2.2 MeV, allowing for a precise measurement of the 
v e survival probability of 8 B-i/s in the MSW transition 
region. About 8xl0 2 counts per year are expected. 

Search for time- variations. The enormous amount 
of solar neutrino events collected in LENA will offer 
the possibility to search for temporal modulations in 
the neutrino flux arriving at Earth, especially regard- 
ing the 7 Be signal. Various processes that might cause 
such processes have been suggested: Apart from the 
annual modulation induced by the eccentricity of the 
terrestrial orbit, changes in the survival probability of 
solar v e s might be induced by fluctuations of the so- 
lar matter density and magnetic field [208] . or by the 
transit of terrestrial matter |209j before reaching the 
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detector. Even solar neutrino production rates might 
vary in the course of the solar cycle of about 11 years 
210^, or might be subject to short-term variations cor- 
related to the oscillation of the solar core temperature 
due to helioseismic waves |211j . 

Currently, the best limits on periodical v flux vari- 
ations arise from the Super-Kamiokande and SNO ex- 
periments, excluding modulation amplitudes of more 
than 10% in the 8 B-i/ signal [HQ 1212] • The investi- 
gated range of modulation periods extends from the 
order of hours to years. However, due to the ^10 4 
7 Be-z/ events per day available in LENA, the sensitiv- 
ity to modulations at low amplitudes is expected to 
be far greater: The MC analyses performed in |213j 
point to a 3cr discovery potential for amplitudes as 
low as 0.5%, covering a period range extending from 
tens of minutes to a hundred years or more. This will 
allow to probe the high-frequency regions associated 
to helioseismic g-modes, but also to test the tempo- 
ral uniformity of solar fusion processes on long time 
scales. 

4.4 Indirect dark matter search 

Dark matter particles might be abundantly present in 
the Universe and able to annihilate (decay) efficiently 
into Standard Model particles, in particular neutrinos, 
in regions where they are highly concentrated. We 
consider these annihilations (decays) in the galactic 
halo and show how LENA could be used to set general 
limits on the dark matter annihilation cross section 
and on the dark matter lifetime. 

4.4.1 Introduction 

With the next generation of neutrino experiments 
we will enter the era of precision measurements in neu- 
trino physics. These detectors, and specifically LENA, 
thanks to their great capabilities, might also be used 
to test some of the properties of the dark matter (DM). 

DM is copiously present in the Universe, having 
been produced in its very first instants. In the sim- 
plest case, the DM particles were in equilibrium in the 
Early Universe thermal plasma, decoupling when their 
interactions become too slow compared with the ex- 
pansion of the Universe. After decoupling, a thermal 
distribution remains as a relic which constitutes the 
DM we observe today. For the simplest assumption of 
thermal freeze-out, which holds in most of the models 
of DM, the annihilation cross section required to re- 
produce the observed amount of dark matter is given 
by (<tav) — 3x 10~ 26 cm 3 /s. Subsequently, within the 
framework of cold DM, structure forms hierarchically 
with DM collapsing first into small haloes and even- 
tually giving rise to larger ones, as galaxies and clus- 
ters of galaxies. Large concentrations of DM emerge, 
for example in the center of galaxies, such that the 
DM particles could annihilate efficiently and produce 
detectable fluxes of Standard Model (SM) particles, 



such as photons, neutrinos, positrons and antipro- 
tons. Such particles could be produced also if DM 
is not stable but decays with a lifetime longer than 
the age of the Universe in order to be present to- 
day. Among these particles, neutrinos are the least 
detectable ones. Therefore, if we assume that the only 
SM products from the DM annihilations (decays) are 
neutrinos, a limit on their flux, conservatively and in a 
model-independent way, sets an upper (lower) bound 
on the DM annihilation cross section (lifetime). This 
is the most conservative assumption from the detec- 
tion point of view, that is, the worst possible case (see 
|214j for a discussion on the implications of other de- 
cay modes for DM decays). Any other channel (into 
at least one SM particle) would produce photons and 
hence would give rise to a much more stringent limit. 
The bounds so obtained are on the total annihilation 
cross section (lifetime) of the DM particle and not only 
on its partial annihilation cross section (lifetime) due 
to the annihilation (decay) channel into neutrinos. 

In this section, and following and reviewing the ap- 
proach of |214H217j , we consider this case and evaluate 
the potential neutrino flux from DM annihilation (de- 
cay) in the whole Milky Way, which we compare with 
the relevant backgrounds for detection. In such a way, 
we obtain an estimate of the sensitivity on the DM 
annihilation cross section (lifetime) by LENA. 

4.4.2 Searching neutrinos from MeV DM 

For energies below ^200 MeV, information on the 
direction of the incoming neutrino is very poor if the 
detection is via interactions with nucleons, as is the 
case in LENA. Thus, we take the flux averaged over 
the entire galaxy. The differential neutrino or antineu- 
trino flux per flavor from DM annihilation (decay) av- 
eraged over the whole Milky Way is given by 

^- = P k (E v , m x ) R Q p k Q J avg , k , (4.2) 

where m x is the DM mass, po = 0.3GeVcm~ 3 is a 
normalizing DM density at R Q = 8.5 kpc (the distance 
from the Sun to the galactic center), and J aV g,k is 
the average over the whole galaxy of the line of sight 
integration of the DM density (for decays, k = 1) or 
of its square (for annihilations, k = 2), which is given 

by 

Javg,k = p 1 fc / / p(r) k dld(cosip'), 

Z K G Po J -I JO 

r = - 2lR Q cosip' +l 2 , 

*max = s] (R 2 vir - sin 2 il>R%) + R G cos V, (4.3) 

i? v ir being the halo virial radius. Commonly used pro- 
files [2181 220 tend to agree at large scales, although 
they may differ significantly in the inner part of galax- 
ies. The overall normalization of the flux is affected 
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by the value of J aV g,k, scaling as p k . For instance, for 
different profiles 218 -220 , astrophysical uncertainties 
can induce up to a factor of ~100 (~6) for annihila- 
tions [2161 I217j (decays 1214) ). For concreteness, in 
what follows we present results using the Navarro, 
Frenk and White (NFW) profile [219], with J avg i = 2 
HH] and J avg>3 = 5 [HI]. 

All the dependence on the particle physics model is 
embedded in Vk as 



V-2 



1 dN x 1 

3 dE v rn x T x 
1 dN? (<j a v) 
3 dE v 



2 m 2 



for decays and 



for annihilations, 



where the neutrino or antineutrino spectrum per flavor 
is given by 



dNx 
~dE v 
dN 2 
dE v 



5(E V - m x ) 



for decays and 
for annihilations. 



The factor of 1/3 arises from the assumption that the 
annihilation or decay branching ratio is the same for 
the three neutrino flavors. This is not a very restrictive 
assumption, for even when only one flavor is produced, 
a flux of neutrinos in all flavors is generated by the 
averaged neutrino oscillations between the source and 
the detector. 

The neutrinos produced in DM annihilations (de- 
cays) travel to the Earth where they can be revealed 
in present and future neutrino detectors. Importantly, 
the signal is monoenergetic, allowing to distinguish it 
from backgrounds continuous in energy. The number 
of signal neutrino events is given by 



Af ~ cr dot (f) N t 



argct 



te, 



(4.4) 



where the detection cross section <7dct needs to be eval- 



uated at E v 



,/2) for annihilations (de- 



cays), the total flux of neutrinos or antineutrinos is 
given by </>, AT target indicates the number of target par- 
ticles in the detector, t is the total time-exposure, and 
e is the detector efficiency for this type of signal. From 
Eq. (4.4), and assuming the annihilation cross section 



required to reproduce the observed amount of DM, 
(a A v) = 3 x l(r 26 cm 3 /s, (or a lifetime r x - 10 24 s) 
we expect a few events for an exposure of 1 Mt-yr, re- 
quiring large detectors such as LENA. 



4.4.3 MeV Dark Matter search in LENA 



It has been shown |214l I217j , particularly for Super- 
Kamiokande, that already present large neutrino de- 
tectors severely constrain the DM properties, namely 
annihilation cross section and lifetime and, depend- 
ing on the DM profile assumed, exclude an important 
part of the parameter space (see Fig. 1 in |217) . Fig. 1 
in [Hi] and also Fig. 1 in [22T] ). 



ID 

9 
8 

\ 7 
S 

? 6 

>» 

o 

2 5 

I 3 

Z 

2 

1 





DM Signal + Background 
Background 



50 60 
E . [MeV] 



FIG. 20: Expected signal in the LENA detector, in 
Pyhasalmi, after 10 years of running for two values of the 
DM mass, m x = 20 (40) MeV and m x = 60 (120) MeV 
for annihilations (decays) for {uav} = 3xl0~ 26 cm 3 /s 
(r x = 8.9xl0 23 s and r x = 2.7 x 10 24 s, respectively for 
each mass). Dashed lines represent the individual contri- 
butions of each of the three types of background events in 
this type of detector (reactor v e , DSNB and atmospheric 
neutrinos). The solid lines represent the backgrounds plus 
the expected signal from DM annihilation (decay) in the 
Milky Way. Taken from [217] . 



Here we describe the analysis performed in |217j for 
the physics reach of the LENA detector. The excel- 
lent background rejection allows for a significant im- 
provement on present bounds. At the energies of in- 
terest, few tens of MeV, the inverse beta-decay cross 
section (i> e p — > ne + ) is by two orders of magnitude 
larger that the v — e elastic scattering cross section. 
The advantage of the LENA detector is the excellent 
energy resolution and the fact that the inverse beta- 
decay reaction can be clearly tagged by the signal in 
coincidence of the positron annihilation followed by a 
delayed 2.2 MeV photon, which is emitted when the 
neutron is captured by a free proton. Thus, the only 
relevant backgrounds for these events come from reac- 
tor, atmospheric and diffuse supernova v e interacting 
with free protons in the detector: The flux of reac- 
tor i> e 's below ~ 10 MeV represents a background by 
orders of magnitude higher than the expected neu- 
trino flux from DM annihilations (decays) [3J. The 
Diffuse Supernova Neutrino Background (DSNB), al- 
though not yet detected, might potentially represent a 
background in the interval ~10-30MeV. Atmospheric 
zVs constitute the dominant background in the en- 
ergy range above ~30 MeV. The normalization of the 
flux depends on the location of the detector [222H225] , 
more specifically on the geomagnetic latitude, varying 
roughly within a factor of 2 [3]. Moreover, NC reac- 
tions on carbon will create a background that is not 



yet well determined. See Sect. |4.2| for a closer discus- 
sion of these backgrounds. 



In Fig. 20 (taken from |217j ). the expected sig- 
nal and background spectra are shown for LENA in 
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Pyhasalmi after 10 years of data taking. Other lo- 
cations return similar results. The assumed scintilla- 
tor mixture is 20% PXE (G 16 H 18 ) and 80% Dodecane 
(C12H26) for a fiducial volume of 50xl0 3 m 3 , which 
amounts to 3.3 x 10 33 free protons. The rates and spec- 
tra are calculated using a gaussian energy resolution 
function of width [3 

o-lena = 0.10 MeV y/E/MeV. (4.5) 

Two values for the DM mass are depicted: m x = 
20 (40) MeV and m x = 60 (120) MeV for annihilations 
(decays) for (a A v) = 3xl0- 26 cm 3 /s (r x = 8.9xl0 23 s 
and t x — 2.7xl0 24 s, respectively for each mass). In 
the case of low values of the DM mass, even with the 
small rate predicted, a rather easy discrimination be- 
tween signal and background could be possible. For 
higher values of the masses, the energy of the initial 
neutrino cannot be precisely reconstructed from the 
measured positron energy. Therefore, the signal at 
these masses has a spread over an interval of ~10 MeV, 
degrading the sensitivity of the detector. 

4.4.4 Conclusions 

Determining the DM identity and its properties is 
one of the fundamental questions to be answered in 
the future in astroparticle physics. In regions of the 
Universe where DM is highly concentrated, such as the 
center of galaxies, DM particles can annihilate (decay) 



efficiently producing observable fluxes of SM particles. 
Large neutrino detectors might be able to observe the 
neutrino so produced providing bounds or measure- 
ments of the DM mass and annihilation cross section 
(lifetime). It should be noted that neutrinos are the 
least detectable particles of the SM and therefore pro- 
vide the most conservative bounds on DM annihila- 
tions (decays). Importantly the DM neutrino signal is 
mono-energetic, allowing for an enhanced discrimina- 
tion between signal and continuous backgrounds. 

LENA would be particularly suited to these searches 
thanks to the large size, the excellent energy resolu- 
tion and the good background discrimination. For DM 
masses in the few tens of MeV, it could observe a sig- 
nal if these MeV particles exist and the annihilation 
cross section is the one required to reproduce the ob- 
served amount of DM (or its lifetime r x ~ 10 24 s). In 
particular, the LENA detector would have the capa- 
bility to find a positive signal at ^2cr in a large part 
of the mass window of interest. A null signal in LENA 
would indicate that, if DM particles with mass ~10- 
100 MeV exist, then they must live longer than ~10 24 s 
and they were not produced thermally or the annihila- 
tion cross section at freeze-out was velocity-dependent. 
A positive signal would imply that DM is constituted 
by particles with masses in the tens of MeVs, would 
measure its mass and would determine the cross sec- 
tion which was relevant at DM freeze-out in the Early 
Universe (or its lifetime), for a given halo profile. 
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5 Terrestrial neutrino sources 



Not only astrophysical bodies, but also our own planet 
is a source of low-energetic neutrinos (Sec. 5.1 ). More- 
over, there is a number of anthropogenic neutrino 
sources intense enough at low energies to allow for neu- 
trino oscillation physics: Antineutrinos from nuclear 
reactors (Sec. |5.2[ ) may allow for a high-precision mea- 
surement of the "solar" neutrino mixing parameters. 
Neutrinos from strong radioactive electron capture 
sources may provide a unique opportunity to investi- 
gate flavor oscillations on a very short baseline, provid- 
ing sensitivity to the mixing angle 13 and especially 



v & disappearance into sterile neutrinos (Sec. 5.3). Al- 



ternatively, low-energy neutrinos generated by a pion 
decay at-rest beam might offer sensitivity to #13 and 
especially the CP-violating phase 5cp (Sec. 5.4 1. Dot- 
ing the scintillator target might allow to search for 
neutrinoless double-beta decay (Sec. 5.5). Finally, a 
large low-energy detector is sensitive to many other 
rare processes connected to non-standard interaction 



of neutrinos or exotic decay processes (Sec. 5.6) 



5.1 Geoneutrinos 

Geoneutrinos (geo-i/s) are D e s produced inside the 
Earth during /3-decays of naturally occurring radioac- 
tive elements. They are direct messengers of the abun- 
dances and distribution of radioactive elements within 
our planet, information strongly constraining all geo- 
chemical and geophysical models. Geoneutrinos have 
been successfully detected by the liquid-scintillator ex- 
periments KamLAND and Borexino. However, the ge- 
ological information contained in these measurements 
is still limited, mostly because of low statistics. A 
multi-kT detector like LENA, featuring the radiopu- 
rity already achieved by Borexino, would address sev- 
eral questions of geological importance. This section 
presents the event and background rates expected for 
LENA (both in Pyhasalmi and Frejus), and projects 
the precision at which the total geo-z/ flux as well as 
the U/Th ratio could be measured. 

5.1.1 Introduction 

Geo-^s originate from the /3-decays of radioactive 
elements in the Earth's crust and mantle, predomi- 
nantly from 40 K and nuclides in the chains of 238 U 
and 232 Th. These neutrinos probe direct informa- 
tion about the absolute abundances and distribution 
of these radioactive elements inside the Earth. Their 
measurement quantifies the radiogenic contribution to 
the total heat flux of the Earth, constraining geochem- 
ical and geophysical models of the planet. This infor- 
mation provides constraints on the many and complex 
processes that operate inside the Earth, including the 
generation of the Earth's magnetic field, mantle con- 
vection, and plate tectonics. In addition, determining 
the absolute abundances of refractory elements (i.e. U 



and Th) in the planet provides insight into its origin 
and formation. 

Estimates of the Earth's surface heat flux 
emerge from temperature gradient measurements from 
^40 000 drill holes distributed around the globe. Us- 
ing these data, geophysical models typically conclude 
that the present surface heat flux is 47±2TW [226] . 
This conventional view has been challenged by an al- 
ternative flux estimate of 31±1TW [227] . A signifi- 
cant contributor to this heat flux comes from the heat 
producing elements, K, Th and U, with its flux pro- 
portion dependent upon their absolute abundance in- 
side the Earth. The many models that describe the 
composition of the Earth come from cosmochemical, 
geochemical and geophysical observations and predict 
a range of abundances and distributions of these ele- 
ments [228H23T] . 

The Earth has a silicate shell, the Bulk Silicate 
Earth (BSE), surrounding a metallic core, with the 
core being an iron-nickel mixture (with proportions 
set from cosmochemical constraints) often considered 
to contain negligible quantities of Th and U |232[ 1233] . 
The BSE describes the primordial, non-metallic Earth 
condition that followed planetary accretion and core 
separation, prior to its differentiation into a man- 
tle, oceanic crust, and continental crust. Elements 
excluded from the Earth's core are referred to as 
lithophile and those that accreted onto the Earth in 
chondritic proportions are the refractory elements. 

Chondritic meteorites are undifferentiated samples 
with refractory element abundances in equal propor- 
tion and record a high temperature condensation char- 
acteristic of the cooling nebular. An important guide 
to predicting planetary compositions is given by the 
compositional match between chondritic meteorites 
and the solar photosphere on a one to one basis, over 
5 orders of magnitude for the non-gases, based on an 
equal atom abundance of silicon. 

Thorium and uranium are refractory lithophile ele- 
ments and contribute equally ^80 % of the total radio- 
genic heat production of the Earth, while the remain- 
ing fraction is due to 40 K, a volatile element (assuming 
Th/U~4 and K/U~10 000) . During mantle melting 
and because of their chemistry and size, K, Th and 
U are quantitatively partitioned into the melt and de- 
pleted from the mantle. Thus, the continental crust, 
has over geologic time, been enriched in these elements 
and has a sizable fraction (about half) of the planet's 
inventory, producing radiogenic power of 7.3±2.3 TW 
(2a) [23"t| . 

The range of BSE models predicting the Th, U, and 
K abundances (Tab.fxj) translates to radiogenic heat 
contributions of 12-30 TW, and thus allow other possi- 
ble heat sources to make up the total surface heat flux. 
Additional heat might might originate from accretion, 
gravitational contraction, latent heat from phase tran- 
sitions, or from a (rather exotic) nuclear reactor in the 



42 



core/core-mantle boundary (CMB). Systematic errors 
in both geochemical and geophysical models are not 
very well known and the validity of several assump- 
tions on which they are based is not proven. Thus, ob- 
servations of the planetary geo-v flux will yield trans- 
formational insights into the Earth's energy budget. 

Typically, based on geophysical calculations, pa- 
rameterized convection models of the mantle require 
higher radiogenic heat contributions (~70 % of the to- 
tal heat flux) in order to describe the Earth's cooling 
history in terms of a balance of forces between ther- 
mal dissipation and mantle viscosity. In contrast, geo- 
chemical models using cosmochemical and geochemi- 
cal observations predict the BSE abundance of U and 
values for Th/U and K/U, 4 and 1.4x 10 4 , respectively, 
with an uncertainty of ~10%. Consequently, a geo- 
chemist's view of the Earth predicts that its budget of 
heat producing elements in the BSE are up to a factor 
of ^3 lower than the models predicted by geophysi- 
cists. Thus, the relative contribution of the radioac- 
tive power to the total planetary heat flux is poorly 
known. 

The first ideas regarding geo-z^s are from the sixties 



238 but only recently the first experimental results 
from large volume scintillator detectors are available. 
KamLAND and Borexino have observed geo-neutrinos 
with similar precision [2031 1239H241 ] . Borexino has 
significantly greater radiopurity and lesser reactor flux 
than KamLAND, allowing Borexino to make the mea- 
surement in 5% of the exposure of KamLAND. In ad- 
dition, both experiments have placed limits (<3TW) 
on the potential contribution of a putative geo-reactor 
deep in the Earth's interior. 

In spite of these first successful experimental re- 
sults, it has not been possible yet to neither discrim- 
inate among several predictions concerning the radio- 
genic heat production nor measure Th/U ratio of the 
observed energy spectra. Even the very low back- 
ground measurement performed by Borexino has lim- 
ited power of geological predictions due to slowly ac- 
cruing statistics. Several future experiments, as for 
example SNO+ project in Canada, have among their 
aims geo-v measurements. However, a real break- 
through in this field would come with a very large 
volume detector at 50 kt scale, like LENA. 

In liquid scintillator detectors, the D e are detected 
via the inverse beta decay, with a kinematic threshold 
of 1.8 MeV (Sect. 4.2 1. The characteristic time and 
spatial coincidence of prompt e + and delayed neutron 
events offers a clean signature. Since 1 kt of liquid scin- 
tillator contains about 10 32 free protons (the precise 
value depending on the chemical composition) and the 
exposure times are of order of a few years, the events 
rates are conveniently expressed in terms of a Terres- 
trial Neutrino Unit (TNU), defined as one event per 
10 32 target protons per year. 

The geo-z^ flux produced from U and Th inside the 
Earth is some 10 6 /cm 2 s. Due to neutrino oscillations, 
the anti-neutrino flux arriving at the detector in the 



electron flavor will be smaller than that produced: for 
our calculation we have considered an asymptotic sur- 
vival probability (Pee) — 0.57 following the best fit 
obtained in |242j . Only a small fraction (about 5%) 
of D e from the 238 U and 232 Th series are above the in- 
verse beta-decay reaction threshold, while those from 
40 K decays are below this threshold. Geo-z^s originat- 
ing from different elements can be distinguished - at 
least in principle - due to their different energy spectra, 
as only vs from the uranium chain contribute at ener- 
gies E e + > 2.25 MeV. The exact spectrum depends on 
the shapes and rates of the individual decays within 
U and Th chains, and on the abundances and spa- 
tial distribution of U and Th in the crust and in the 
mantle. 

The geo-v signal spectrum extends to E e + ss 
2.6 MeV. However, v e from nuclear power plants repre- 
sent a background for geo-z^ detection, E e + extending 
up to ~10MeV. In the following the expected geo-z' 
and reactor v signals at Pyhasalmi and Frejus sites 
are discussed. The potential of the LENA project to 
achieve geologically interesting results is discussed as 
well. 

5.1.2 The geoneutrino signal 

Different calculations for geo-z' production have 
been presented in the literature [243-246]): all models 
rely on the geophysical 2° x 2° crustal map of |247j and 
on the density profile of the mantle as given by the Pre- 
liminary Reference Earth Model (PREM) [248]. For 
the calculation of geoneutrino signal we adopt the val- 
ues of U and Th abundance recommended in |249| for 
the sedimentary layers and the values reported in [234] 
for the upper, middle and lower crust (Tab. XI). The 



la uncertainties for the upper and middle crust are 
from [234] . For the lower crust, we adopt an uncer- 
tainty indicative of the spread of published values. 

The composition and the circulation inside the 
Earth's mantle is the subject of a strong and so far 
unresolved debate between geochemists and geophysi- 
cists: assuming that a spherical symmetry holds and 
U and Th abundances do not decrease with depth, 
the extreme predictions for the signal are obtained 
by placing U and Th in a thin layer at the bottom 
and distributing it with uniform abundance over the 
mantle. For a fixed total U mass in the BSE model, 
m(U) = 0.8xl0 17 kg, and a fixed ratio of the elemen- 
tal abundance Th/U = 3.9 |232j . the contribution to 
the geo-z^ signal from the crust and the mantle is ob- 
tained by using the proximity argument presented in 
250 : the minimal (maximal) contributed flux is ob- 
tained by placing uranium and thorium as far (close) 
as possible to the detector. 

Our prediction for geoneutrino signal is obtained 
by the mean of these extremes, assigning an error 
so as to encompass both of them. For the cen- 
tral value of BSE model, the expected signal from U 
and Th at Pyhasalmi is 51.3±7.1TNU. At Frejus, it 
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Authors 






«(U) [ng/g] 


o(Th)/a(U) 


ff M (U + Th) [TW] 


Turcotte and Schubert 


(2002) 


[235] 


31 


4.0 


19 


Anderson 


(2007) 


230 


28 


4.0 


17 


Palme and O'Neill 


(2003) 


236 


22 


3.8 


12 


Allegre et al. 


(1995) 


237 


20 


3.9 


11 


McDonough and Sun 


(1995) 


228 


20 


3.9 


11 


Lyubetskaya and Korena^ 


;a (2007) 


229 


17 


3.7 


7 


Javoy et al. 


(2010) 


[2"3T] 


12 


3.5 


3 



TABLE X: Uranium content a(U) in the Bulk Silicate Earth, the Th/U ratio and the radiogenic heat production in the 
mantle (Hm) due to U and Th according to different authors. 



Reservoir 



Units 



o(U) 



a(Th) 



Sediments pg/g 1.68 ± 0.18 6.91 ± 0.8 
Upper Crust pg/g 2.7 ± 0.6 10.5 ± 1.0 
Middle Crust pg/g 1.3 ± 0.4 6.5 ± 0.5 
Lower Crust pg/g 0.6 ± 0.4 3.7 ± 2.4 

Oceanic Crust pg/g 0.1 ± 0.03 0.22 ± 0.07 

TABLE XI: U and Th mass abundances in the Earth's 
reservoirs [23411249] . 



Location 


Signal 1-10 MeV Sij 


mal 1-2.6 MeV 




[TNU] 


[TNU] 


Pyhasalmi 


70.9±3.8 


20.8±1.1 


Pyhasalmi* 


145.9±7.7 


37.3±1.9 


Frejus 


554±29.4 


145±7.7 



TABLE XII: Expected reactor v e signal, ('future Finnish 
reactors are taken into account). The corresponding energy 
spectra are shown in Fig.|22| 



is 41.4±5.6TNU: the accuracy of about 14% corre- 
sponds to "3cr". In a target mass of 44 kt, correspond- 
ing to 2.9xl0 33 free protons, we expect a eeo-v signal 

the ex- 
geo-vs at 
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of the order of 10 3 events/year. In Fig. 
pected signal S(U + Th) from U and T 
Pyhasalmi and Frejus is shown as a function of the 
radiogenic heat production rate H (U + Th). 

For a given total uranium mass in the Earth, to(U), 
corresponding to a fixed radiogenic heat production 
i/(U + Th), the minimal and maximal signals are pro- 
vided by the terrestrial models consistent with avail- 
able geochemical and geophysical observational data 
and by proximity argument |250j . The uncertainty 
band is wide because the signal is dominated by con- 
tribution from the crust: a refinement of the reference 
model taking into account the regional contribution is 
appropriate. Considering that some 50 % of the signal 
from the crust originates from a region within 200 km 
from both detectors, a better geological and geochemi- 
cal description of the regions surrounding the detectors 
is needed for a more precise estimate of the geoneu- 
trino signal. A comprehensive study of the crust near 
the detector is essential for resolving global models. 

5.1.3 Reactor neutrino background 

The expected reactor v e flux was calculated based 
on the same assumptions as described in |240j . The 
neutrino oscillations parameters (Am 2 2 = 7.65 -10 -5 
eV 2 ; sin 2 6 12 = 0.304) from Ref. [251] were used. The 
monthly load factor of year 2009 was considered for all 
493 world nuclear reactors [252 . The expected reactor 
v e signal at Frejus and Pyhasalmi and shape of the 
oscillated spectrum are given in Tab. |XIl| and Fig.[22] 
respectively. 



For Pyhasalmi site also the case of possible future 
reactors was considered, assuming a typical 80 % load 
factor. In particular, the Olkiluoto-3 power plant is 
under construction and should be operating in 2013 
with 4.3 GW thermal power. At the same site, 360 km 
from Pyhasalmi, the construction of Olkiluoto-4 reac- 
tor with power up to 1.8 GW was approved. In addi- 
tion, the construction of additional reactor with power 
up to 4.9 GW is under approval (Pyhajoki site 130 km 
from Pyhasalmi is among possible sites). 

5.1.4 Determining the geoneutrino flux 

LENA, thanks to its large volume, would be a real 
breakthrough in the field of geo-^ detection. Within 
the first year, geologically significant results could 
be obtained. Independently from the final location, 
within the first year LENA would measure the total 
geo-f flux at the level of few percent, by far more 
precise than the current experiments as Borexino or 
KamLAND could reach. 

We simulated the expected energy spectrum at 
Pyhasalmi (with and without future Finnish reactors) 
and Frejus, based on 2.9xl0 33 target protons and the 
expected reactor i/ e and geo-z^ fluxes described above. 
A light yield of 400 photoelectrons/MeV at the upper 
limit of the achievable range was assumed. The ex- 
pected spectra were convoluted with the consequent 
energy dependent resolution. 

LENA aims to reach, and possible exceed the radio- 
purity of Borexino detector. Therefore, the no- 
background approximation is reasonable, since Borex- 
ino has shown that the final v e spectrum is almost 
background free (the total background is less than 2 % 
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FIG. 21: The expected geo-v signal at Pyhasalmi (left) and Frejus (right) as function of radiogenic heat due to U and 
Th in the Earth H(U + Th). The area between the red and blue lines denotes the region allowed by geochemical and 
geophysical constraints. The green region is allowed by the BSE model according to [232] . 
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FIG. 22: The expected reactor v e signal at Pyhasalmi (left) and Frejus (right). The red spectrum (left) is when future 
reactors Finnish reactors are taken into account. Such spectra would be measured by an ideal detector with 10 32 free 
target protons in one year. 



of the total u e spectrum |240j ). 

The shape of the expected spectra is shown in 



Fig. 23 The chondritic U/Th ratio was assumed. The 
final precision of the geo-v flux measurement which 
could be reached after 1, 3, and 10 years is given in 
Tab. |XIli| Systematic uncertainties on the expected 
reactor background flux are not considered here and 
might limit the sensitivity if they surpass the statisti- 
cal error of the measurement. 



Live time Pyhasalmi Pyhasalmi* Frejus 



lyr 


3% 


4% 


6% 


3yrs 


2% 


2% 


3% 


lOyrs 


1% 


1% 


2% 



TABLE XIII: Expected precision in the measurement of 
the total geo-^ flux (*future Finnish reactors are taken into 
account). Details in text. 



5.1.5 Potential to measure the U/Th ratio 

The precision at which the U and Th fluxes, as well 
as their ratio, could be measured at Pyhasalmi and 



Frejus sites is summarized in Tab. XIV In these cal- 
culations, the assumptions taken into account are the 
same as in the previous section. However, instead of 
the single geo-neutrino component (with the U and 



Th ratio fixed to the chondritic value) the two indi- 
vidual contributions of U and Th chains are left free 
in the fit. 



The Pyhasalmi site is strongly preferred for this 
measurement. An example of a possible 5-yr energy 
spectrum of reactor and U and Th geo-^s is shown 
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FIG. 23: The expected oscillated visible energy spectrum at Pyhasalmi (left, future reactors in Finland are considered) 
and Frejus (right) for 1 year statistics and 2.9 xlO 33 target protons, considering a light yield of 400 pe/MeV. The 
contribution of the reactor v^s is shown by the filled orange area while that of geo-z/s by the dashed blue line. The yellow 
area isolates the contribution of the geo-i^s in the total signal. 



Location Live time U flux Th flux U/Th 





[yrs] 


[%} 


[%] 


[%} 


Pyhasalmi 


1 


6 


12 


17 




3 


3 


8 


10 




10 


2 


4 


5 


Pyhasalmi* 


1 


7 


14 


21 




3 


4 


8 


11 




10 


2 


4 


G 


Frejus 


1 


14 


25 


35 




3 


9 


12 


20 




10 


4 


7 


11 



TABLE XIV: Expected precision in the measurement of 
the U and Th geo-u flux and in the U/Th ratio. (* future 
Finnish reactors are taken into account). Details in text. 



in Fig. 24 together with the U vs Th contour plot 
resulting from the fit. 



5.1.7 Backgrounds 

LENA aims to achieve a similar or even better ra- 
diopurity level than realized in Borexino. Here, we 
briefly describe the background sources relevant for u e 
detection. 

9 Li- 8 He: 9 Li (T 1/2 = 178 ms) and 8 He (T 1/2 = 
119 ms) are /3 _ -neutron emitters, that are produced 
by cosmic muons crossing the detector. Borexino 
measured 15.4 events/(100 t-yr) |240j . which scales to 
about 1500 events per year in LENA, considering a 
reduced muon flux by a factor of 5. This background 
can be reduced to about 1 event per year in LENA, 
if a 2 s cut after every detected muon is applied. As 
9 Li and 8 He are produced close to the muon track, 
it is possible to reduce the dead time from ^6 % to 
~0.1 %, if only a cylinder with 2 m radius around the 
muon track is vetoed. 



5.1.6 Directionality 

In the inverse beta decay reaction the neutron is 
scattered roughly in forward direction. Thus, it is 
possible to obtain directional information about the v e 
event by measuring the displacement between the neu- 
tron and the e + event |253j . As the average displace- 
ment is with 1.9±0.4cm |254] rather small compared 
to the neutron and e + position reconstruction uncer- 
tainty, the v e direction can only be reconstructed with 
a large uncertainty. But, by analyzing a large number 
of geo-i' events, it is still possible to extract informa- 
tion about the v e angular distribution |254j . From 
the angular distribution of the geo-f events, the dif- 
ferential radial distribution of terrestrial radio-nuclides 
could be determined |255j . This would be important 
to differentiate between the geo-i'S coming from the 
mantle and those from the crust. 



Fast Neutrons: Cosmic muons that pass the detec- 
tor can produce fast neutrons. These neutrons have a 
large range and can reach the Inner Vessel of LENA 
without triggering the muon veto. In the scintillator 
they can mimic P e events, as they give a prompt sig- 
nal due to scattering off protons and a delayed signal 
caused by the neutron capture on a free proton. The 
fast neutron background rate was analyzed with a MC 
simulation [50]. In the geo-^ energy region, less than 
10 events per year are expected. Compared to the 
expected signal this background is negligible. 

13 C(a,n) 16 0: Neutrons can also be produced in 
the scintillator by 210 Po a decays and subsequent 
13 C(a,n) 16 reactions. If the radio-purity level of 
Borexino is reached in LENA, about 10±1 events per 
year are expected in LENA |240j . 




FIG. 24: Left: Expected 5-year visible energy spectrum for 2.9 xlO 33 target protons and light yield of 400 pe/MeV at 
Pyhasalmi site (future Finnish reactors not considered). The contribution of the reactor £- e s is shown by the filled orange 
area while the yellow area isolates the contribution of the geo-^s in the total signal. The data were generated with the 
chondritic U/Th ratio while in the fit the contribution from U (blue line) and Th (red line) were fit individually. Right: 
Corresponding contour plots (1,2,3,4, and 5 a C.L. outwards from the best fit black point, respectively) for the absolute 
U and Th number of events resulting from the fit. Each of the three solid lines corresponds to a fixed U/Th ratio. The 
central line corresponds to chondritic U/Th ratio, while the two external lines correspond to this ratio changed by ±20 %. 



5.2 Reactor neutrinos 

Experiments with reactor anti-neutrinos have a long 
and successful tradition in neutrino physics. Kam- 
LAND was the first reactor neutrino experiment to ob- 
serve a deficit in the flux, confirming the Large Mixing 
Angle MSW solution to the solar neutrino problem. 
Consequently, the experiment performed the most pre- 
cise measurement of the oscillation parameter Am^. 
A series of experiments at distances of several meters 
to ~lkm to the reactor core(s) has led to the current 
upper bound on the mixing angle #13, dominated by 
the result of the CHOOZ experiment 256 . A new 
generation of experiments with a multi-detector setup 
composed of detectors near and far (~lkm) to the 
cores aim to measure or constrain further the value of 
#13. The use of liquid scintillator detectors for reactor 
neutrino detection is a perfectly established technique 
and regardless which location is finally chosen for the 
realization of LENA, there will be a measurable reac- 
tor neutrino signal. 

While anti-neutrinos from nuclear reactors form a 
considerable background for the detection of geoneu- 
trinos and the DSNB in LENA, their signal also offers 
the opportunity to perform a high-statistics study of 
neutrino oscillation effects, especially to improve the 
knowledge on the parameters that drive solar neu- 
trino oscillations. Each reactor provides a high in- 
tensity, isotropic source of anti-neutrinos with a well- 
known initial spectrum, resulting from f}~ decays of 
fission products ( 235 U, 238 U, 239 Pu, 241 Pu being the 
four main isotopes) . The overall emitted anti-neutrino 



spectrum is computed 7 from measurements of beta 
spectra at ILL for 235 U, 239 Pu, and 241 Pu and the- 
oretical calculations for 238 U (see e.g. |258j for a poly- 
nomial parametrization) . During reactor operation, 
the abundance of 235 U decreases, while that of 239 Pu 
and 241 Pu increases. If the fuel evolution of all reac- 
tors is known, this burn-up effect can be taken into 
account. For the sensitivity studies presented here, 
a typical averaged isotopic composition is used. The 
contribution of stored spent fuel elements to the de- 
tected signal is considered negligible, as they contain 
mainly long lived emitters with a low Q value. 

In the LENA detector, anti-neutrinos are detected 
through the inverse /3-decay process on free protons 
v e + P — > e + + n (energy threshold of 1.8 MeV) with 
well known cross section, followed by neutron capture. 
Experimentally, the clear signature of the coincidence 
signal, formed by the prompt positron signal followed 
by the delayed neutron capture in spatial correlation 
can be used for powerful reduction of accidental back- 
ground. The visible energy E v i s of the prompt event 
is related to the energy of the incident neutrino E„ e 
by E vis = E 9e - m n +m p + m e . 

In |259j , the possibility of a high precision measure- 
ment of the solar mixing parameters Am 2 ^ and sin# 12 



New calculations of reactor anti-neutrino spectra from these 
measurements, including information from nuclear databases, 
have been performed in 12571 . While the shapes of the spectra 
and their uncertainties are comparable to that of the previous 
analysis, the absolute flux normalization is shifted by about 
+3% on average. 
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has been investigated, assuming a LENA type liquid 
scintillator detector that is located at Frejus. Here, the 
reactor neutrino flux is highest compared to other sites 
under consideration. Fig. [22] shows on the right the 
spectrum from inverse beta decays expected at Frejus. 
A dominant part of the total flux (67 %) is provided 
by the four nearest reactors within a distance of up 
to 160 km in Switzerland and France. As the authors 
of |259| point out, their distances are located between 
the first and the second survival probability minimum, 
and hence spectral information should provide a pow- 
erful tool to measure the oscillation parameters. 

The large flux originating from French and Swiss nu- 
clear power plants corresponds to a rate of of 1.7xl0 4 
inverse beta decay events per year in a fiducial mass of 
44 kt, two orders of magnitude larger than the Kam- 
LAND event rate. A threshold of 2.6 MeV on the visi- 
ble prompt energy was applied to eliminate the signal 
from geoneutrinos for these numbers. 

In this scenario, 3<t uncertainties below 3 % on 
Am|i and of about 20% on sin #12 could be obtained 
based on 1 year of exposure. After 7 years of data 
taking, the 3<r uncertainties would diminish to 1% in 
Aniji and 10% in sin0i 2 , respectively. While this 
would mean only a moderate improvement compared 
to present-day accuracies in the case of sin<9 12 , the 
uncertainty in determining the value of Am^ would 
decrease by almost an order of magnitude 260 . 

In case of the Pyhasalmi site, the total event rate 
coming from currently operating reactors is lower by 
almost on order of magnitude. A third reactor core at 
the Olkiluoto plant is under construction since 2005, 
and the permission for a fourth was approved in 2010. 
Approval for an additional reactor in a site at 130 km 
distance to Pyhasalmi is under discussion. The fu- 
ture situation with new reactors would correspond 
to a roughly doubled expected reactor neutrino flux. 
Therefore, the Pyhasalmi site will still be the preferred 
for the detection of geoneutrinos and the DSNB due 
to the lower reactor neutrino background. Neverthe- 
less, one can expect a useful total number of reactor 
neutrino events accumulated over the operation time 
of LENA, improving the determination of neutrino os- 
cillation parameters. 



5.3 Neutrino oscillometry 

An extended liquid-scintillator detector LENA offers 
the opportunity for neutrino oscillometry [261] . Based 
on a monoenergetic v e source, the characteristic spa- 
tial pattern of v e disappearance can be detected within 
the length of detector. Radioactive nuclides under- 
going electron capture produce monoenergetic neutri- 
nos: Sufficiently strong sources of more than 1 MCi 
activity are produced at nuclear reactors. In the three- 
flavor scenario, the investigation of the mixing param- 
eters 0i3 and Am^ 3 are the most promising due to the 
short oscillation length i 13 . Moreover, oscillometry is 
a unique tool to probe the existence of oscillations into 



a fourth sterile neutrino. LENA can be considered as 
a versatile tool for neutrino oscillation measurements 
at short baselines. 

5.3.1 Introduction 

The most precise and unambiguous way to detect 
neutrino oscillations is a determination of the oscil- 
lation pattern in the distance-dependent flux of the 
given neutrino flavor over the entire oscillation length. 
Since the oscillation length is proportional to the neu- 
trino energy, neutrino oscillometry would require a de- 
tector hundreds or even thousands of kilometers long 
if used with the present or proposed neutrino beams 
that feature energies of several hundred MeV. As this 
is unrealistic, all beam experiments aiming at neutrino 
oscillations consider just a single or at most a two- 
point measurements instead of the full oscillometric 
approach. Also when using reactor neutrinos, the dis- 
tance from the source to the first minimum is about 
2 km — still beyond the technological and financial 
constraints for a detector. To be able to perform neu- 
trino oscillometry using a realistic-size detector like 
LENA (100m long), one needs a strong source of mo- 
noenergetic neutrinos with an energy of a few hundred 
keV. Comparable sources have already been produced 
by neutron irradiation in nuclear reactors [262 264 . 
Neutrino oscillometry potentially provides a compet- 
itive and considerably less expensive alternative to 
long-baseline neutrino beams. 

5.3.2 Detection principle 

In liquid scintillator, v e s at sub-MeV energies are 
detected by the recoil electrons from elastic neutrino- 
electron scattering. Any decrease in the detection rate 
along the detector that exceeds the geometric factor 
will give, for the first time, a continuous (oscillometric) 
measure of flavor disappearance. 

LENA is well suited for an oscillometric measure- 
ment due to its large height h of 100 m, the low de- 
tection threshold of ~200keV, and the considerable 
fiducial volume of ~35 kt in this energy region (com- 



pare Sect. 4.3). Assuming a light yield of 200pe/MeV 
for LENA, the expected position sensitivity is ~25 cm 
at 500 keV electron recoil energy [6] . The energy res- 
olution will be ~10% in this region [265] , 

As the cross-sections for v-e scattering are tiny, a 
very strong neutrino source has to be used to provide 
adequate statistics. Fortunately, there is a variety of 
radionuclei decaying via electron capture (EC). Since 
EC is a two-body process, the emitted electron neu- 
trino is monoenergetic and carries most of the transi- 
tion energy. Tab. |XV| lists EC isotopes featuring suit- 
able Q-values to produce monoenergetic neutrinos of 
a few hundreds of keV and with half-lives of a few 
months allowing for convenient handling. 

Sources of this kind have been produced in the past 
by neutron irradiation at nuclear reactors: Usually, 
a lighter (A-l), stable isotope of the same element 
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is exposed to the intense neutron flux generated in- 
side a reactor. For the calibration of the GALLEX 
experiment |262j . a 51 Cr source of an initial activ- 
ity of 62PBq (1.7 MCi) was produced by plac- 
ing 36 kg of metallic chromium, enriched in 50 Cr, at 
the core of the Siloe reactor in Grenoble (35 MW ther- 
mal power) for a period of 23.8 days. In principle, 
the 36 kg batch of enriched chromium is still available 
and could be reused for LENA. Assuming an activ- 
ity of 5 MCi, about 1.9 xlO 5 neutrino events would be 
expected. This already enormous statistics could be 
further increased either by repeating the cycle of neu- 
tron activation and measurement runs in LENA or by 
a further increase in source activity. 



5.3.3 Short baseline neutrino oscillations 



The survival probability of electron neutrinos in a 
short baseline experiment can be approximated as 



Pr, 



1 - sin 2 2% • sin 2 (7r • l/Ly), (5.1) 



as long as the mixing effects are clearly disentangled 
due to different oscillation length Ly. The length Lij 
can be written to 



= 2.48 m • 



EL eV J 



McV Am 



(5.2) 



In the three-flavor scenario, the short baseline is L 2 3 ps 
Li 3 . Assuming the value Am§ 2 = 2.5 x 10 -3 eV 2 for 
the mass squared difference that can be derived from 
the global oscillation analysis, the baseline (when ex- 
pressed in meters) is approximately equal to the neu- 
trino energy in keV: 



L 13 [m] « £„[keV], 



(5.3) 



about 3 % of the solar oscillation length L 12 . However, 
if a four-flavor scenario including the sterile neutrino 
of the Reactor Antineutrino Anomaly (RAA) is con- 
sidered, the shortest baseline is due to Am 2 x > 1.5 eV 2 
268 . In this case, the oscillation length L 14 for a EC 
source experiment is ~lm. The number of events in 
an differential volume d^ in the cylinder can be writ- 
ten in the following form [264] : 

dN(£) = -^n e a{Eu)p{E v ,l^m 2 26 ij )dV(e), (5.4) 

where N v is the v e source intensity, £ is the distance of 
the detection region from the source, n e is the detector 
electron density (n e — 3 x 10 29 m~ 3 for an LAB-based 
scintillator), and o{E v ) stands for the v e -e scattering 
cross-section at the neutrino energy E v . The detection 
probability p is a function of the distance to the source, 
the mixing angle Oij and the oscillation length Lij-. 

p(£) = 1 - X (Eu) sin 2 2% ■ sin 2 (^/L lJ ), (5.5) 



where x(E u ) takes into account the effect of the other 
flavors. For sterile neutrinos, x(E^) — 1- The integral 
number of events Ni nt can be deduced from Eq. (5.4). 
It can be presented in the form: 



N tilt = N (l-g(L ij ,h) 



sin 2 26»., 



(5.6) 



where Nq is the expected event number without oscil- 
lations, while the fraction of "disappearing" neutrino 
events is a function of the oscillation probability and 
the geometric factor g that depends on the fraction of 
the oscillation length contained inside the detector of 
can be used as base for a sensitivity 
etection of neutrino mixing. 



height Eq. (5.6 



estimate for the 



5.3.4 Physics case for oscillometry 

Neutrino oscillometry offers an elegant way to ad- 
dress a number of questions related to neutrino oscil- 
lations: a precise determination of the mixing angle 
#i3 and of the oscillation length L13, confirming the 
results of the global analysis. Beyond the standard 
oscillation picture, oscillometry will be very sensitive 
to z/ e oscillations into sterile neutrinos on the eV mass 
scale predicted by the RAA: Based on the observed 
rate deficit in reactor and radiochemical neutrino ex- 
periments, the existence of this fourth neutrino has 
been recently proposed in [268J . An oscillometric mea- 
surement in LENA will allow a precise determination 
of the associated mixing parameters 6*14 and L14. 

Mixing parameters 6*13 and L13. For a precise de- 
termination of #13, an advantage of the short baseline 
oscillometry is the absence of matter effects. These 
effects cause a degeneracy in the determination of the 
oscillation parameters in long-baseline beam experi- 
ments (Sect. 6.2). 

The oscillometric approach to determining #13 is 
a measurement of the differential rate dN/d£ of v-e 
scattering events as a function of the distance £ from 
the neutrino source. Due to the still relatively large 
oscillation length of several hundred meters, such a 
measurement would require a strong EC source and 



multiple measurements: Fig. 25 shows AN/d£ for the 



case of a 5x55 days measurement campaign based on 
a 51 Cr source of an initial activity of 5 MCi. A detec- 



tion threshold of 200 keV is assumed (Sect. 5.3.5 ). The 
rates have been normalized to the full solid angle. The 
dashed lines indicate the statistical ler uncertainties on 
the differential rate, assuming a bin width of 10 m: For 
large I, these uncertainties increase substantially due 
to the geometric decrease in the detected event rate 
with £ 2 . 

Alternatively, the mixing angle can be determined 
by the integral number of events via Eq. ( 5.6 ) : Monte- 



Carlo calculations return N — 3.8x 10 4 /MCi and 
g = 0.65% for 51 Cr (55 days), or N = 2.9xl0 4 /MCi 
and g — 2.3% for 75 Se (160 days). The achieved sen- 
sitivity is only a function of source strength and the 
number of measurement runs: Fig. [26^ shows 90 % ex- 
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Nuclide Ti/2 [d] Q e [kcV] E v [kcV] (BR) £^ max [kcV] Material v intensity [Bq] 



aY Ar 


35 


814 


811 (100%) 


617 


4U Ca, Ar 


8.3xl0 lb 


5i Cr 


28 


753 


747 (90%) 


560 


so Cr 


2.3xl0 16 


75 Se 


120 


863 


450 (96%) 


287 


Se 


l.lxlO 14 


113 Sn 


116 


1037 


617 (98%) 


436 


Sn 


8xlO n 


145 Sm 


340 


616 


510 (91%) 


340 


Sm 


2xl0 12 


i69 Yb 


32 


910 


470 (83%) 


304 


Yb 


l.lxlO 15 



TABLE XV: Potential EC v & sources that can be produced by neutron irradiation in nuclear reactors. The half-life Tim, 
the Q-value of the reaction, the energy E v of the neutrino line and the corresponding branching ratio BR, as well as the 
maximum electron recoil energy -B e ,max are shown. The achievable neutrino source intensities have been estimated for 
1 kg batches of the irradiated elements, assuming natural isotope abundances and a 10-day irradiation with a neutron 
flux of 5x 10 14 n/cm 2 s. Neutron-capture cross sections were taken from [267] . 




FIG. 25: Differential v-e scattering event rate for a measurement 
campaign based on a Cr source installed on top of LENA (5 x 55 d, 
5MCi). The dashed lines indicate the statistical uncertainties (Iff) 
assuming a bin width of 10 m. A correction for the solid angle has 
been applied, resulting in an increase of uncertainties with distance. 
The assumed oscillation amplitude is sin 2 2#i3 = 0.17. 



FIG. 



26: Upper limits for sin 2 20 13 (90% GL.) 
as a function of the initial source strength 
and measurement repetitions. Results for both 
51 Cr and 75 Se are shown, considering statistical 
uncertainties only. 



elusion limits for oscillations for both isotopes. 75 Se 
reaches by far better limits due to the lower neutrino 
energy and therefore larger value of g. In this way, 
sensitivity to sin 2 2#i3 sa 0.1 could be reached by five 
runs with a 3.5 MCi Se-source. However, the required 
exposure will further increase when uncertainties in- 
troduced by the subtraction or suppression of back- 
grounds are considered (Sect. |5.3~5 ). 

In principle, also the oscillation length L13 could 
be determined by oscillometry, although — even for 
75 Se — the spatial oscillation pattern is only partially 
contained within the extensions of LENA. The result 
could be compared with the neutrino energy that is 
usually well known or that can be measured indepen- 
dently very precisely [269] . For 51 Cr, the neutrino en- 
ergy is presently known with a precision of 0.03 %. 
Since Eq. (5.3| is valid if the global-analysis value of 
Am| 2 = 2.5 x 10 -3 eV 2 is used, this comparison will 
be helpful for assessment of the global analysis itself. 



Sterile neutrinos, B44 and L 14 . The LENA detec- 
tor provides unique sensitivity for the possible fourth 



(sterile) neutrino that is introduced by the RAA. Since 
the new neutrino v s is sterile, its existence will man- 
ifest in a disappearance of v e s into i/ s , the ampli- 
tude being governed by the mixing angle #14. The v e 



survival probability P ee (£, L14) is given by Eq. (5.1) 



The best fit values for the mixing parameters are 



sin 2 20 14 = 0.16 and A 



vn. 



> 



1.5 eV 2 [H] (see also 



According to Eq. (5.2), the oscillation length L 



14 

should be rather short, L14 < 1.24 m for the case 
of 51 Cr. Therefore, the oscillation v e -H> v s could 
be observed several times within the first 10 m from 
the source. This opens an excellent possibility for 
direct oscillometry. It is worthwhile to note here 
that oscillation lengths for active and sterile neutri- 
nos, L13 = 742m and L 14 = 1.24m (both for 51 Cr) are 
fully disentangled and can be derived independently. 

The differential event number dN/d£ as obtained 
from Eq. ( 5.4 ) is depicted in Fig. 27 assuming the best- 



fit RAA mixing parameters and a single 55-days run 
with a 51 Cr source. Statistical uncertainties for a bin 
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width of 1 m are far smaller than the disappearance 
amplitude. 

Like for #13, we determine the sensitivity to the 
amplitude sin 2 29^ by the integral event number in 



LENA, using Eq. (5.6 1. As the oscillation is fully con- 



tained within the detector, g reaches the maximum 
value of 50%. The resulting sensitivity as a function 



of source strength and runs is shown in Fig.|28]for 51 Cr: 



The high sensitivity of LENA is clearly demonstrated: 
a single run with a 5 MCi source would be sufficient to 
exclude the best-fit value of the RAA. Similar results 
are expected for 75 Se. 

The RAA analysis only gives a lower limit for 
Amy > 1.5 eV 2 , and therefore an upper limit to the 
oscillation length L 14 < 1.24m (for 51 Cr). As long 
as Lu is large compared to the spatial resolution of 
about 25 cm, a precise determination of this parameter 
can be expected. Therefore, the sensitivity for Lu will 
vanish for Am 2 4 w 7.5 eV 2 . Due to the lower energy 
and therefore lower spatial resolution of 75 Se-i^ recoils, 
51 Cr seems the better candidate for this search. 



5.3.5 Experimental uncertainties 



While Sect. 5.3.4 describes the optimum results for 



oscillometry achievable in LENA, an actual experi- 
ment will suffer from a number of uncertainties re- 
ducing the sensitivity. In the following, we discuss 
two main aspects, the uncertainty of the initial source 
strength and the subtraction of background events 
from the signal rate. 

Source activity. In the calibration campaigns of 
GALLEX/GNO and SAGE experiments that used 
strong sources based on 51 Cr (—2 MCi) and 37 Ar 
(0.4 MCi), great care was given to an exact determi- 
nation of the source activity |262|, 1263) . Various meth- 
ods were used, ranging from precision measurements 
of source weight and heat emission to direct count- 
ing of decays in aliquots of the sources. The greatest 
accuracy reached for 51 Cr was 0.9% [262], 0.4% in 
case of 37 Ar 263 . While this uncertainty will not 
play a dominant role for the detection of large oscil- 
lation amplitudes as in the case of #14, it has consid- 
erable influence if the expected effect is of the same 
order of magnitude, i.e. for # 13 : However, since L 13 
is known from global analysis, a precise measurement 
of the neutrino rate in the first 10— 20 m of LENA can 
be exploited as a normalization for the search in the 
remaining volume, provided the number of events in 
the near volume exceeds ~10 4 . 

Background subtraction. Independent of the used 
isotope, solar neutrinos pose an irreducible back- 
ground for all oscillometric measurements. 7 Be-i/s will 
be detected at a rate of ~0.5 counts per day and ton, 
featuring a maximum recoil energy of 665 keV, only 
slightly above the spectral maximum of 51 Cr. In ad- 
dition, radioactive impurities inside the scintillation 
volume have to be considered: 14 C sets the energy 
threshold of detection to —200 keV, while traces of the 



isotopes 85 Kr, 210 Po and 210 Bi dissolved in the scin- 
tillator will cause background contributions over the 
whole energy range of the source signals [201] . 

In the 7 Be analysis of Borexino |201j . most of the 
background contributions are eliminated by pulse- 
shape discrimination and a spectral fit to the signal 
region, separating the neutrino recoil shoulder from 
background spectra. A similar technique could be 
applied in LENA, the efficiency depending on the 
achieved photoelectron yield and pulse shaping prop- 
erties. This analysis will be aided by the fact that 
the EC source can be removed from the detector, pro- 
viding an exact measurement of the background rates. 
Nevertheless, the £ 2 decrease of the signal rate will 
mean that the background rate will dominate in the 
far-region of the detector, considerably enhancing the 
signal rate uncertainties. This does not assail the 
search for sterile neutrinos that mainly concentrates 
on the parts of the detector closest to the source. How- 
ever, it is a serious issue for ^-experiments in which 
the oscillation signature is limited to the far-region. 
The feasibility of an oscillometric search will depend 
on the availability of strong sources, background con- 
ditions and the efficiency of spectral separation. 

5.3.6 Conclusions 

Thanks to its low energy detection threshold 
(-200 keV) and considerable length (-100 m), LENA 
is exceptionally well suited to perform determina- 
tion of neutrino oscillation parameters. The needed 
electron-capture source emitting high-intensity mo- 
noenergetic and low-energy neutrinos can be produced 
by neutron irradiation in a nuclear reactor: currently, 
MCi-sources of 51 Cr and 75 Se seem the most promising 
candidates. The disappearance of electron neutrinos 
can be monitored over the full length of the detector 
by the neutrino-electron scattering event rate. How- 
ever, radioactive background and the signal of solar 
7 Be neutrinos will reduce the accuracy in the far re- 
gion of the detector. The resulting oscillometric curve 
as well as the integral event number potentially allow 
for an accurate determination of the mixing angles #13 
and 6*i4 as well as the associated oscillation lengths 
£13 and L14. 

LENA will achieve great sensitivity for the sterile 
neutrinos predicted by the RAA [268 : The best-fit 
mixing parameters could be conclusively tested by a 
single run with a 5-MCi 51 Cr-source. Also L14 can be 
determined precisely, provided Am 2 4 is not too large. 
However, a search for oscillations driven by #13 will 
be far more demanding: Multiple runs with strong 
sources as well as excellent detector performance and 
background conditions would be required to reach a 
sensitivity in sin 2 2#i3 that is competitive to current 
reactor v e and long-baseline experiments. 

Although not discussed here in detail, a similar ex- 
periment can be conducted with a strong antineutrino 
source, using e.g. 90 Sr( 90 Y) source with a spectral end- 
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FIG. 27: Differential v-e scattering event rate for a 55-day run 
with a 5MCi 51 Cr -source installed on top of LENA, including a 
correction for the solid angle. The first 10 m are shown. The dashed 
lines indicate the statistical uncertainties (la) assuming a bin width 
of 0.1m. RAA best-fit mixing parameters are used. 
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FIG. 28: Upper limits for sin 2 20 u (90% GL.) 
as a function of the initial source strength and 
measurement repetitions. Results for 51 Cr are 
shown, considering statistical uncertainties only. 



point of 2.3 MeV. While the range of v e energies will 
be limited to the range from 1.8 to 2.3 MeV due to the 
threshold of the inverse beta decay, this will require an 
analysis resolved both in space and energy to compen- 
sate the slight differences in oscillation length. How- 
ever, exploiting the final-state coincidence will corre- 
spond to a strong reduction of background. This will 
allow to perform the search for sterile neutrinos also 
in the v sector |271j . For #13 and L13, the oscillation 
baseline will be too long to show significant disappear- 
ance within the detector. 



5.4 Pion decay at-rest experiment 

The DAEDALUS (Decay At-rest Experiment for 6 C p 
studies At the Laboratory for Underground Science) 
concept 272 proposes a neutrino oscillation experi- 
ment on three comparatively short baselines of 1.5, 8 
and 20 km. Neutrinos are created by charged pions 
decaying at rest, which are in turn produced by high- 
power synchrotrons. Via the decays 
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monoenergetic Vn as well as spectra of v e and are 



generated (Fig. 29). The energies range to a kine- 
matic maximum of ~50 MeV, matching the relatively 
short oscillation baselines. The neutrinos propagate 
from three locations at different distances to a single 
large- volume detector. The sought-for signal is the ap- 
pearance of v e from — > D e oscillations driven by the 
small mixing angle f5 13 . Most importantly, the rela- 
tive rates observed for the medium and far baselines 
depend on the size of the CP-violating phase r5cp- 
The original proposal foresees the LBNE Water 
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FIG. 29: Energy distribution of neutrinos in a n decay 
at-rest beam [2721. 



Cherenkov detector(s) for the v e appearance measure- 
ment. However, LENA features an intrinsic capabil- 
ity for the identification of inverse beta decay events, 
offering excellent background discrimination for this 
channel. The necessary cross-calibration of the neu- 
trino intensities from the three 7r-decay sources will be 
possible via neutrino-electron scattering and charged- 
current reactions on 12 C. Therefore, a combination of 
the DAErSALUS and LENA programs seems a promis- 
ing alternative (or extension) to the discussed long- 
baseline neutrino beam scenarios (Sect. 6.2 1. However, 



detailed calculations on the event and background 
rates as well as the expected sensitivity are needed. 



5.5 Neutrinoless double-beta decay 

The huge amount of instrumented mass provided by 
LENA might open the possibility of a large neutrino- 
less double-beta (0^2/3) decay experiment, based on 
136 Xe dissolved in the liquid scintillator. The solu- 
bility in organic liquid scintillators of 136 Xe at atmo- 
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spheric pressure is about 2% in weight, allowing po- 
tentially an experiment with 200 tons of active mass 
or more. The energy resolution is a very crucial pa- 
rameter in 0^2/3 experiments. Therefore, a more dense 
PMT coverage might be required, at least in the cen- 
tral region of the detector. Further studies are needed 
to asses the real sensitivity of such an experiment. 
Nevertheless, this seems the only realistic way for 0^2/3 
experiments at the 100 1 scale, which arguably would 
be able to attack the normal hierarchy region of neu- 
trino masses. 

5.6 Search for other rare processes 

The unique expected characteristics of the detector 
such as ultra low background and large target mass 



provide the possibility of search for others low energy 
rare processes |273) . They are searches for nucleon de- 
cay into invisible channels N — > 2v and NN — > 2v 
|274j . solar neutrino decay v H — > v L + 7 |275j . heavy 
neutrino emission in 8 B-beta decay [276] and high en- 
ergy solar axions [2771 [278 . Pauli Exclusion Principle 
in nuclei 279 and electron stability [280J can be tested 
at more high level. The precise measurements of solar 
neutrino-electron elastic scattering and excitation of 
nuclear levels by neutral currents allow to study non 
standard neutrino interactions 128144283] . 
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6 GeV physics 



While the emphasis of the LENA physics program is 
on low-energy neutrinos (E < 100 MeV), the experi- 
ment can also contribute to several aspects of neutrino 
and particle physics associated to GeV energies. Ac- 
tually, the search for proton decay into kaon and an- 
tineutrino was one of the first items considered to play 
an integral part in the LENA concept, since the visi- 
bility of the kaon's energy deposition in the scintillator 
highly increases the detection efficiency substantially 
in comparison to water Cherenkov detectors (Sec. 6.1 1. 

In the last years, it also became evident that liquid- 
scintillator detectors will be a serious option for the 
use as a far detector in a long-baseline neutrino beam 
experiment, and for the investigation of atmospheric 
neutrino oscillations. An overview of possible neutrino 
beam experiments and the information that could 
be won from atmospheric neutrinos are outlined in 
Sees. |6.2| and |6.3[ respectively. 



the SU(5) model predict a longer proton-decay life- 
time with a larger uncertainty, typically from 10 30 to 
years 



6.1 Nucleon decay search 

Due to the large target mass and the intended long 
measurement time, LENA offers the opportunity to 
search for nucleon decays. Currently, the best lim- 
its on proton lifetime are hold by Super-Kamiokande 
HH1I1H5|, and it seems not likely that LENA will sub- 
stantially improve the limit for p — > 7r°e + . However, 
the sensitivity for the decay mode p — > K + v is an 
order of magnitude larger than in water Cherenkov 
detectors. Moreover, the search in LENA is expected 
to be background-free for about 10 years, allowing to 
set a limit of t p > 4x 10 34 yrs (90 % C.L.) if no event is 
observed after this period 286 . This already probes 
a significant fraction of the proton lifetime range pre- 
dicted by SUSY theories [257J [555] ■ 

6.1.1 Theoretical predictions 

In the standard model of particle physics, protons 
are stable. This is a consequence of the baryon number 
(B) conservation which has actually been introduced 
empirically into the model. It is interesting to realize 
that there is no fundamental gauge symmetry which 
generates the conservation of B. For this reason, the 
validity of B-conservation can be considered as an ex- 
perimental question. However, several theories beyond 
the standard model actually predict an instability of 
the proton: 

GUT SU(5). In the minimal Grand Unified The- 
ory SU(5), M x ~ 10 15 GeV/c 2 , the predicted life- 
time is T p _ i . Tr Q e + = 10 29 years. The first generation 
of large water-Cherenkov detectors motivated by this 
prediction observed no evidence of proton decay in the 
p — > 7r°e + mode and therefore ruled out the model. 
The lifetime of the proton largely depends on the mass 
scale of the super-heavy particles mediating the de- 
cay process (X and Y bosons). Further extensions of 



10 36 



GUT SO(10). The proton lifetime predicted 
by the SO(10) extension of the SU(5) model is 
around 10 32±1 yrs for non-supersymmetric models and 
10 3 yrs if Supersymmetry is included 



SUSY SU(5). In the minimal supersymmetric SU(5) 
model, the dominant decay modes of the proton in- 
volve pseudo-scalar bosons and anti-leptons [291J: 



K+v, n+V, K°e+, Ky+, 7r°e 



(6.1) 



where the relative strengths depend on the specific 
exchange of the SUSY particles involved. However, 
in most of the models the proton-decay channel p — > 
K+v is favored [2571 HHFJ [Ml ■ The predictions 
concerning the lifetime of the proton are in the order 
of 10 33 to 10 34 years [287ll288j . 

6.1.2 Detection mechanism 

Within the target volume of LENA, about 1.6xl0 34 
protons, both from carbon and hydrogen nuclei, are 
candidates for the decay. This number has been cal- 
culated for PXE (Sect. |2~3[ ), which is assumed as ref- 
erence in the following. As all decay particles must be 
contained inside the active volume, the fiducial volume 
is about 5% smaller. In the case of LAB, the proton 
number will be further reduced by about 6 % due to 
its lower density (see also Tab.[l]). 

In the case of protons from hydrogen nuclei 
(-0.25 xlO 34 protons in the fiducial volume of LENA), 
the proton can be assumed at rest. Therefore, the pro- 
ton decay p —> K+v can be considered as a two-body 
decay problem, where K+ and v always receive the 
same energy. The energy corresponding to the mass 
of the proton, m p — 938.3 MeV is thereby given to 
the decay products. Using relativistic kinematics, it 
can be calculated that the particles receive fixed ki- 
netic energies, the antineutrino 339 MeV and the kaon 
105 MeV. 

The antineutrino escapes without producing any 
detectable signal. However, the large sensitivity of 
LENA for this decay channel arises from the visibil- 
ity of the ionization signal generated by the (kinetic) 
energy deposition of the kaon. A water Cherenkov de- 
tector is blind to this signal as the kaon is produced 
below the Cherenkov threshold in water; only the sec- 
ondary decay particles are visible, greatly reducing the 
detection sensitivity. 

In LENA, the prompt signal of the decelerating kaon 
is followed by the signal arising from the decay parti- 
cle^): After t k + = 12.8 ns, the kaon decays either by 
K+ ->• /x+i^ (63.43%) or by K+ ->• tt+tt (21.13%). 
In 90% of these cases, the kaon decays at rest 293 . 
If so, the second signal is again monoenergetic, either 
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corresponding to the 152 Me V kinetic energy of the 
/i + or 246 MeV from the kinetic energy of the tt + and 
the rest mass of the tt° (which decays into two gamma 
rays creating electromagnetic showers) . A third signal 
arising from the decay of the muon will be observed 
with a large delay (t^+ = 2.2 us). A more detailed 
discussion can be found in [294] . 

If the proton decays inside a carbon nucleus 
(^1.2xl0 34 protons in the fiducial volume), further 
nuclear effects have to be considered. First of all, 
since the protons are bound to the nucleus, their ef- 
fective mass will be reduced by the nuclear binding 
energy Ef,, 37 MeV and 16 MeV for protons in s-state 
and p-states, respectively. Secondly, decay kinemat- 
ics will be altered compared to free protons due to 
the Fermi motion of the proton. The Fermi mo- 
menta in carbon have been measured by electron scat- 
tering on 12 C |295j . The maximum momentum is 
about 250MeV/c. A range for the effective kinetic 
energy of the kaon has been derived by Monte Carlo 
simulations: (25.1 — 198.8) MeV for the s-state and 
(30.0 - 207.2) MeV for the p-state [253] . 

In any case, the experimental signature of the pro- 
ton decay in LENA is not substantially affected by 
nuclear effects or the kaon decay mode: A coincidence 
signal arising from the kinetic energy deposited by the 
kaon and from the delayed energy deposit of its decay 
particles will be observed. 

6.1.3 Background rejection 

The main background source in the energy range of 
the proton decay are atmospheric muon neutrinos v^. 
Via weak charge-current interactions, these create fj, 
inside the detector, a substantial fraction in the energy 
range relevant for the proton decay search. Moreover, 
additional kaons can be produced in deep inelastic 
scattering reactions at higher v^. In the following, the 
arising background rates and possible rejection cuts 
are shortly outlined. For a thorough discussion, see 
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Muon events. The rate of muon events from atmo- 
spheric neutrinos in the relevant energy range can be 
derived from Super-Kamiokande measurements [293 . 
At Pyhasalmi, the rate corresponds to 1190.4 Vn- 
induced muons per year [294] . 

In order to distinguish the real proton decay signals 
from muon background events, a pulse shape analysis 
can be applied. MC simulations show that the kaon 
deposits its energy within 1.2 ns, leading to a fast but 
resolvable coincidence with the kaon decay products 
after t k + (Sec. 6.1.2 1. A typical time profile is shown 
in Fig.[30| This double signature can be used to dis- 
criminate atmospheric Vn events as long as the kaon 
decay is sufficiently delayed to produce a discernible 
double signal, i.e. the delay is large compared to the 
time resolution of the detector. 

In the analysis presented in [294] . signal and back- 
ground events are discriminated via the signal rise 
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FIG. 30: Signature of the proton decay into kaon and an- 
tineutrino in LENA, based on the analog sum of all PMT 
channels. The prompt signal is generated by the decelera- 
tion of the kaon. After a delay of (in this case) 14 ns, the 
second peak is caused by its decay particles [294] . 



time. No time-of-ffight correction is applied. Based 
on 2x10 proton decay and muon events in the rele- 
vant energy regime, an analysis cut can be defined that 
rejects all muons and retains a detection efficiency of 
e p sa 65 % for proton decay. The sensitivity e p is an 
order of magnitude larger than the one obtained in the 
Super-Kamiokande analysis, corresponding to a simi- 
lar increase in the proton lifetime limit. 

The corresponding background rejection efficiency is 
at least e M > 1— 5x 10~ 5 . This results in an upper limit 
of ^0.05 muon events per year that are misidentificd 
as proton decay events. 

The prominence of the double-peak signature gen- 
erated by the kaon and decay particles depends on the 
performance of the liquid scintillator. The two crucial 
parameters are the optical transparency, especially the 
scattering length, and the fast fluorescence time con- 
stant n discussed in Secs. |2.3| and |3.1| Moreover, there 
is also a dependence on the position of the event within 
the detection volume. The studies performed in [294J 
used a fixed value of T\ — 3.5 ns (similar to PXE) and 
varied the light propagation parameters as well as the 
vertex position. Assuming reasonable values for ab- 
sorption and scattering lengths, e v is changing only 
on a level of a few percent. An LAB-based scintilla- 
tor might feature a slightly lower efficiency due to the 
larger t\ — 5.2 ns. A study performed using a pes- 
simistic value of t\ = 6 ns returned an efficiency of 
e p = 58%fM\. 

Hadronic event. In case of charged current reactions 
of atmospheric v^s at larger energies, hadrons can 
be produced along with the final state muon. These 
events are dangerous if they are able to mimic the dou- 
ble signature of the proton decay. While this is not 
the case in pion and hyperon production, interaction 
modes creating an additional kaon in the final state 
may be mistaken as signal events [297] . In principle, 
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these events can be discriminated by the additional 
decay electron of the muon created in the CC reac- 
tion. However, this signal is sometimes covered by the 
muon signal itself: Monte Carlo simulations return an 
upper limit of 0.06 irreducible background events per 
year for this channel. 

6.1.4 Proton decay sensitivity 

Based on the efficiencies of the rise time cut, the 
sensitivity of LENA for the proton decay search can 
be determined. The observed activity due to proton 
decays is given by the expression: 
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where e p — 0.65 is the efficiency, N p — 1.45x10 
the number of protons in the fiducial volume, t m is the 
measurement time and r v is the lifetime of the proton. 

If the proton lifetime corresponds to the current best 
limit to this channel by Super-Kamiokande, (t v = 2.3- 
10 33 yrs) 285 , about 40.7 proton decay events will be 
observed in LENA in a measurement time of ten years. 

For establishing a new lower limit on proton life- 
time, the number of background events observed over 
the measurement time is the main issue. Combin- 
ing the expected background rates from atmospheric 
neutrino-induced muon and kaon production, a rate 
of 0.11 background events per year or 1.1 events in 10 
years can be obtained. This result is an upper limit 
on the expected background rate |294j . 

In case there is no signal observed in LENA within 
these ten years, the lower limit for the lifetime of the 
proton will be placed at r p > 4xl0 34 yrs at 90% C.L. 
If one candidate is detected, the lower limit will be 
reduced to r p > 3xl0 34 yrs (90% C.L.), featuring a 
32 % probability that this event is due to background 

6.1.5 Conclusions 

The MC studies carried out in |294j determine the 
efficiency for the proton decay search in LENA to 
^65%. Based on this, a new lower limit for the pro- 
ton lifetime of t > 4xl0 34 yrs (at 90% C.L.) can be 
reached if no proton decay event is observed within ten 
years. The high efficiency is based on the distinct pulse 
shape of the proton decay mode p — > K + v in LENA. 
Since the values predicted by the favored theories for 
the proton decay in this channel are of the order of 
the value resulting from this analysis |287l [288 , it is 



obvious that LENA measurements would have a deep 
impact on the proton decay research field. 

LENA might also provide relevant sensitivity levels 
to other nucleon decay channels. While the analysis 
presented here is independent of the tracking capabil- 
ities of the detector (Sect. 3.3), in others (e.g. p — > 
7r e + ) the possibility of reconstructing the decay ver- 
tex might be necessary to discriminate background sig- 
nals. However, these aspects require further studies. 



6.2 Long-baseline neutrino beams 

Accelerator-based neutrino beam experiments might 
prove to be the only viable way to determine the value 
of the CP-violating phase S CP in the leptonic sector 
and the neutrino mass hierarchy. The oscillation base- 
lines discussed today range from hundreds to thou- 
sands of km, corresponding to GeV neutrino energies. 
In spite of its focus on low-energy neutrinos, LENA 
might serve as a far detector for such an experiment. 
We review general properties of future beam experi- 
ments and discuss both a conventional neutrino beam 
to Pyhasalmi and a beta-beam to Frejus in the context 
of the tracking capabilities of LENA. 

6.2.1 Concept and goals 

During the last decade we witnessed drastic changes 
in our understanding of neutrinos. A number of ex- 
periments have shown that neutrinos violate lepton 
flavor through oscillations and therefore must have 
mass. But despite the large progress many fundamen- 
tal questions remain unanswered. Some of these ques- 
tions can be answered by sending an artificial neutrino 
beam over a long distance (several 100 km) to LENA. 
This is the subject of this chapter. The questions that 
can be answered with a neutrino-beam are: 



What is the value of the last unknown mixing 

?1 3 ' 



angle 9y 8? 



• What is the hierarchy of the neutrino masses 
(sign of 6*23 )? 

• Is the mixing angle 023 maximal? 

• Do neutrinos violate the CP symmetry? 

A number of other experiments are trying to find the 
value of . These are reactor neutrino experiments 
(DoublcChooz, DayaBay, Reno) and long-baseline ex- 
periments (Minos, T2K). First indications of a value 
between 5 and 10 degrees have been found. The value 
is important to see the pattern of neutrino mixing and 
to understand the sensitivity of most other measure- 
ments. Even if measurements are achieved in the near 
future, LENA will still provide one of the most precise 
measurement of 0i3 and similarly of 023 • To our cur- 
rent knowledge 023 is consistent with maximum mix- 
ing for 45°. Here precision is important to understand 
whether 023 is exactly maximal or 'just' accidentally 
close to 45°. Furthermore we can determine the sign 
of sin 023 and therefore the mass hierarchy. 

But the most important goal from these measure- 
ments is the search for CP-violation in neutrino oscil- 



There are first indications for a non-zero value of $13 from 
the T2K, Minos, DoubleChooz experiments. Even if these 
indications turn out to be valid, LENA will still provide one 
of the most precise measurements of 
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lations. Today we know that two of the three mix- 
ing angles are substantially larger in neutrino mix- 
ing compared to the quark mixing, allowing for much 
larger CP-violation in the lepton sector. The relevant 
quantity - Jarlskog's determinant - is approximately 
4 • 10~ 5 for the quarks and 0.028 sin <5for the leptons, 
if we assume a value of 5° for 613, not far below the 
current limits. The CP- violating phase 5 is unknown 
today. Its measurement will be the prime goal. 

Moreover, matter effects modify the survival proba- 
bilities of v and v for beams over very long distances. 
The sign of this change can in principle be exploited to 
determine the neutrino mass hierarchy. However, the 
combined effects of CP violation and matter might be 
hard to disentangle. 

Artificial neutrino beams are produced from the de- 
cay of certain unstable particles, emitting neutrinos in 
their decay. High intensity beams of the mother parti- 
cle are produced and directed towards LENA, creating 
a more or less collimated neutrino beam. The neutrino 
beam travels through the earth while all other beam 
particles are absorbed. The technologies to produce 
neutrino beams may be split into three classes. 

In so-called conventional neutrino beams a high 
intensity proton beam is directed onto a target. In 
the induced hadronic interactions pions are produced 
which subsequently decay as tt — > [iv^. Only the 
charged pions contribute. Magnetic horns focus pions 
of one charge in the direction of LENA and defocus 
the other charge. By switching the polarity of the 
horns one may choose between neutrinos and 
anti-neutrinos v^. The beam is broad in energy and 
has a small contamination from electron-neutrinos. 



Beta-beams emerge from the decay of radioactive 
ions which are accelerated and stored in a storage 
ring with a straight section pointing towards the 
detector. Accelerating f3 + emitters produces a beam 
of u e s and fi~ emitters D e &. Producing these ions in 
sufficient quantity is a technological challenge. An 
accelerator complex is needed to accelerate the ions 
to high energies. The beam is well focused, broad in 
energy and pure in flavor. 

A neutrino factory produces pions in the technol- 
ogy of the conventional neutrino beam. But now the 
muons from the decay of the pions are captured, re- 
formed into a beam and accelerated. From the de- 
cay of negative muons (/jT —> e~h> e v^) a neutrino 
beam with two flavors is created. Accelerating pos- 
itive muons creates the opposite flavors. The beam is 
higher in energy, somewhat less focused and pure in 
its two flavors. 

The first two technologies are described in more de- 
tails in the following sections. The neutrino factory is 
not pursued further as it needs a magnetized detector 
to distinguish v from v interactions. It is unlikely that 



LENA will be magnetized because of the negative ef- 
fect on the PMT performance. There might be the 
possibility to use recoil neutrons and protons to iden- 
tify v and v. However, further MC studies are needed 
on this aspect. 



6.2.2 Conventional neutrino beam 

For the beam source, the most evident candidate is 
SPS at CERN, producing a beam of 400 GeV protons. 
The current maximum power is 300 kW, a 400 kW up- 
date is intended. Using the planned PS2 with SPS will 
also permit to increase the proton power to 1.2 MW 
(by 2016). In case of other future upgrades in the pro- 
ton production chain, other options may be available. 
To produce neutrinos of 3-5 GeV by pion (and kaon) 
decay we need at least 20 GeV protons. Larger en- 
ergies may produce more flux, but on the other hand 
the high-energy tail may induce more neutral current 
background. 

Neutrinos are produced via decays of pions and 
kaons whose fractions are typically 10 pions for one 
kaon. The decay modes are: 



K+ 



(99.98%) 
(0.01%) 

(63.4%) 
(0.0015%) 



7r°e+i/ e (5%) 
ttV 



(3%) 



Neutrinos are also produced by the muon decay, 
though this is mainly background: 



e + v e v b 



Independent of the initial proton energy, the tt spec- 
tra are always peaked at just below 500 MeV (Fig. 31 ). 
The flux depends quite linearly on both proton energy 
and luminosity. The resulting v spectra feature a max- 
imum energy of E max — 0.43i? w or E mayi = 0.95£j<- , re- 
spectively. The shape of the v spectrum in forward di- 
rection depends also on the magnetic focusing system 
for tt/K's. It determines the number and spectrum of 
the pions entering the decay pipe. Optimizations and 
simulations for the beam are being made elsewhere 
298, 299J . Fig. 32 depicts some sample spectra from 



those simulations. 

Using the CNGS beam as an example, the beam 
is 97% pure, with small admixtures of v e (1%), v e 
(0.1 %) and (2 %). The would feature equal but 
opposite contamination. The beam will be run in two 
phases, one with positive focus, (CP + ) and other with 
negative focus (CP - ). Asymmetrical running times 
might be used, e.g. 2 yrs (u) and 6 yrs (p) , because of 
the lower cross section of vs. 
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FIG. 31: Simulated pion spectra (flux/proton energy), for 
different proton energies from |300| . To get 4 GeV neu- 
trinos we need at least 10 GeV pions and consequently 15 
GeV protons, though at least 20 GeV proton beam would 
be preferred. 5 GeV proton beam will give neutrinos of 
1-2 GeV and below. 



accelerator 




SPS PS2 




old 


upgrade with PS2 


beam enegy [GeV] 


400 


400 400 50 


pot [10 19 /y] 


7.6 


11 33 


E p N pot [10 22 GeVpotyr] 


3 


4.4 13.2 


Beam power [MW] 


0.3 


0.4 1.2 0.4 



TABLE XVI: Assumed properties of potential accelerators 
to be used as the neutrino source: The currently running 
SPS (old), a possible upgrade of SPS, SPS combined with 
PS2 (planned for 2016), and the PS2 on its own. 



6.2.3 Beta-beams 

For beta-beams, the choice of isotopes depends on 
the baseline. Potential ions are listed in Tab. lXVlTl 
6 He and 18 Ne are considered in case of the short base- 
line from CERN to Frejus. In contrast, 8 Li and 8 B 
are the best choice for a large baselines such as CERN 
to Pyhasalmi as both isotopes feature relatively high 
Q-values of about 13 MeV. The maximum neutrino en- 
ergy is given by the relativistic 7 factor at production 
times Q. For example, with the SPS (450 GeV p + 
energy), the baseline for a detector at the first oscil- 
lation maximum is 1100 km and 2100 km for the two 
isotopes. The baseline from CERN to Pyhasalmi is 
2300 km. It is obvious that a beta beam to Pyhasalmi 
would be technically very demanding, and we primar- 
ily consider the CERN-Frejus baseline in the following. 

Fig. [33] shows the layout of a beta-beam facility at 
CERN [299] . Radioactive ions are produced as neutral 
gas, ionized in an ECR source and accelerated. The 
acceleration starts with a LINAC and a rapid cycling 
synchrotron and continues with the existing PS and 



type isotope 


Z 


A A/Z T 1/2 


Qfs 


(K) 


(El ab ) 






[sec] 


[MeV] 


[MeV] 


[GeV] 


(/3+) 8 B 


5 


8 1.6 0.77 


13.9 


7.37 


4.15 


18 Ne 


10 


18 1.8 1.67 


3.4 


1.86 


0.93 


OS") b He 


2 


6 3.0 0.81 


3.5 


1.94 


0.58 


8 Li 


3 


8 2.7 0.84 


13.0 


6.72 


2.27 



TABLE XVII: Potential isotopes for the creation of a beta- 
beam to LENA. (EZ) is measured in the rest frame of the 
decaying isotope. 



SPS machines. Finally, the ions are injected into a 
decay ring with a straight section to the detector. Ions 
are continuously injected into the decay ring which is 
running at fixed energy. The intensity goal is to direct 
in the order of 10 18 neutrinos to the detector per year. 

Even for short baselines, the beta-beam concept 
poses technological challenges. The biggest is the pro- 
duction, collection and ionization of a sufficient num- 
ber of isotopes. A number of different concepts are 
discussed. The ISOL method is considered in the EU- 
RISOL study, production in a ring was proposed in 
|301j . and for some isotopes direct production with a 
deuteron beam is possible [302 . Other challenges are 
the injection and storage of such a large number of ions 
in the decay ring especially for rings with high 7 and 
the collimation of the decay losses in the accelerators. 

6.2.4 Synergies and perspectives 

At the time LENA will start data taking, the value 
of 0i3 might be already established. If it is not too 
small, LENA might stand a good chance to discover 
CP-violation by determining the phase 5cp in the 
PMNS matrix. In a long-baseline oscillation experi- 
ment, (5cp exhibits itself by different oscillation prob- 
abilities for v and v. 

Conventional Beam to Pyhasalmi. In a conven- 
tional Vfx/Vfx beam, V (v^ — > v e ) will be different from 
'P {pn - > Ve)- To discover this difference, the beam is 
operated for a certain amount of time with vs and 
then a matching amount with Ps. For Pyhasalmi, the 
beam will have to travel along a relatively long base- 
line (> 1000 km). In this case, oscillation probabilities 
in the far detector are changed by the Earth matter 
effect: In fact, v e and v e are affected differently due 
to their different interaction cross sections with the 
electrons contained in terrestrial matter. The sign of 
the changes depends on the neutrino mass hierarchy. 
This is a disadvantage in view of a clear determina- 
tion of (5cPj but also offers the opportunity to dis- 
cover both CP violation and mass hierarchy at the 
same time. The possible degeneracy might be resolved 
by using two different baselines or a common analysis 
with other experiments. 

The conventional beam is an appearance experiment 
v„ — > v e at the far detector. As presented in Sect 
LENA features excellent flavor identification and 
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FIG. 32: Recent neutrino spectra from a conventional beam by Longhin. Left column: neutrino mode; right coulumn: 
anti-neutrino mode. The configuration in the upper plots is optimized for long-distances (> 2000 km) and the lower plots 
for medium-distance (1000 km). The flux of each flavour is indicated on the plot in units of neutrinos per 100m 2 and 
year at a distance of 100 km from the target and as a relative fraction of the total flux in percent. 



ter than 5 % energy reconstruction for energies above 
1 GeV. However, backgrounds due to NC ir° produc- 
tion play an important role in determining the sensi- 
tivity for Sep and the mass hierarchy. Further studies 
on the discrimination of this background are necessary. 

Beta-Beam to Frejus. Alternatively if one is 
mainly interested in Sqp or 6*13 turns out to be small, 
a beta-beam over the short distance from CERN to 
Frejus might prove the better option. For a beta- 
beam, 813 is found by the appearance signal of u^s 
at the far detector, and Sep by the comparison of the 
oscillation probabilities V (y e — > v^) and V (D e —> v^). 
Therefore, it is necessary to reliably isolate the weak 
Vn signal from the large number of v e events. The re- 
construction studies presented in Sect. 3.3 indicate an 
excellent rejection efficiency for quasi-elastic v e events 
and a reliable vertex reconstruction if the energy 
of the final state muon exceeds 200 MeV. However, 



the discrimination efficiencies for NC/CC backgrounds 
producing a charged pion in the final state has not 
been evaluated yet. 

In case of a combination of conventional and beta 
beam, the availability of beams of v e /v e and v^v^ al- 
low to test even more fundamental symmetries: T- or 
even CPT-invariance. T-invariance can be tested by 
comparing for example V (y^ — > v e ) with V (y e — > v^). 
If CPT is conserved, a discovery of CP-violation must 
be accompanied by break-down of T-invariance. The 
CPT-symmetry can be tested in the comparison of 
V (y^ — > v e ) with V (y e — > v^). It might be reason- 
able to start with either conventional or beta-beam to 
search for CP-violation, and to add the second beam 
if CP-violation is observed. 

Extensive studies on the discovery range of CP- 
violation with LENA are still missing. A few beam 
configurations have been simulated with GLoBES 
|3TJ31 13D2. See for example [HE HI]. The existing 
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EURISOL scenario 




FIG. 33: Conceptual layout of a beta-beam facility at CERN from the EURISOL design study. 



simulations show that for not too small values of #13 a 
substantial fraction of the parameter space (typically 
60 to 80 %) of the CP-violating phase 8 can be covered. 
Typically the sensitivity is good for values of sin 2 2#i3 
down to 0.01 and then starts to diminish below. 



6.3 Atmospheric neutrinos 

Based on the current status of the tracking studies pre- 



sented in Sect. 3.3 LENA also offers the opportunity 
to investigate atmospheric neutrinos. The large vol- 
ume substantially extends the sensitivity of previous 
large-volume liquid-scintillator detectors to the multi- 
GeV region, filling the energy-gap between previous 
underground experiments and high energy neutrino 
telescopes which become sensitive above 10 GeV. 

Compared to Water-Cherenkov detectors a good en- 
ergy resolution is expected up to 20 GeV and even 
above. We therefore expect a good measurement of 
the flux and angular spectrum of atmospheric neutri- 
nos up to a few tens of GeV. 

These measurements will depend on the ability to 
identify the neutrino flavor in charged current interac- 
tions and to identify simultaneously the neutral cur- 
rent interaction rate for the total flux normalization. 
The separation of neutral currents from is e is based 
on the separation of 7r° from electrons, while the dis- 
crimination of charged pions from muons must rely on 
the identification of the it decay. The performance 
of LENA regarding these issues is still under study 
(Sect.|3~3|. 

The phenomenology of atmospheric neutrinos is rich 
in the multi-GeV region in particular with respect to 
oscillations |307|, 1308] . The measurement of the direc- 
tion and energy of the incoming neutrino will estab- 
lish a long baseline experiment with variable baselines 



from a few tens of kilometers for neutrinos from above 
to more than 12 000 km for vertically up-going neutri- 
nos originating from air showers on the other side of 
the Earth. 

For vertical muon neutrinos the survival proba- 
bility of oscillates with a broad minimum at about 
20-25 GeV and increases back to unity above that 
|307j . Towards the horizon this 1 st oscillation shifts 
to smaller energy down to 1 GeV and the higher or- 
der minima at lower energy do likewise. The good 
energy resolution and statistics of LENA may allow 
to measure the zebra-shaped patterns of alternating 
oscillation minima and maxima with unprecedented 
resolution. This will allow to measure #23 & n d Am 2 ^ 
to a high precision and hence, to probe the oscilla- 
tion hypothesis in a previously not tested parameter 
region. 

Correlated to the disappearance of muon neutrinos, 
we expect the appearance of tau neutrinos. Again, 
this would provide a unique tool to study the parame- 
ter space and verify the oscillations of neutrinos. The 
ability to utilize the v T detection channel strongly de- 
pends on the ability of LENA to identify and separate 
tau neutrinos, which do, however, appear at a sub- 
stantially higher rate than e.g. backgrounds from at- 
mospheric v e or neutral current interactions: Above 
a few GeV, the ratio i? Me of to v e fluxes increases 
from il„ e 2 at 1 GeV to R^ e s» 5 — 10 at a few tens 
of GeV. 

Very interesting structures in the zenith dependent 
oscillation probabilities appear, if also matter oscil- 
lations are taken into account |308l I309| . For nadir 
angles 8 < 33° the neutrinos have travelled through 
the core of the Earth and a strong resonance pattern 
appears, e.g. with maximum disappearance of v e at 
about 3 GeV. For neutrinos not crossing the core the 
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FIG. 34: Flavor conversion probabilities for neutrinos (left panel) and antineutrinos (right panel) for cos (9 = —0.7 versus 
the (anti) neutrino energy. Both the v e — > v a (p e — > u a ) and — > v a (p^ — > u a ) probabilities are plotted. We use normal 
hierarchy of v masses and a vanishing CP phase along with the best-fit oscillation parameters (see text for details). 



effect abruptly shifts to larger energies. The ampli- 
tude of the above structures depend on (9 13 and differs 
for neutrinos and anti-neutrinos. The dependence on 
the CP phase and S is small. A change in hierarchy 
switches the resonance pattern from neutrinos to an- 
tineutrinos resulting in different event rates between 
the normal and inverted hierarchies. The difference 
arises because of a smaller flux and cross-section for 
antineutrinos than neutrinos [307] . A complete study 
of the performance of LENA regarding neutrino mass 
hierarchy measurements is under development. 

Fig. 34 extracted from [310J, depicts the conversion 
probabilities for neutrinos (left panel) and antineutri- 
nos (right panel). The Preliminary Reference Earth 
Model |248] is used for the density profile inside the 
Earth. The assumed input v mixing parameters are: 
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2.4 x 10~ 3 eV 2 , 

8 x 10~ 5 eV 2 , 

0.31, 
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sider normal v mass hierarchy only and cos 9 - 
The upper and lower plots correspond to v e 



(v e — » v a ) and 



fa -> v a ) 



The CP violating phase is set to 8 = 0. We con- 



: -0.7. 
-> v a 
conversions, re- 
spectively. The anti-neutrino conversion probabilities 
are not affected by matter in case of normal v mass hi- 
erarchy and vacuum conversion formalism apply. For 
this value of the nadir angle, i/'s do not pass through 
the Earth's core. Conversions mostly take place in the 
mantle with an average density of (p) ~ 5gcm~ 3 . The 
dip at ~6GeV for P^-^ in Fig.[34| (left panel, upper 
plot) corresponds to the 1-3 or high MSW resonance 
energy E H = Am 2 3 cos 20 13 /(2v / 2G F (p)) « 6 GeV. 
The width of the dip is 2tan29 1 ^E H rj 3.6 GeV. At 
energies ^> Eh, the conversion probablities are domi- 
nated by vacuum oscillation. A complete study of the 
performance of LENA regarding neutrino mass hier- 
archy measurements is under development. 



A high statistics measurement with good energy- 
and angular resolution and flavor identification as it is 
anticipated by LENA has the opportunity to use at- 
mospheric neutrinos as a new tool for science ranging 
from precision neutrino physics to an improved under- 
standing of the Earth's interior. 
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7 Conclusions 

Liquid scintillator is a very attractive detection tar- 
get for the next generation of large-volume neutrino 
observatories: 

• Availability. The organic liquids that serve as 
the primary materials of the scintillator arc used 
in very large quantities in industry. Therefore, 
they are both economically produced and easily 
available. Due to the large market, their indus- 
trial re-use is easy. 

• Performance. The light yield of organic liquid 
scintillators is roughly 50 times higher than the 
light emission in water by the Cherenkov effect. 
It was demonstrated in Borexino that extremely 
high levels of radiopurity can be reached in liquid 
scintillators. This provides the opportunity to 
search for very rare events at the energy scale of 
natural radioactivity (down to 200 keV), in par- 
ticular for low energy neutrino astronomy and 
neutrino geology. 

• Detection channels. Besides electrons and 
protons, organic liquids offer 12 C (and 13 C) as 
target material for neutrino interactions. The 
multitude of interaction channels is the base for 
a spectroscopic measurement differentiating be- 
tween the flavors of neutrinos and anti-neutrinos. 
This might be crucial for investigating the com- 
plex neutrino signature of a Supernova explo- 
sion. In such an event, neutrino- matter as well 
as neutrino-neutrino interactions might give rise 
to extensive swaps and distortions of the initial 
flavor spectra. 

• GeV tracking. Recent investigations of the 
photon arrival times for GeV particles in a ho- 
mogenous large-volume scintillation detector in- 
dicate an unexpected accuracy of directional 
information and particle identification. These 
parameters are very well determined for track 
lengths that exceed several tens of centimeters, 
corresponding to particle energies of several hun- 
dred MeV. 

• Versatility. The detection sensitivity in a large- 
volume scintillation detector covers an energy 
range reaching from sub-MeV energies to the 
scale of several GeV, providing access to a large 
range of topics in neutrino physics, geology and 
astronomy. 



In this work, we have presented a broad variety of sci- 
entific research areas that can be addressed by a liquid- 
scintillator detector. In the context of a European 
Large Infrastructure for astro particle physics, LENA 
is put forward as a viable and cost effective alternative 
for a next-generation neutrino detector. The fields of 
research enclose low energy neutrino astronomy as well 
as elementary particle physics, which can be accessed 
by the investigation of natural neutrino sources and 
includes also nuclcon decay search. If LENA is used 
as far detector in a next-generation neutrino beam ex- 
periment, this will allow for a unique investigation of 
neutrino oscillation parameters as well as CP-violation 
in the lepton sector. 

Profound expertise has been obtained in construc- 
tion and operation of the presently running liquid- 
scintillator detectors KamLAND and Borexino. Their 
successes in neutrino physics and astronomy reflect the 
technological maturity. The results of the LAGUNA 
design study as well as of the specific design stud- 
ies investigating a possible realization of LENA in the 
Finnish Pyhasalmi mine indicate a time frame of 8 
to 10 years for an executive design and the detector 
construction. 
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